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We review fundamental issues arising in the exact solution of the one-dimensional Hubbard 
model. We perform a careful analysis of the Lieb-Wu equations, paying particular attention to 
so-called 'string solutions'. Two kinds of string solutions occur: A strings, related to spin degrees 
of freedom and k-A strings, describing spinless bound states of electrons. Whereas A strings were 
thoroughly studied in the literature, less is known about fc-A strings. We carry out a thorough 
analytical and numerical analysis of fc-A strings. We further review two different approaches to the 
thermodynamics of the Hubbard model, the Yang- Yang approach and the quantum transfer matrix 
approach, respectively. The Yang- Yang approach is based on strings, the quantum transfer matrix 
approach is not. We compare the results of both methods and show that they agree. Finally, we 
obtain the dispersion curves of all elementary excitations at zero magnetic field for the less than 
half-filled band by considering the zero temperature limit of the Yang- Yang approach. 

PACS numbers: 71.10.Fd, 71.10.Pm, 71.27. +a 

Key words: Hubbard model, strongly correlated electrons, thermodynamics, excitations 



Contents 



Introduction 



II 



Bethe ansatz for the Hubbard model 



A Eigenfunctions and eigenvalues 



B SO (4) symmetry 



6 

6 

8 

9 

D Completeness of the Bethe ansatz| 10 



C Discrete Takahashi equations 



III 


Lieb-Wu equations for a single down spin (I) Graphical solution 






A All charge momenta k„ real 






B fc-A two string 






C Summary . . . 





11 

12 



[IV Lieb-Wu equations for a single down spin (II) Self-consistent solution 






A Zcroth order and the discrete Takahashi equations 




B First order corrections 




C Summary 





18 



20 



V Lieb-Wu equations for a single down spin (III) Numerical solution 



22 

A Numerical mcthodl 22 



B 



C 



Numerical solution for two electrons 23 

1 Equations for the fc-A string 23 

2 Numerical solutions for N = 2 23 



N umerical solution for three electrons 27 

Equations for the fc-A string| 27 



Numerical examples of fc-A string solutions . . . 27 

A complete list of solutions for N = 3 and M = 1 29 



* Address after March 1, 1999: Department of Physics, Ochanomizu University, 2-1-1 Ohtsuka, Bunkyo-ku, Tokyo 112-8610, 
Japan 



1 



D Summary 



VI Thermodynamics in the Yang- Yang approach and excitation spectrum in the infinite volume| 



A Takahashi's Thermodynamic Equations 



B Thermodynamics 



C Zero Temperature Limit 



D Ground state for a less than half-filled bandj 



E Excitations for a less than half-filled band 



F Excitations in the half-filled bancj 



G Relation to the root-density formalism 



H Summary 



VII Thermodynamics in the quantum transfer approach 



A The classical counterpart] 



B Diagonalization of the Quantum Transfer Matrix 



C Non-linear integral equations 



D Analytical solutions of the integral equations 



1 Strong-coupling limit 



2 Free-Fermion limit 



E Numerical Results 



3 Low-temperature asymptotics 



4 High-temperature limit] 



F Equivalence of string based thermodynamics with QTM approach 



VIII Conclusions 



APPI NDIXES 



Derivation of the wave function 



General setting 



Two electrons 



Three electrons 



N electrons 



Inhomogeneous Heisenberg model 



1 Algebraic Bethe Ansatz 



2 Explicit expressions for the eigenstates 



The spectrum of three electrons with one-down spin for L — 6 



3 



I. INTRODUCTION 



The Hubbard model was introduced as a simple effective model for the study of correlation effects of d-electrons 
in transition metals (TJQ] (see also ||). It is believed to provide a qualitative description of the magnetic properties 
of these materials and the Mott metal-insulator transition Q . Despite of its appealing conceptual simplicity rigorous 
results for the Hubbard model are rare. The dimension of the underlying lattice is a crucial parameter. Two of the 
most important theorems valid for arbitrary lattice dimension are due to Nagaoka |jj and to Lieb ||. Nagaoka's 
theorem states, that creating a single hole in a half-filled connected lattice for infinitely repulsive interaction renders 
the ground state ferromagnetic. Lieb's theorem is valid for arbitrary finite repulsion. It states that on a bipartite 
lattice at half- filling the ground state has spin S = h\\B\ — |^4||, where \B\ (\A\) is the number of sites in the B 
(A) sublattice. For reviews of rigorous results about the Hubbard model in arbitrary dimensions see Some 
simplifications occur in the limit of infinite lattice dimension [p|-ji"l|. However, most exact results have been obtained 
for the one-dimensional lattice. This is because a complete set of eigenf unctions of the Hubbard Hamiltonian is only 
known for this case. The one-dimensional Hubbard model was solved by Lieb and Wu 0]. They used the nested 
Bethe ansatz discovered in ]L^ , |l3[ ], 

There exists a vast literature on the Hubbard model. Some of the most significant results have been collected in 
two reprint volumes. The volume [|5| gives a general overview for the years 1963-19900 For a review of the history 
of the exact solution in one dimension including a rather exhaustive list of references until about 1992 we refer the 
reader to 

Most of the references listed in [|~7j are based on the seminal 1968 paper (lj] by Lieb and Wu. In this paper the 
problem of diagonalizing the Hamiltonian was reduced to solving a set of coupled nonlinear equations known as the 
Bethe ansatz or Lieb-Wu equations. Lieb and Wu calculated the ground state energy of the system. They showed 
that the model at half-filling is an insulator for arbitrary positive value of the coupling U. In other words, they showed 
that the half-filled model undergoes a Mott transition at critical coupling U = 0. 

In 1972 Takahashi proposed a classification of the solutions to the Lieb-Wu equations fl8| , which is commonly 
referred to as 'Takahashi's string hypothesis'. Analogous classifications are used in all models solvable by the Bethe 
Ansatz method (see e.g. 19 2(J for the case of the Heisenberg model and ||l[] for the case of the Anderson model, 
which bears certain similarities with the Hubbard model). The string hypothesis is the basis of many subsequent 
publications. In the paper |18|, Takahashi used it to obtain a set of nonlinear integral equations that determines the 
thermodynamics of the Hubbard model. By solving these equations in some limiting cases, he was able to calculate 
the low temperature specific heat in |2^] . 

At the beginning of the 80's Woynarovich resumed the study of the excitation spectrum of the Hubbard model 
[2||26| which was started ten years earlier |27],|28| . He gave a detailed analysis of the charge excitations at half-filling 
[23] and was the first to study gapped, spin-singlet charge excitations at half- filling. These involve the first examples 
of excitations which in the sequel will be called fc-A-string excitations p3. In his article (2^] Woynarovich presented 
the explicit form of the Bethe ansatz wave function (see equations (^)-^) below). 

Since the publication of the reprint volume |l^ ] there have been several interesting developments. Two of the 
present authors |^9|,[50| showed that the excitation spectrum at half-filling in the absence of a magnetic field is given 
by the scattering states of only four elementary excitations. Two of them carry charge but no spin, the other two 
carry spin but no charge. These elementary excitations are called holon, antiholon and spinon with spin up or down, 
respectively. They form the fundamental representation of SO(4). In the same articles |2^,Q the S'-matrix of the 
four quasiparticles was obtained. Thus there is now a complete and satisfactory picture on the level of elementary 
excitations of spin-charge separation in the one-dimensional Hubbard model at half-filling. Spin and charge degrees of 
freedom can be excited separately, but the corresponding quasiparticles do interact, albeit weakly. The interaction is 
seen in the non-triviality of the 5-matrix. Spin-charge separation is one of the most interesting properties of the one- 
dimensional Hubbard model. Recently experimental evidence was found for the existence of spin-charge separation 
in quasi one-dimensional materials |]3l| , j32"|| . 

Another interesting recent development is the calculation of the bulk thermodynamic properties of the Hubbard 
model within the quantum transfer matrix approach |53|,[34| . In contrast to the traditional approach [^8 35 3^], the 
quantum transfer matrix approach leads to a finite number of non-linear integral equations that determine the Gibbs 
free energy. This enables high-precision numerical calculation of thermodynamic quantities such as the charge and 



^ee also jTJ 
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spin susceptibilities over the entire range of doping, temperature and magnetic field, ft further opens the interesting 
perspective of calculating the correlation length at arbitrary finite temperature. 

It was shown in the papers |57] [59| how to use a pseudo particle approach in order to obtain transport properties 
(optical conductivity) of the one-dimensional Hubbard model. 

There was also progress in the understanding of the algebraic structure of the Hubbard model. Shiroishi and 
Wadati showed p0| , that the i?-matrix, which was constructed earlier by Shastry p|-p3[, and which underlies the 
integrability of the Hubbard model, satisfies the Yang-Baxter equation. Martins and Ramos Q,[l5| were able to 
construct a variant of the algebraic Bethe ansatz for the Hubbard model. They obtained the eigenvalue of the transfer 
matrix of the two-dimensional statistical covering model (see also jf6|). This result was later used in the quantum 
transfer matrix approach to the thermodynamics J34| . Another interesting algebraic result was the discovery of a 
quantum group symmetry of the Hubbard model on the infinite line. The Hamiltonian is invariant under the direct 
sum of two Y(su(2)) Yangians p7| . The relation of these Yangians to Shastry's i?-matrix was clarified in p8| , p)[ , 
where it was also shown that the eigenstates of the Hubbard Hamiltonian on the infinite interval, at zero density 
transform like irreducible representations of one of the Yangians. 

The purpose of this article is to give a pedagogical introduction to the Bethe ansatz solution of the one-dimensional 
Hubbard model and at the same time to fill some gaps in the previous literature. We present a detailed account 
of the Bethe ansatz solution for periodic boundary conditions and of the thermodynamics of the model. There 
are two approaches to the thermodynamics. The approach of Takahashi |T^,|35||36| relics on a string hypothesis for 
the Hubbard model and is a natural generalization of Yang and Yang's thermodynamic Bethe ansatz for the delta 
interacting Bose gas [po| . The second approach J34| is built on a lattice path integral formulation of the partition 
function. We compare both approaches and discuss their specific advantages. Special attention is given to an aspect 
which, although fundamental, was largely ignored in the previous literature, namely fc-A strings. fc-A strings are 
spin-singlet bound states of electrons. 

The Bethe ansatz for the one-dimensional Hubbard model |l4| gives the eigenf unctions and eigenvalues of the 
Hubbard Hamiltonian parametrized by two sets of quantum numbers {kj} and {A/}, which are solutions of the Lieb- 
Wu equations (see formulae (0), (@) below). The kj and A/ are called charge momenta and spin rapidities, respectively. 
The Lieb-Wu equations have finite and infinite solutions kj, A;. They should be considered separately, because they 
have different occupation numbers. Every finite solution kj (or A;) can be occupied only once. If two finite kj (or 
A;) coincide, the wave function vanishes (see formulae (pi]), ( p"2| ) and below). By contrast, infinite kj (or A;) can be 
generally occupied more than once. Later we shall explain that the multiplicities of occupation of the infinite fc's and 
A's are given by the dimensions of the representations of corresponding su(2) symmetry algebras. In order to study 
this carefully the 'regular' Bethe ansatz was defined in |5lfl . 'Regular' means that all fc's and A's are finite. They may 
be real or complex. 

Takahashi's string hypothesis jl8| is a statement about the structure of the regular solutions of the Lieb-Wu 
equations in the thermodynamic limit. Except for the real solutions (all kj and all A; real) there are solutions 
involving complex kj and A;. The complex momenta and rapidities occur in two kinds of configurations, which are 
symmetric with respect to the real axis. These configurations are called strings. There are A strings involving only 
spin rapidities and k-A strings, which involve charge momenta as well. The A strings can be interpreted as bound 
states of magnons, whereas the fc-A strings describe spin singlet bound states of electrons. 

The A strings in the Hubbard model are similar to the A strings in the isotropic Heisenberg spin chain, which have 
been extensively studied in the literature [^2| , ^9|j5^j5^] . Much less attention has been given to the fc-A strings, which 
are peculiar to the Hubbard model. They play an important role at half-filling, where the fc-A two string enters the 
calculation of the phase shift in the holon-holon scattering |29|]3C| l. Holons are the lowest lying charge excitations of 
the Hubbard model. At half-filling they have a gap. Below half-filling they are gapless, however, whereas all fc-A 
strings lead to gapped excitations in the thermodynamic limit (see section VI). Hence, fc-A strings below half- filling 
do not contribute to the low energy properties of the model. What is their physical significance then? They do 
contribute to the high temperature thermodynamic properties of the Hubbard model. This is, in fact, the context in 



which they were first introduced |18|. On the other hand, fc-A strings are interesting as a curious kind of excitations, 
which does not exist in more simple integrable models. 

According to the string hypothesis the fc-A strings approach certain ideal configurations as the number of lattice 
sites becomes large. We call these configurations 'ideal fc-A strings'. They are characterized by to complex A's and 
2m complex fc's. The A's involved in an ideal fc-A string have common real part A', and their imaginary parts are 
(to — f)^, (to — 3)^, . . . , —(to — 1)^. In the repulsive case {U > 0) the fc's are given as 

fc 1 = 7r - arcsin(A' + to^) , 



■5 



k 2 = arcsin(A' + (m-2)i^) , 

k 3 = 7T - k 2 , 



fc 2m - 2 = arcsin(A'-(m-2)f ) , 
k 2m - 1 = it - k 2m ~ 2 , 

k 2m = 7T - arcsin(A' - m^) . 

This means that Re sin(fc a ) = A'. Hence, the sin(fc a )'s and the A's form a string in the complex plane. 

In the thermodynamic limit, when the strings become ideal, the variables describing their width can be eliminated 
from the Lieb-Wu equations. We call the resulting set of equations discrete Takahashi equations. 

Some reformulation of the string hypothesis may be necessary before it will be possible to achieve a rigorous 
mathematical proof. Yet, the string hypothesis has passed many tests, and there is no doubt by now that it describes 
the physics of the Hubbard model correctly. We summarize our understanding of the issue and present several 
important tests and consequences of the string hypothesis. We shall mostly concentrate on the Hubbard model below 
half-filling, since the case of half- filling was treated elsewhere ^|3(J . 

Let us outline the plan of this article. 

In section II we summarize known basic results about the Hubbard model. This section contains a review of the 
Bethe ansatz solution of Lieb and Wu [Q . We present the wave function in the form given by Woynarovich |2^] and 
discuss its discrete symmetries, namely the symmetries under permutations of electrons and quantum numbers and 
the particle-hole and spin-reversal symmetries. This leads to the notion of regular Bethe ansatz states. We proceed 
with explaining the SO(4) symmetry |55|-|59[|, which is characteristic of the Hubbard model. Then we introduce the 
discrete Takahashi equations. SO (4) symmetry and discrete Takahashi equations are the prerequisites for the proof 



of completeness of the Bethe ansatz 1 60 , which is reviewed at the end of the section. 

Section III comprises a rigorous analytical study of the Lieb-Wu equations for one down spin. For the sake of 
pedagogical clarity we mostly focus on the cases of two and three electrons. This is our first test of the string 
hypothesis. The result is positive: k-A strings do exist as solutions of the Lieb-Wu equations. They become ideal in 
the thermodynamic limit. Furthermore, the counting of solutions implied by the discrete Takahashi equations agrees 
in this case with the counting obtained directly from the Lieb-Wu equations. 

Section IV is devoted to the self-consistent solution of the Lieb-Wu equations for three electrons and one down spin. 
Self consistency arguments underly the derivation of the discrete Takahashi equations. In the simple case considered 
in this section we can be more explicit. We calculate the deviations of the rapidities and momenta (solutions of the 
Lieb-Wu equations) from their ideal positions (solutions of the discrete Takahashi equations). It turns out that these 
deviations vanish exponentially in the thermodynamic limit. This is our second positive test of the string hypothesis. 

In section V we complement the analytical considerations of the previous sections with numerical results. Numerical 
data based on the Lieb-Wu equations are compared with data, which were obtained independently of the Bethe 
ansatz by direct numerical diagonalization of the Hamiltonian. We find perfect agreement between the two numerical 
methods. The energy levels obtained by the two methods agree within a numerical error of 0(1CP 15 ). Our numerical 
study confirms the completeness of the Bethe ansatz and the correctness of the counting of the solutions implied by 
the string hypothesis. This is our third successful test of the string hypothesis. 

In section VI we review Takahashi's approach to the thermodynamics of the Hubbard model fll8|j . We show 
that the dressed energies of all fc-A strings (bound states of electrons) follow from Takahashi's integral equations 
(thermodynamic Bethe ansatz equations) in the zero temperature limit. We can actually do better. Starting from the 
thermodynamic Bethe ansatz equations and passing to the zero temperature limit we obtain a complete classification 
of all elementary excitations at zero magnetic field, below half-filling. 

Takahashi derived his equations in order to calculate thermodynamic quantities such as the specific heat or charge 
and spin susceptibilities for the Hubbard model |lj,|2^,^l|,^2| . Nowadays there is an independent method to calculate 
these quantities, which does not rely on strings. It is called the quantum transfer matrix method (6^Q,Q. In section 
VII we compare the results of both methods and find that they agree well. The string hypothesis also passes this 
significant test. If the string hypothesis would miss one of the elementary excitations, it should be visible in the 
thermodynamics. 

Section VIII contains a brief conclusion and a list of interesting open problems. 
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In appendix A we present a derivation of the Bethe ansatz wave function for the Hubbard model. It can be seen 
from the derivation that charge momenta and spin rapidities do not have to be real. To every solution of the Lieb-Wu 
equations there corresponds a well defined periodic wave function. This is particularly the case for the fc-A string 
solutions. 

In appendix B we derive the algebraic Bethe ansatz solution of the inhomogeneous isotropic Heisenberg model, 
which is needed to construct the Bethe ansatz wave function in appendix A. 

Appendix C contains the tables of our numerical data for three electrons and one down spin. 



II. BETHE ANSATZ FOR THE HUBBARD MODEL 



A. Eigenfunctions and eigenvalues 

The Hamiltonian of the one-dimensional Hubbard model on a periodic L-site chain may be written as 

L L 

^ = -EE ktvCj+i,* + 4+1,**,*) + u E( n JT - h)( n oi - 1) ■ (!) 

J=l<r=T,J. J'=l 

c~j~ a and Cj. CT are creation and annihilation operators of electrons in Wannier states, and periodicity is guaranteed by 

setting CL+i,a = ci,a- Tij tl7 — Cj a Cj.a is the particle number operator for electrons of spin a at site j , U is the coupling 
constant. The eigenvalue problem for the Hubbard Hamiltonian (Q) was solved by Lieb and Wu H| using the nested 
Bethe ansatz . The Hubbard Hamiltonian conserves the number of electrons N and the number of down spins 
M. The corresponding Schrodinger equation can therefore be solved for fixed N and M. Since the Hamiltonian is 
invariant under particle-hole transformations and under reversal of spins fy}| , we may set 2M < N < L. We shall 
denote the positions and spins of the electrons by Xj and <jj, respectively. The Bethe ansatz eigenfunctions of the 
Hubbard Hamiltonian ([j]) depend on the relative ordering of the Xj. There are N\ possible orderings of the coordinates 
of N electrons. Any ordering may be related to a permutation Q of the numbers 1, . . . , N through the inequality 

1 < Xqx < xq 2 < • • • < xq N < L . (2) 

This inequality divides the configuration space of N electrons into iV! sectors, which can be labeled by the permuta- 
tions Q. The Bethe ansatz eigenfunctions of the Hubbard Hamiltonian (|l|) in the sector Q are given as 

( N \ 

ip(xi,...,x N ;a 1 ,...,a N )= 2J sign(PQ) <Pp(<tqi, a QN ) exp I k Pj x Q A . (3) 

Pes N \ j=i ) 

Here the P-summation extends over all permutations of the numbers 1, . . . , N . These permutations form the symmetric 
group Sn- The function sign(<5) is the sign function on the symmetric group, which is —1 for odd permutations and 
+1 for even permutations. The spin dependent amplitudes (pp((7Qi, . . . , <jqn) can be found in Woynarovich's paper 
[ p3| . They are of the form of the Bethe ansatz wave functions of an inhomogeneous XXX spin chain, 

M 

<Pp(<tqi, ■ ■ ■ ,oqn) = ^2 A(\ nl ,...,\ irM )Y[F P (\ nl ;y l ) . (4) 

ttESm 1=1 



Here Fp(X; y) is defined as 



F q s 1 tt A-sinfcpj -iU/A 

P[ ,V > A-sinfcp.„ +i[//4 J-J- A-sinfcp, +1C//4 ' { > 



sin kp y + iU/A r~ A — sin kpj + iC//4 
and the amplitudes A(Xi, . . . , Am) are given by 



A(M,...,X M )= Xm - X "_\ iU/2 . (6) 

l<m<n<M 



y,j in the above equations denotes the position of the jth down spin in the sequence <tqi, . . . , <jqn- The y's are thus 
'coordinates of down spins on electrons'. Below we shall illustrate the notation through an explicit example. 
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The wave functions (|J) are characterized by two sets of quantum numbers {kj} and {A;}. These quantum numbers 
may be generally complex. The kj and Aj are called charge momenta and spin rapidities, respectively. The charge 
momenta and spin rapidities satisfy the Lieb-Wu equations 



>■ >■ If A ' W '' U ' ' j = 1 N (7) 



^ A; - sin kj - ill/ 4 yr A; - A m - iU/2 

H A,-sin/L+k7/4 H Xi - X m + iU/2 ••>-»* • W 

A derivation of the wave function (||) and the Lieb-Wu equations (Q), (|8|) is presented in appendices A and B. 

The wave functions (^) are joint eigenfunctions of the Hubbard Hamiltonian (Q) and the momentum operator^ with 
eigenvalues 

N U ( N \ 

E = -2^cos kj + —{L - 2N) , P = ^ kj mod 2ir . (9) 

j=i \i=i / 

The 'coordinates of down spins' yj which enter (Q) depend on (<7i, . . . , o\zv) cm<i on (xi, . . . , xn). The follow- 
ing example should help to understand the notation. Let L = 12, N = 5, M = 2, and let, for example, 
(a*,..., a*) = (7,3,5,1,8), (ai,...,a 5 ) = (TTTU)- Then x 4 < x 2 < x 3 < x x < x 5 , i.e. Q = (4,2,3,1,5). It 
follows that (xqi,. . .,x Q5 ) = (1,3,5, 7,8) and (a Q i, . . -,a Q5 ) = (ITTTD- Thus y x = 1, y 2 = 5. 

Whenever it will be necessary, we shall indicate the dependence of the wave functions (^) on the charge momenta 
and spin rapidities by subscripts, iji = V , fei,...,fc JV ;Ai,...,A M ■ Let us consider the symmetries of the eigenfunctions under 
permutations, 

i>(x P1 ,...,xp N ;(Tp 1 ,...,ap N ) =siga(P)tp(x 1 ,...,x N ;a- 1 ,...,a N ) , P e S N , (10) 
i>kp 1 ,...,k PN ;X 1 ,...,x M = sign(P)V'fe 1 ,...,fc J v;A 1 ,...,A M ! p e S N , (11) 

^k u ...,k N ;Xpx,...,X PM = 1pk 1 ,-,k N ;X 1 ,..,,X M i P ^ S M ■ (12) 



Equation ( |10| ) means that the eigenfunctions respect the Pauli principle. ( |1 1| ) and (|TJ) describe their properties with 
respect to permutations of the quantum numbers. They are totally antisymmetric with respect to interchange of the 
charge momenta kj, and they are totally symmetric with respect to interchange of the spin rapidities A/. Hence, in 
order to find all Bethe ansatz wave functions we have to solve the Lieb-Wu equations (0), (||) modulo permutations 
of the sets {kj} and {A;}. The kj's have to be mutually distinct, since otherwise the wave function vanishes due to 

§. In fact, the Az's have to be mutually distinct, too. This is called the 'Pauli principle for interacting Bosons' (see 
). We would like to emphasize that there are no further restrictions on the solutions of ([?]), (||). In particular, the 
spin and charge rapidities do not have to be real. 

Bcthe ansatz states on a finite lattice of length L that have finite momenta kj and rapidities A;, a non- negative 
value of the total spin (N — 2M > 0), and a total number of electrons not larger than the length of the lattice (N < L) 
are called regular (cf. pi] , page 562). 

There exist two discrete symmetries of the model which can be used to obtain additional eigenstates from the regular 
ones [ pT| . The Hamiltonian is invariant under exchange of up and down spins. This symmetry allows for obtaining 
eigenstates with negative value N — 2M of the total spin from eigenstates with positive value of the total spin. This 
symmetry docs not affect the number of electrons. Thus, its action on regular states does not lead above half filling. 
States above half filling (N > L) can be obtained by employing the transformation Cj a — > ( — l) J 'c^, Cj a — > (— iycj a , 
a =1, 1, which leaves the Hamiltonian ([[]) invariant, but maps the empty Fock state |0) to the completely filled Fock 
state | Tl)- 



2 For a proper definition of the momentum operator see appendix B of IK 
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B. SO(4) symmetry 

The Hubbard Hamiltonian (|l]) is invariant under rotations in spin space. The corresponding su(2) Lie algebra is 
generated by the operators 

L L L 

j=l 3=1 3=1 

[C,C + ] = -2C 2 , [C,C Z ] = C , [C + ,C] = ~C f • 

For lattices of even length L there is another representation of su(2), which commutes with the Hubbard Hamiltonian 
[HI 0' This representation generates the so-called 77-pairing symmetry, 

L L L 

j'=l j'=l j'=l (- 14 ) 

[jj, ryt] = -2ry 2 , [77,77*] =77 , M, »7 Z ] = -?7 t ■ 

The generators of both algebras commute with one-another. They combine into a representation of su(2)0su(2). 

The 77-pairing symmetry connects sectors of the Hilbert space with different numbers of electrons. The operator 
rj\ for instance, creates a local pair of electrons of opposite spin and momentum it. Hence, in order to consider the 
action of the ^-symmetry on eigenstates we write them in second quantized form. 

L 

\k 1 ,...,k N ;\ 1 ,...,\ M ) = 2J i 3 ki,...,k N ;X 1 ,...,x M (xi,...,XN;ai,...,a N )c^ uai ...c^ N>aN \0) , (15) 

a:i,... ) xjv=l 

where a\ = . . . = um =1 and <tm+i = • • ■ = ojv =T- It is easily seen that 

(C +rf)\ki,..., k N - Ai, . . . , A M ) = (M - f . . . , k N ; Ai, . . . , A M ) • (16) 

Here M — ^ is integer, since L is even. Therefore the symmetry group generated by the representations (|ll|), ( fl4| ) is 
SO(4) rather than SU(2)xSU(2) [g. 

It was shown in |51| that the regular Bethe ansatz states are lowest weight vectors of both su(2) symmetries ( |l3| ) 
and®, 

C\kx, . . . , k N ; Ai, . . . , Am) = , 77^1, . . . , k N ; Ai, . . . , Am) = . (17) 

This is an important theorem. It was the prerequsite for the proof of completeness (see section D) of the Bethe ansatz 
for the Hubbard model in [pC|. The proof of ( |l7j) is direct but lengthy jplj . ( and 77 are applied to the states (|l5|), 
and the Lieb-Wu equations (j? ) , (§) are used to reduce the resulting expressions to zero. We would like to emphasize 
that the proof of is not restricted to real solutions of the Lieb-Wu equations. It goes through for all solutions 
corresponding to regular Bethe ansatz states including the strings. 

Since the two su(2) symmetries (|l3|), (|lj) leave the Hubbard Hamiltonian (Q) invariant, additional eigenstates which 
do not belong to the regular Bethe ansatz can be obtained by applying £T and to regular Bethe ansatz eigenstates. 
Since 

C z |fci,...,/cAr;Ai,...,A M ) = (Af- f )|fc 1 ,...,fc A r;Ai,...,A M ) , (18) 
77*|fci,...,A;jv;Ai,...,AM) = \{N - L)\kx, . . . ,k N \\x, . . . ,\ M ) , (19) 

a state |fci, . . . , fczv; Ai, . . . , Am) has spin |(JV — 2M) and 77-spin i(L — N). The dimension of the corresponding 
multiplet is thus given by 

diniM.JV = (N-2M+ 1)(L - N + l) . (20) 

The states in this multiplet are of the form 

\k 1 ,...,k N ;X 1 ,...,\ M ;a;(3) = (tf) a (rff\h, . . . , k N ; A x , . . . , A M ) , (21) 

where a = 0, . . . , N — 2M and (3 = 0, . . . , L — N. Note that states of the form (^lj) can be obtained from regular 
Bethe ansatz states with N > N, M > M by formally setting some of the charge momenta and spin rapidities equal 
to infinity p3Mp0|. 
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C. Discrete Takahashi equations 



Let us now formulate Takahashi's string hypothesis Jul more precisely: All regular solutions {kj}, {A/} of the 
Lieb-Wu equations (Q), (§h consist of three different kinds of configurations. 

(i) A single real momentum kj. 

(ii) m A's combining into a A string. This includes the case m — 1, which is just a single A Q . 

(iii) 2m k's and m A's combining into a fc-A string. 

For large lattices (L ^> 1) and a large number of electrons (N ;§> 1), almost all strings are close to ideal, i.e. the 
imaginary parts of the fc's and A's are almost equally spaced. 
For ideal A strings of length m the rapidities involved are 

A^ = A™ + (m-2j+l)f . (22) 

Here a enumerates the strings of the same length m, and j = 1, . . . , m counts the A's involved in the ath A string of 
length m. A™ is the real center of the string. 

The fc's and the A's involved in an ideal fc-A string are (for U > 0) 

k x a — 7r - arcsin(A'™ + m^-) , 

k 2 a = arcsin(A'™ + (m-2)f ) , 

fc^ = vr — k 2 a , 



: (23) 
fc2™- 2 =arcsin(A™-(TO-2)f) , 

h 2m-l _ _ iSLm-2 
ot ot > 

k 2 a m = it - arcsin(A'™ - mf ) , 

and 

A™ =A'r + (m-2j + l)f . (24) 
Again m denotes the 'length' of the string, a enumerates strings of length m, and j counts the A's involved in a given 



string. A'™ is the real center of the fc-A string. The branch of arcsin(x) in (23) is fixed as — n/2 < Re(arcsin(a;)) < n/2. 

The string hypothesis assumes that almost all solutions of the Lieb-Wu equations (^) , (^) are approximately given 
by (E3)-(p^) with exponentially small corrections of order 0(exp(— SL j), where 5 is real and positive and depends on 
the specific string under consideration. 

Using the string hypothesis inside the Lieb-Wu equations (Q), (||) and taking logarithms afterwards, we arrive at 
the following form of Bethe ansatz equations for strings, which we call discrete Takahashi equations 



n— 1 a—1 



N-2M' /Kin ■ i \ oo M' m /A'" _ A' m \ 

L[arcsi n (A': + nf)+arcsin«-nf)] = 27 rJ-+ £ f ° ! ' +EE 9 - • ( 27 ) 

j=l ^ I ' m=l /3=1 \ / / 

Here we assumed L to be even. Zj, J™, and J'™ are integer or half-odd integer numbers, according to the following 
prescriptions: ij is integer (half odd integer), if ^ m (M m + M' m ) is even (odd); the J™ are integer (half odd integer), 
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if N — M„ is odd (even); the J'™ are integer (half odd integer), if L — (N — M' n ) is odd (even). M n and M' m are 
the numbers of A strings of length n, and fc-A strings of length m in a specific solution of the system (|25|)-(|2^). 
M' = y^ ^L] nM' n , is the total number of A's involved in fc-A strings. The integer (half-odd integer) numbers in 
2a)- (E7p have ranges 



L T L 



if °° 

\ m— 1 



M m - 1 



nm lvJ -m 




N - 



2M' - t nm M' m 



(28) 
(29) 

(30) 



where t r 



2 min(m, n) 



The functions 9 and m „ in ( 25 )- (|27|) are defined as 6{x) = 2arctan(a;), and 



@ram («^) 



26> 



20 (? 



20 



|ra - m| + 2 
■■■ + 29 



20 



x 



(-) 

\2nJ 



m-2 
) , if n = 



n + m 



, if n =/= m, 



(31) 



m. 



In terms of the parameters of the ideal strings total energy and momentum (|9j) are expressed as 

oo Mi 



P = 



N—2M' 



^ ^ (2 Re arcsin (A^ + n^) - (n + 1)tt) 



n— 1 a—1 



mod 2ir 



N — 2M' 



ML 



E=-2 ]T cos(fc,)+4^^Re V /l-(A'S+nf) 2 + ^(i-2iV) 



(32) 



(33) 



Equations (|25|)-(|30D can be used to study all excitations of the Hubbard model in the thermodynamic limit. They 
are the basis for the derivation of Takahashi's integral equatio ns pjj ], which determine the thermodynamics of the 
Hubbard model (see sections VI and VII). Applications of (|25|)-(|30|)~are usually based on the following assumptions. 



(i) Any set of non-repeating (half odd) integers Ij, J", J' subject to the constraints ([28|)-(|30|) specifies one and 



only one solution {kj}, {A™}, {A'™} of equations (|25|)-(|27|). 

(ii) The solutions {fc 7 j, {A^}, {A'^} of (p5|)-(|27|) specified by a set of non-repeating (half odd) integers Ij, J™, J'™ 
subject to (|2^)-(po|) are in one-to-one correspondence to solutions of the Lieb-Wu equations (0), (||). 

(iii) For large L and V almost every solution {kj}, {A;} of the Lieb-Wu equations (Q), (|8|) is exponentially close to 
the corresponding solution {kj}, {A™}, {A'™} of the discrete Takahashi equations, which means that the strings 
contained in {kj}, {A;} are well approximated by the ideal strings determined by {kj}, {A™}, {A'™}. 



D. Completeness of the Bethe ansatz 

The proof of completeness of the Bethe ansatz given in |3(| is based on assumptions (i) and (ii) above. Similar 
assumptions were proved for other Bethe ansatz solvable models ]65| . Note that assumption (ii) does not mean that 
the classification of the solutions of the Lieb-Wu equations ([?]), (|8|) into strings is actually given by (|25|)-(|3"o|). There 
may be a redistribution between different kinds of strings. This phenomenon was observed in a number of Bethe 
ansatz solvable models and was carefully studied by examples [ ^2p^j60| ] (see also section V.B). It turned out that the 
redistribution did in no case affect the total number of solutions of the Bethe ansatz equations. 

Using (i) and (ii) above, the proof of completeness reduces to a combinatorial problem based on (^8|)-(^o|) 



From (|2§|)-(30) we read off the numbers of allowed values of the (half odd) integers Ij, J™, J' a in a given configuration 



{M n }, {M' n } of strings. These numbers are 
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(i) L for a free kj (not involved in a fc-A string), 

(ii) N — 2M' — Y2m=i tnm,M m for a A string of length n, 

(iii) L-N+ 2M' - Em=i tnmM' m for a fc-A string of length n. 

The total number of ways to select the Ij, J™, J'™ (recall that they are assumed to be non-repeating) for a given 
configuration {M n }, {M' n } is thus 

,< , r, / ,s / ^ \ tt {N-2M'-Y°° , t nm M m \ ^ (L- N + 2M'-Y°° , i nm M' \ , , 

n({M„},{M' n })=( iV _ 2M ,)n( M* ) ■ (34) 

Hence, the number of regular Bethe ansatz states for given numbers ./V of electrons and M of down spins is 

n ieg (M,N)= Yl n({Mn},{K)) > (35) 

{m„},{m;} 

where the summation is over all configurations of strings which satisfy the constraints N — 2M' > and M — 
Ym=i m {M m + M^). Finally, the total number of states ( pl| ) in the SO (4) extended Bethe ansatz is 

rhat(L) = «rc g (Af, N) dim M , N = ^ cg (M, N)(N - 2M + l)(L - N + 1) , (36) 

M,N M,N 



where the sum is over all M, N with < 2M < N < L. The sums ( |35| ) and (|36| ) were calculated in |60[|. It turns out 
that 

n tot (L)=4 L , (37) 

which is the dimension of the Hilbert space of the Hubbard model on an L-site chain. 
Let us list again the essential steps that led to the above proof of completeness: 

(i) Impose periodic boundary conditions. 

(ii) Take Woynarovich's wave function (§-(§). 

(iii) Define the Bethe ansatz in the narrow sense of regular Bethe ansatz (see below ([l2|)). This eliminates infinite 
fc's and A's whose multiplicities are not under control. 

(iv) Prove the lowest weight theorem ([l7]). Then gluing back solutions with infinite fc's and A's is equivalent to 
considering the multiplets (J2l|). 



(v) The multiplicities of occupation of infinite fc's and A's are given by the dimensions ( p0[ ) of the multiplets (|21|). 

(vi) Use Takahashi's integers (pq)-(pO|) for counting of the regular Bethe ansatz states. 

III. LIEB-WU EQUATIONS FOR A SINGLE DOWN SPIN (I) - GRAPHICAL SOLUTION 

In this section we study the Lieb-Wu equations (^), (||) in the most simple non-trivial case, when there is only one 
down spin, M = 1. For pedagogical clarity some emphasis will be on the most instructive cases N = 2 and N = 3. 
These are the cases which we also studied numerically (cf. section V). Some of the analytical calculations, however, 
are presented for general N, simply because the general arguments are simple enough and enable treating the cases 
N = 2 and N = 3 to some extent simultaneously. 

The Lieb-Wu equations for N = 2 and M = 1 were studied before in appendix B of [^0|. There the emphasis 
was on the redistribution phenomenon mentioned in section II. D. For U = the Hubbard Hamiltonian turns into a 
free tight binding Hamiltonian, and there is no bound state of electrons (fc-A string) left. It is therefore clear that 
bound states decay as the coupling becomes weaker. In appendix B of |6C| ] it was shown that each time a fc-A string 
disappears from the spectrum at a certain critical value of the coupling U > 0, a new real solution emerges. Here we 
take a slightly different point of view. We fix U and study the solutions for large finite L. It turns out that there 
is no redistribution phenomenon for the most simple fc-A strings consisting of two complex conjugated fc's and one 
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real A as L — > oo. These strings always exist for large enough finite L, and their number is in accordance with the 
counting implied by Takahashi's discrete equations ( p5|)-(|30| ). In this respect the fc-A strings of the Hubbard model 
are different from the A strings in the XXX spin chain |j4|. 
For M = 1 the Lieb-Wu equations (g), @ read 

ji, t A - sin kj — if/4 . , T ,„_. 

e 3 =1 • , , , . J = 1,...,JV , (38) 



A - sin kj + iU/A 

sin kj - ill /A 
A - sin kj + ill/ A 



Yj A — sin fcj - if/4 



Equation ( |39"| ) can be replaced by the equation for the conservation of momentum, 

e i(fei+...+fejv)L _ 2 _ r^Q\ 



38J) and (|40|) follow from ( |38[ ) and (39) and vice versa. Let us take the logarithm of (|40|) and solve ( J3q ) for sinfcj — A. 
We obtain the equations 

(fci + . . . + fc^) mod 2-7T = — — , ra = 0, . ..,L-1 , (41) 

L 

sinfc, -A= ^ctgf^ ) , j = l,...,JV , (42) 



which are equivalent to fl38j) and (p9j) but more convenient for the further discussion. 



A. All charge momenta kj real 



The equation 



a U fqL 
sin 9- A = — ctg I — 



(43) 



is easily solved graphically for g as a function of A (see figure 1). It has at least L branches of solutions qe(A) belonging 
to the interval < q < 2ir. There is at least one branch with (£ — 1)^ < q < t^-. Yet, for ^ < q < ^ there may 
be more than one solution in the interval [(£ — 1)^,£^], if U or L is too small. We have the following uniqueness 
condition, 



— 1 = mm oUsmo > max o„ — ctg — = max 

Q<q<2n H 0<q<2ir H A \ 2 J 0<q<2% 



UL 



qL 
2 



UL 
~8~ 



(44) 



which is equivalent to 



L > 



U 



(45) 



We call this condition Takahashi condition. In the following we will assume the Takahashi condition to hold. 
Equation ( ^3| ) defines A as a function of q. We have 



dA 

— = cos q 
dq 



UL 



■ sin 



qL 



Then, using (|4q), ^ > 0. Hence, all branches qe(A) of solutions of ( |43] ) are monotonically increasing, 
can summarize the properties of the solutions of equation (f43|) in the following lemma. 



(46) 

§M<M > o. We 



Lemma III.l If L > fj, equation (U5l) has exactly L branches of solutions qe(A), which have the properties 
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(<) 
(it) 
(Hi) 



(^-1)^<^(A)<^ 
dq e (A) 



dA 



> 



(47) 



lim qe (A) = (i-l)^ 

A— >— oo _L 



lim g^(A) = 

A^oo 



' L 



(iii) can be seen from figure 1. 











sin(q) - A 


-A 




2it \ 
L \ 


4it \ 
L \ 


L 









FIG. 1. Sketch of equation (^). 



Lemma III. 1 is sufficient to characterize and count the real solutions of ( pSf ) and (|3^) (or equivalently (^l|) and 
We are seeking for solutions of ( |38| ) and (|39|) modulo permutations, where all kj are mutually distinct. Choose 
N branches qe 1 (A) < . . . < q^ N (A) of solutions of ( (43|) , and define 

Q(A) = % (A) + ... + g^(A) . (48) 
Then kj = q^ (A), j = 1, . . . , TV, and A solve (||) and (§|), if and only if 



2tt 

Q(A) mod 2vr = m — . (49) 
L 



Now, using Lemma III. 1 , 



2-7T 2ir 
lim Q(A) = (£ 1 + ...+£ N -N)— , lim Q(A) = (h + . . . + t N )— . (50) 

A — * — OO Li A — >oo _L 

Furthermore, d< ^j^ > 0. Thus there are precisely N — 1 values A Q , a = 1, . . . , N — 1, which satisfy (|49| ) for a given 
choice l\ < . . . < In of branches of equation ( |43| ) (recall that we exclude A = ±oo). We summarize our result in the 
following 

Lemma III. 2 Let L > |j. TTien i/iere are precisely (^)(N — 1) solutions with all kj real to equations $3f\), ^3^)- To 
every choice of N mutually distinct intervals Tj = [(£j — l)^-,£j^-], j = 1, . . . , N, £j ^ £k, there correspond N — 1 
solutions with kj £ Tj. These solutions are characterized by N — 1 different values of A. 

Let us show that our result coincides with Takahashi's counting (|2^)-(|3(]). We have N real fc's and one A string 
of length 1. Since there is no fc-A string, there is no J' n a to specify. M' — 0, Mi = 1, Mj = Q for j > 1. Thus 
^ 1 < 5 (N — tu — 1) = \ (N — 2) , and there are N — 1 possible values of J 1 . The number of different sets {7i , . . . , 7jv} 
follows from @ as (£). This means that Takahashi's counting predicts a total number of (AT — 1) real solutions, 
which is in accordance with Lemma III. 2 as long as the Takahashi condition (H^) is satisfied. In the special cases 
N = 2 and N = 3 we find (^) and 2(g) solutions, respectively. 
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B. fe-A two string 



Let us consider equation fi2) in the case that two of the kj's are complex conjugated and the others are real. We 
may set fc_ = k\ = q — i£, k + = k^ = q + i£ with real q and real, positive £. It follows from that 



sin(q + i£) = sin(g)cosh(£) + icos(g)sinh(£) = A 

or, if we separate real and imaginary part of this equation, 

U sin(qi) 



U sin(gi) — isinh(£L) 
4 cosh(£L) — cos(qL) 



sin(<7)cosh(£) = A + 
cos(<7)sinh(£) : 



4 cosh(£L) — cos(qL) 

U sinh(£L) 
4 cosh(£L) — cos(qL) 



(51) 

(52) 
(53) 



Note that for £ = equation fl53|) is satisfied identically and ( |52| ) turns into ((43|). 
Let us consider the two-electron case N — 2 first. Then, by equation j4l|), 

q = m- , m = 0,...,2L-l . (54) 
L 

Equation (|53| ) determines £ as a function of g = JTifj, and equation (|5^) determines A. We have qL = mir. Hence, 
sin(qL) = 0, cos(qL) = (— l) m , and (J52|) and (53) decouple into 



sinh(£) 



A = sin(g)cosh(£) 

U sinh(££) 



4cos(g) cosh(£L) - (-1)' 



(55) 
(56) 



This decoupling is a peculiarity of the two particle case, which makes it more simple than the general case. We can 
discuss equation (p6[) graphically (see figure 2). 




FIG. 2. Sketch of equation <M 



Let us define 



f(0 



sinh(£L) 



th 



cth(^ 



cosh(£L) - (-l) m 

/(£) > for all £ > 0. Hence, ( [So] ) can have solutions for positive £ only if 



for m odd 
for m even 



(57) 
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<jr 37T 7T 

- < q < — , for U > , - < |g - tt| < tt , for [/ < . (58) 



For the sake of simplicity let us concentrate on the repulsive case U > 0. It follows from (p7|) that (|56| ) always has 
exactly one solution £ m for every even m which satisfies ((5^). The condition for a solution with odd m to exist is that 
the derivative of the right hand side of equation (|56|) is larger than the derivative of the left hand side of equation 
as £ approaches zero from the right, i.e. 



UL 



8cos(g) 



> 1 . (59) 



This is satisfied for all q with 5 < q < tF, if and only if 



L>- . (60) 

Again we have found the Takahashi condition fl45|). If the Takahashi condition is satisfied, then there is one and 
only one fc-A two string solution for every m satisfying ( |54| ) and (|5^), and we can easily count these solutions. Their 
number as a function of L is different for odd and even L, respectively. 

Lemma III. 3 Let L > y , U > 0, N = 2. Equations ( pte[ ) and ( p?^ ) have k-A two-string solutions only if the 
momentum q of the center of the string is in the range ^ < q < 4f . The allowed values of q in that range are 
quantized as q = m^. For every q — m,j- there is one and only one k-A two string. The total number of k-A two 
strings is L for L odd and L — 1 for L even. 

Comparing our result with the prediction of Takahashi's counting (p^)-(|30|) we find again agreement. Now there is 
no free kj and no A string, thus no Lj and no J™. Furthermore, M[ = M' — 1, and Mj = for j > 1. Thus 
\J rl \ < h(L — *n — 1) = h{L — 2), which means that there are L — 1 possible values of J' 1 . This agrees with our 
Lemma since L was assumed to be even in (|2^)-(|3C|). 

Note that limL^oo /(£) = 1, pointwise for all £ > 0. This suggests the notion of an ideal string determined by the 
equations 

U 7T 

A = sin(g)cosh(f) , sinh(£) = - 4cQg ^ > 1 = m j ■ (61) 



Replacing ( |55| ) and ( pq ) by (61) means to replace the curves denoted by 'm odd' and 'm even' in figure 2 by the 
dashed line. Clearly, solutions of (|55|), ( |56"| ) are in one-to-one correspondence to solutions of (|6l]), if the Takahashi 
condition is satisfied. We can formulate the following corrolary of Lemma III. 3. 

Lemma III. 4 Let L > y, N = 2. Let ni^ be a sequence of integers < < such that hin^^oo thl J = q, 
■§ ^ 1 ^ According to Lemma III. 3 this defines a sequence £ mL of solutions of (pflj. This sequence has the limit 

lim U = ~ arsinh ( ) , (62) 

i.e. in the thermodynamic limit all k-A two strings are driven to the ideal string positions. 

Proof: (62) follows from (pa), since the sequence £ mi/ is bounded from below. Let us prove the latter statement. 



Assume the contrary. Then there is a subsequence £ mL . of ^ mL which goes to zero, and it follows from (56) 
that lmij^oo /(£ Ki ) = 0. This means (i) is odd for all sufficiently large j, and (ii) lim^oo sinh(£ mi Lj) = 
linij^oo £„ lL . Lj = 0. We conclude that 



lim — = um ~T~ = = um 7~~ — = — o — "TT 7 (63) 

j-> oo smh(£ TOi . L 3 ) j—oo Lj j->oo 8 C0S ( mL ^ £2.) 8cos(g) 



which is a contradiction. Thus, the lemma is proved. 

Let us now proceed with the case N — 3. We have to solve the following system of equations (cf. (|l]), (^2|), (|52"|), 
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2tt . x 

2q + k 3 = m— , m = 0, ...,3L-1 , (64) 

Lj 

Mh)-A=jctg(^\ , (65) 

A = 8in(g)cosh(0 - t UTT^ , T\ . ( 66 ) 

4 cosh(£L) - cos(qL) y ' 

cos( g )sinh(g) = -- — — V ^ ; . - . (67) 



4 cosh(£ L) — cos(qi) 
We choose a branch of solution qe(A) of (|65| ) and insert it into (^34j). This yields 



L 

Using this result, ( |66| ) and ( |67| ) turn into 

8/ (A) 



COS i 



m- — . (68) 



A = sin ( - ) C oshff ) - - sinfa (A)L/2) 

r L 2 J COS ^ j 4 cos(®(A)L/2) - (-l)-cosh(^) ' [ ' 

( mf _ am) sinh(0 = _^ __ . (to) 

V 2 / vw 4 cosh(fi)- (-l)™cos(#(A)£/2) v 1 



This is a system of two equations in two unknowns £ and A. In contrast to the case N = 2, which was considered 
above, these equations do not decouple. 

We shall first consider the solutions ^, m (A) of equation (|70|). For this purpose we define 

sinh(|L) { qe(*)L \ , U 

■M?) = u7cT\ ' a =(- 1 ) C0S o > b = 7 77TV • (71) 

cosh -a V 2 J 4cos ( m *_«^ 



With these definitions equation (^FJ) turns into 

Binhfc) = &/«(£) , (72) 

Note that a is real and \a\ < 1. Hence, > for all positive £, and a necessary condition for J72] ) to have a 

solution is 6 > 0. This means that 



7T 7T %(A) 37T 7T 

- < m - < — , for [/ > , - < 

2 ~ L 2 ~ 2 ' ' 2 ~ 



7T %(A) 

m 7r 

L 2 



< 7T , for U < . (73) 



Let us concentrate again on the case U > 0. The inequality for the case [7 > in ( |73| ) holds for all real A, if and 
only iff<m-£<^-l (cf. Lemma III.l). Equation @ can be easily solved graphically for fixed a (\a\ < 1) 
and b > (see figure 3). The function / a (£) has the following properties, (i) a < 0: / a (£) is monotonically increasing 
(/£(0 > 0) and concave < 0), /„(0) - 0, #(0) = ^ > § and lim^/^) = 1. (ii) a > 0: /„(£) has a 

single positive maximum and a single positive turning point ^\ < / o (0) = 0, /„(0) = > -j and 
lim^oo f a (0 — 1. These properties are sufficient to conclude that ( |72] ) has a unique solution ^ m (A) for all real A, if 
and only if ^ < m — ^ < 4^ — 1 (recall that U > 0) and the Takahashi condition L > ^ is satisfied. Note that / a (£) 
as a function of a interpolates between the two branches of the function /(£), equation (^). These two branches are 
the two dashed lines envclopping the functions / a (£) in figure 3. 
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H(x> 








a>U . 


sh(x) 


b 






a<0 x^ 7 ^ 






X 



FIG. 3. Sketch of equation 



Lemma III. 1 and Lemma III. 3 allow us to understand the behaviour of the solutions £^ m (A) of ( |70| ) as A — > ±oo. 
Lemma III. 1 applied to ( |70| ) implies 

cos ((m-t+ 1)1) sinh(0 = -^ —^^^L__ , for A — > — oo , (74) 

eos ((m - 0£) sinh(0 = , for A ^ oo . (75) 

Using Lemma III. 3 we conclude that 

lim &, m (A) = Zm-i+i , lim 6,m(A) = £ m -t , (76) 



where £ m is the unique solution of equation (56). This solution exists (cf. ( |58|) and recall that U > 0), if and only if 
■j < m < Hence, the range of validity of© is restricted to 

_ < m _ i < 2_ _ i . (77) 

Let us insert £^ m (A) into fl69|). Then ( ]69| ) turns into 

A = «, m (A) , (78) 

where S^, m (A) is defined by 

^ • / tt <w(A)\ ,/ t ,.u ^ sin(%(A)L/2) 
g<|W (A) = sm [ra T - cosh(fc, m (A)) - T c08(g<(A)L/2) _ ( _ 1)mc0Bh( ^ fn(A)L) • < 79 ) 

Using (|76| ) we obtain the asymptotics of g^ m (A), 

lim ge, m (&) = sin ((to - £ + l)f ) cosh(£ m _ £+:L ) , lim ge, m (A) = sin ((m - cosh(£ TO _^) . (80) 

A— >-oo A— >oo 

We see that <^ jTO (A) has finite asymptotics for A — ► ±oo. Since gi^ m (A) is continuous in A, we arrive at the conclusion 
that there exists a solution A^ m of ([78h for every pair (£, to), which satisfies w7f). Hence, we have shown the following 



Lemma III. 5 Let L > ^, U > 0, N = 3. Equations (j*?j[) and |?^ ) ftave solutions consisting of one fc-A £wo string 
and a single real k only if the momentum q of the center of the string is in the range ^ < q < 4^ . For every choice 
of branch of the real momentum k% (cf. figure 1), there exist L — 1 fc-A two strings for odd L and L — 2 fc-A two 
strings for even L. This gives a total number of L(L — 1) solutions of this type for L odd and of L(L — 2) for L even, 
respectively. 
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Let us compare with Takahashi's counting ( p8[ ) -(p0| ). We have a single free kj and no A string, which means that 
there is one Ij = I and no J™. It follows from (|28|) that I may take L different values. Furthermore, M[ = M' = 1, 
and Mj — for j > 1. Thus |J' | < \{L — 3), which leads to L — 2 possible values of J' 1 . Takahashi's counting 
therefore gives L(L — 2) solutions with one fc-A string and one real fc. This is in aggreement with our Lemma. 
We are now ready to state the following generalization of Lemma III. 4, 

Lemma III. 6 Let L > jj, N = 3. Choose two sequences of integers £l and m^, such that limi_ >00 (mi — = Q> 
f < <Z < Tf- This defines a sequence of solutions &i,mz,(A-) of equation (f?^), which has the limit 

lim & L , mL (A) = - arsinh ( - - ) , (81) 
£-»oo \4cos(g) / 

uniformly in A. Thus all strings corresponding to the sequence £,e L ,m L {A) are driven to their ideal positions. 

Proof: There is no A and no subsequence £f L %mL (A), such that lim^oo ^p L >TOL (A) = 0. This can be seen in similar 
way as in the proof of Lemma III. 4. 

Let us finally note that our considerations for the case N = 3 readily generalize to arbitrary N. For arbitrary N 
we have to consider N — 2 copies of equation (|65|) in the system of equations (|64|)-(|67]). We further have to replace 
<#(A) in j6§| ) by Q(A) = (^(A) + . . . + qg N _ 2 (A) with l\ < . . . < £n-2- The properties of Q(A) (monotonicity and 
asymptotics) then follow from Lemma III. 1 , and all considerations go through as in the case N = 3. 



C. Summary 

In this section we have studied fc-A string solutions of the Lieb-Wu equations (^), (||). We have shown that such 
solutions exist and that they are driven to certain ideal string positions in the limit of a large lattice. We have further 
shown that for a large enough lattice of finite length their number is in accordance with the number of corresponding 
solutions of the discrete Takahashi equations (p5|)-(]27|). 



IV. LIEB-WU EQUATIONS FOR A SINGLE DOWN SPIN (II) - SELF-CONSISTENT SOLUTION 

A. Zeroth order and the discrete Takahashi equations 

In this section we present the self-consistent solution of the Lieb-Wu equations (^) , (0) for the case of three electrons 
and one fc-A string. The Lieb-Wu equations (0), (|j) provide a self-consistent way of calculation of the deviation of 
the strings from their ideal positions. We show that every solution of the discrete Takahashi equations gives an 
approximate solution to the Lieb-Wu equations (0), (@), and we calculate the leading order corrections. These 
corrections vanish exponentially fast as the number of lattice sites L becomes large. 

The Lieb-Wu equations for three electrons and one down spin are 

it t A — sinfc,- — U7/4 „ „ „ 

A — sin kj + ill /4 



3 

Jjl A - sin kj + iU/4 



A - sin kj - iU/4 



(83) 



As in section III. A we may replace equation ( p3| ) by the equation for the conservation of momentum, 

e i(k 1 +k 2 +k 3 )L _ ^ _ rg£\ 



|2p and (g) follow from @ and Q and vice versa. 
Let us follow the usual self-consistent strategy for obtaining a fc-A string solution. As in section III.C we shall use 
the notation fc_ = k\ = q — i£, fc + = q + i£, ^ > 0, i.e. we assume that k\ and k 2 are part of a fc-A string. In order 
to facilitate comparison with the previous literature we shall also use the abbreviations <fi = Re sin(fc + ) = Re sin(fc_) 
and x = I m sin(fc_) = — Im sin(fc + ). We introduce S as a measure of the deviation of the string from its ideal position. 
Then 
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sin(fc_) =<j) + ix = A+ l —+S , sin(k + ) = (f)-ix = A- l —+5 . (85) 



Inserting (ga) into (p3) gives 



^ = i+2? - ^ + = l -Ys • (86) 



We may consider the first equation in ( J85| ) as the equation, that defines S. The second equation in ( pq ) is not 
independent. It is the complex conjugated of the first one and may be dropped for that reason. Similarly, we may 
also drop the second equation in (|S6|). Then we are left with six independent equations, ( |82| ) for j = 3, (|8"i|), and the 
real and imaginary parts of the first equations in (|8^) and (|6|). Note that fci + fc 2 = 2g. Therefore our six equations 
are equivalent to 

2iqL _ A-smfc 3 +kV/4 

A-sinfcs-i^/4 ' [ ' 

e ik 3 L _ A-sinfca -if//4 



A - sin + iJ7/4 ' 
sin(fc_) = A+^-+<5 , (89) 



s = %-r^ ■ 00) 

2 e lfc - L - 1 



Every k-A string solution of the Lieb-Wu equations (|32|), (|83|) gives a solution of equations (|87|)-(|90|) (with real q, £3, 
A and real positive £) and vice versa. 

If L is large and £ = Im k + = 0(1), then 5 is very small and may be neglected in equation (|89|). Then (|9(]) decouples 
from the other equations, which become 



„2ig<°>L _ A -sinfc 3 + i£//4 



A - sin k^'> - iU/4 

^^Ap-sinfcf -iU/4 
A - sin fc 3 + iC//4 

sm(fc W )-A + ^ . (93) 



If, on the other hand, 5 is very small in (|39|), we may neglect it in first approximation and solve (|93j) instead. Now 
( |93| ) implies that ^ < -°- ) = 0(1) (see below). Then, using (^0|), we see that 6 is indeed small for large L. This means the 
assumption 8 be small for large L is self-consistent. 

Let us be more precise. We set k^ — — with real q^ and real, positive Separating ( p3| ) into real 
and imaginary part we obtain two equations which relate the three unknowns q^ ', and Ao, 

A = sin(g(°))cosh(^ 0) ) , (94) 
C (0) - - arsinh ( ^rr) ■ (95) 



4cos(q(°)) 

Let us concentrate on positive coupling U > for simplicity. Since we are assuming that > 0, the range of is 
then restricted to, say, ^ < q^ < 4f by equation (p5|). We obtain the uniform estimate 

£ (0) > arsinh([//4) . (96) 

In order to test, if a solution q(°\ fcg , Aq of ( 91 )- (p3|) is a good approximate solution of (|87])-(90), we use it to 
estimate the modulus of 5 for large L, 

|^| « ^ e -« (0) i < S inh(C//4)L _ (9?) 

2 2 
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The inequality follows fro m j9^ ). We conclude that \S\ becomes very small for large L. For U — 4 and L = 24, 
for instance, the estimate (|97|) gives \S\ < 1.3 • 10~ 9 , whereas the difference between two real k^'s is of the order 
of ^j- = 0.26 (cf. section III.B). For large L it becomes impossible to numerically distinguish between solutions of 
©"© and respectively. 

If we fix the branch of the arcsin as — -| < arcsin(z) < it follows from the inequality ? < < ^ that 

k m = 7T - arcsin (A + f^) . (98) 

Inserting (^8|) into ( |9l|) leads to 

e 2iRe arcsin(Ao+i(7/4)Z- _ Ap - sin fcg ^ — j?7/4 

A - sin 4° 3 + it//4 

We thus have eliminated k^ from the system of equations (|9l])-(|93), and we are left with two independent equations 



2l) and (|9^). These two equations determine the two real unknowns Aq and k^ . Taking logarithms of (92) and (f^) 



3 

we arrive at Takahashi's discrete equatio ns (|2 5| ) , (t27|) for one fc-A string and one real k. 

We have seen that the equations (|9l|)- (p3P~dctcrminc Takahashi's ideal strings. Equations (|j|)-(^8|) on the other 
hand, are equations for non-ideal strings, which solve the Lieb-Wu equations (|82]), (183). Thus, S is a measure for 



the deviation of the strings from their ideal positions. We have further seen that the assumption that S be small is 
self-consistent. In particular, every solution of the equations (|9l|)-(p^), which are equivalent to Takahashi's discrete 
equations, is an approximate solution of equations (87)-(^0|). The approximation becomes extremely accurate for 
large L. 



B. First order corrections 



Inserting any solution of (Pj)-(P3) into (BOJ) we have 



iU 



iU 



(100) 



So the relevant parameter, which controls the deviation of the strings from their ideal positions is e 



Every 

solution of (pl|)-(|93"|) is an approximate solution of (p7|)-(p0|). Let us calculate the leading order corrections. We expect 
them to be proportional to e, 

(101) 
(102) 

(103) 
(104) 



q = 


? (0) 


+ q^e 


+ 0(e 2 ) 


£ = 




K«e 


+ 0(e 2 ) 


fc 3 = 


1.(0) 


+ fc«e 


+ 0{e 2 ) 


A = 


A -+ 


-Aie + 


0{e 2 ) 



The e-expansion of 8 is given in equation (100) above. We also introduce 

4> = sin( g (°))cosh(e (0) ) + 4> {1) e + 0(e 2 ) 

X=^+X {1) e + 0(e 2 ) , 



(105) 
(106) 



since it is the quantities </> = Re sin(fc_) and x = I m sin(fc_) which actually form strings in the complex plane. 

The two sets of variables q, £ and 0, x are n °t independent. Inserting (101) and (102) into the left hand side of 
(105) and (jTOfl) and comparing leading orders in e we find 



)(i) 
,(i) 



3(g(°))cosh(£ (0) ) sin((7( ))sinh(e (0) ) \ / q {1) 
i(g(°))sinh(e (0) ) - cos(g(°))cosh(e (0) ) J \ 



(107) 



Let us insert (100), (105) and (10£) into (J89|) . We obtain to leading order 
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^ =A 1 + ^sm(q^L) 



cos(<? (0) i) 



(108) 
(109) 



i.e. we have already found the leading order correction equation (10S). Next we insert ( |101[ ). ( |103| ) and (104) 
into (p7f) and (Bq) and linearize in e. The resulting equations are 



2g< 



£7 Ax -cos(4° ) )^ 1) 
2i 



* 3 " 2L 



Ao-sin(fcf)) + 

Ax -cos(4 0) )4 1) 
A - S in(4 0) )V + ^ 



LI 

16 



-2<? 



(i) 



(110) 



(111) 



The equation 2q^ + = following from (IK), ( ]l 1 1| ) is, of course, a consequence of momentum conservation. 

Equations (107)-(111) are a system of six linear equations for six unknowns, X ^K an d Aj . Note 

that in the derivation of these equations we have use d all of the equations (^7|)-(p0|). Eq ua tions ( |108] ) and (10E) came 
out of (|89|), ( pp| ) was used to obtain the expansion (100) for <5 in terms of e, and flllO] ), (111) follow from (37) and 
(|88|). Equation (107) is a consequence of the definition of <f> and \- Being a set of linear equations, ( 107 )-( 11_1 ) are 
readily solved, 

A x = -- (sin^L) +tan(q(°))th(^°))cos( g (0) L)) C(L) 



U 



sinOz^L) + tan(gW)th(e (0) ) cos(g (0) L) ) + 0(1/ L) 



2L 
U 



(A -sin(fcf)) + 



cos(fc : 



(0)^ 



C/ 2 (sin(g(°)L) + tan( g (°')th(^°)) cos(g(°>L)) 



4i 



(Ao-sin(fcf)) + 



0{l/L 2 ) 



(i) 



Aj + ^ sin(qWi) = tan(g^)th(e (0) ) cos(g (0) L) + C(l/L) 



[7 



cos(g (0) L) 



(112) 
(113) 
(114) 

(115) 

(116) 
(117) 



The function C(L) in equation (|l!2| ) which gives the explicit form of the 0(1/ L) and 0(1/ L 2 ) corrections in the 
remaining equations is 

-l 

cos( g (°))cosh(£ (0) )(l + tan 2 (g(°))th 2 (£ (0) )) 



C(L) = 



1 



4L 
U 



(A -sin(^ 0) )) + 



2cos(4 0) ) 



1 + G(1/L) 



(118) 



Equations (112 )-( 117 ) give a complete description of the leading order deviation of a non-ideal string from its ideal 
position in the pre sence of on e real fc. The deviations of q and £, which are determined by q^ and follow from 
the equations (107) and ( |116| ), ( |117[ ). In order to see that C(L) — 1 is indeed of order 0(1/ L) on has to use ( |9~i| ) and 



C. Summary 



In this section we have presented a self-consistent solution of the Lieb-Wu equations for the case of three electrons 
and one k-A string. Recall that the existence of these solutions was shown in the previous section. Here we showed 
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that a self-consistent approach naturally leads to Takahashi's discrete equations. We showed that Takahashi's discrete 
equations provide a highly accurate approximate solution of the Lieb-Wu equations in the limit of a large lattice. We 
also showed that there is a natural parameter e = >L that measures the deviation of solutions of the Lieb-Wu 
equations from the corresponding solutions of the discrete Takahashi equations. Employing an algebraic perturbation 
theory we explicitly calculated the leading order deviation in e of a non-ideal fc-A string from its ideal position in the 
presence of a real k. 

V. LIEB-WU EQUATIONS FOR A SINGLE DOWN SPIN (III) - NUMERICAL SOLUTION 

A. Numerical method 

In section III the existence of fc-A string solutions was analytically proven under the Takahashi condition for the 
simplest non-trivial cases with one down spin, M = 1. The deviation of ideal string solutions given by the discrete 
Takahashi equations from the corresponding solutions of the Lieb-Wu equations was evaluated analytically in section 
IV, and it was shown that the corrections vanish exponentially fast as the lattice size L becomes large. To confirm 
these analytical arguments on the existence of fc-A strings, we utilize a complementary numerical approach. Although 
the tractable system size is limited, we can directly obtain fc-A strings and verify the completeness of the Bethe ansatz 
for arbitrary U. We note that numerical solutions for low-lying particle-hole and A string excitations in small finite 
systems can be found in the literature (see e.g. p7[|). As far as we know, however, our data are the first example of a 
numerical study that confirms the completeness of the Bethe ansatz for a small finite system. 

We study again the cases N = 2 and N = 3 for both pedagogical clarity and technical simplicity. We shall employ 
two numerically exact methods, (1) the numerical diagonalization of a real-symmetric matrix using the Householder- 
QR method (we call it Method 1), (2) a numerical method to solve coupled nonlinear equations using the Brent 
method (we call it Method 2). These techniques themselves are conventional. They allow us to obtain numerically 
exact solutions for the Hubbard model. Here we use the term 'numerically exact' to state that our numerical solutions 
give exact numbers except for the inevitable rounding error in the computation.^ 

Our strategy here is the following. 

(i) Obtain all energy eigenvalues with fixed N and M as a function of U using Method 1. This step gives a complete 
list of energy eigenvalues. 

(ii) Obtain numerical (real and/or complex) solutions by solving the Lieb-Wu equations with Method 2. 

(iii) List up all eigenvalues obtained by the two methods and compare them with one another. This step gives a 
confirmation of completeness. 

For our numerical study we used the following form of the Hubbard Hamiltonian, 

L L 

a ( c t+i,« c i,° + c t-i,* c w) + '»<•••"<.. ■ ( 119 ) 

3=1 <r=U 3 = 1 

This form is different from ([!]) by a shift of the chemical potential and by a constant energy shift. For fixed particle 
number N this leads to a shift of the spectrum by ^(2N — L). 

In order to perform a numerical diagonalization we used basis vectors c+ x . . . ct N - M c yi ■ ■ ■ c y M |0) to represent the 
Hamiltonian as a matrix in the sector of fixed N and M. The number of different configurations, x± < . . . < xn-m, 
yi < ... < dm, in this sector is ( N ^ M ) X (^) and determines the size of the matrix. In order to employ method 
2, we rewrote the Lieb-Wu equations into a proper set of real equations, which will be presented in the following 
subsections. Hereafter in this section we assume that L is an even integer. 



3 We will indicate errors by the difference of the left hand side and the right side of each equation in (M), (^), evaluated within 
our numerical treatment. Then it will become clear that in all the equations which we tested the relative error is negligible and 
of the order of the rounding error expected for the double- precision calculation. 
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B. Numerical solution for two electrons 



1. Equations for the k-A string 

Let us discuss numerical k-A two string solutions for the two electron system with one down spin (N = 2, M = 1). 
We shall use the same notation as in section III. The variables to be determined are k + , fc_ and A, where k + is the 
complex conjugate of fc_ and A is real. We may write 

k + = q + iC , k_=q-i£ , (120) 

where < q < 2tt and £ > 0. Then the total momentum is 

P = 2q = mod 2tt . (121) 

L 

This equation restricts the admissible values of q. It is easy to obtain A as a function of q and £. We find A = 
sin(g) cosh(£) (cf. equation (|55|)). Thus we are left with a single equation that determines £. For our numerical 
calculations we wrote it in the form 

— cos(m7r/L) sinh(f) + U/4 , . 
exp (imir + £&) = '[ . - fj -f- . 122 

- cos(m7r/Z)sinh(£) - U/4 



This equation is equivalent to (|56|). Therefore, for positive U, the allowed values of m axe restricted to 

m = £/2 + l,...,3L/2-l (123) 



(cf. equation (J5^) ) . Equation (|l22j) was already studied in appendix B of |60[ . There it was shown that there is a 
redistribution phenomenon as U becomes small. k-A strings corresponding to odd values of m collapse at critical 
values of U given by U m — (8/L)\ cos(m7r/L)|. 



2. Numerical solutions for N = 2 

As a typical example, let us present some numerical fc-A two string solutions for N — 2 (two electrons). For a fc-A 
two string we show the dependence of energy eigenvalues, imaginary parts of charge momenta and the deviation from 
the ideal string positions on U. We put rn = 16 and L = 16 (16 sites). Then we have q = it, i.e. 

k + = 7T + i£ , fc_ = 7T - i£ , A = . 

In the list below the deviation of the string from its ideal position, }^ eal ^ = n j- [ arsinh([7 /4) is measured by 
lmS = sinh(£) - U/4. 

(1) U = 10 

Householder E 
Bethe Ansatz E 

Im<5 

(2) U = 1 

Householder E 
Bethe Ansatz E 

Im<5 



= 10.7703296143355 
= 10.7703296143355 
= 1.64723114637774 
= 1.78994596922166 x 10~ n 

= 4.13148449882288 

= 4.13148449882289 

= 0.255705305537532 

= 8.50098694369150 x 10~ 3 



24 



(3) [7 = 0.1 

Householder E = 4.00668762410970 
BetheAnsatz E = 4.00668762410971 

£ = 5.78176505356310 x 10~ 2 
ImS = 3.28498688384022 x 10~ 2 

(4) U = 0.01 

Householder E = 4.00062917225929 
BetheAnsatz E = 4.00062917225931 

£ = 1.77363435624506 x 10~ 2 
1xo.5 = 1.52372734873120 x 10~ 2 

The string with m = 16 does exist for any U > and is actually the highest level in the spectrum for N = 2. We 
note that the non-ideal string approaches the ideal fc-A string, when U becomes large, even for such a small system. 
Yet, in accordance with our expectations, the ideal string does not provide a good approximation for small U in finite 
systems. 

In table 1 we show a complete list of eigenstates for the case of one down-spin (M = 1, S x = 0) and U — 1.5 for a 
6-site system (L — 6). Note that the value of U — 1.5 is greater than maxU m = 1.1547, or, in the language of section 
III, the Takahashi condition is satisfied. Let us explain the table. 

(i) The 36 eigenstates are listed in increasing order with respect to their energy. 

(ii) S and P denote the spin and momentum of the eigenstate, respectively. 

(iii) The energy eigenvalues obtained by direct diagonalization of the Hamiltonian (the Householder method) and 
by the Bethe ansatz method coincide within an error of O(10 -15 ). 

(iv) The last digit for each numerical value has a rounding error. 

(v) There are 5 fc-A string solutions among the 36 eigenstates, which is consistent with the number L — 1 = 5 
obtained by Takahashi's counting (p8|)-(|30|). 

Let us give more explanations on the table. In fact, it confirms the completeness of the Bethe ansatz as discussed 
in section II. D. We first note that there are 36 eigenstates for the case of two electrons and one down-spin (N = 2 
and M — 1) on a 6-site lattice (L = 6). We recall that the two electrons with one up-spin and one down-spin can 
occupy the same site. The (f) x ft) = 36 eigenstates in table 1 can be classified into the following types. 

(i) 15 eigenstates with two real charge momenta k%, k 2 and one real spin-rapidity A. 

(ii) 5 eigenstates with one fc-A two string. 

(iii) 15 eigenstates with S = 1 and S z — belonging to spin triplets. 

(iv) 1 ^-pairing state. 

Let us consider case (i). In table 1 there are 15 states with real charge momenta which have total spin zero (5 = 0). 
This agrees with our analytical arguments in section III on the number of real, regular Bethe ansatz solutions (cf. 
Lemma III. 2 and below), which should be = 15. Let us also recall that the Ij are integer (half-integer) when 
V M m + M m is even (odd). Thus, for the case M = 1, the Ij should be half-odd integer, which is in agreement 
with the results shown in table 1. 

We now consider case (ii). From the analytic discussion in section III.C, we should have 5 eigenstates with fc-A 
strings. This is in accordance with our numerical data in table 1. Recall that we showed in section III.C (below Lemma 



III. 3) that also Takahashi's counting, using equations (28)-(p0|), leads to the same number of fc-A string solutions 
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No. 


Energy 


S 


P/(tt/3) 


type of solution 


1 


-3.82047006625301 








real Ij = -0.5, 0.5 


2 


-3.00000000000000 


1 


1 


triplet Ij = 0, 1 


3 


-3.00000000000000 


1 


5 


triplet Ij =0,-1 


4 


-2.60959865138515 





1 


real Ij = -0.5, 1.5 


5 


-2.60959865138515 





5 


real Ij = -1.5,0.5 


6 


-2.00000000000000 


1 





triplet Ij = 1,-1 


7 


-1.86256622153075 





2 


real Ij = 0.5,1.5 


8 


-1.86256622153075 





4 


real Ij = —0.5, —1.5 


9 


-1.53909577258373 








real Ij — —1.5, 1.5 


10 


-1.00000000000000 


1 


2 


triplet Ij =0,2 


11 


-1.00000000000000 


1 


4 


triplet Ij =0,-2 


12 


-0.594210897273940 





2 


real Ij = -0.5, 2.5 


13 


-0.594210897273940 





4 


real Ij = -2.5, 0.5 


14 


0.00000000000000 





3 


real Ij = 0.5, 2.5 


15 


0.00000000000000 





3 


real Ij = -2.5, -0.5 


16 


0.00000000000000 


1 


3 


triplet Ij = 1,2 


17 


0.00000000000000 


1 


3 


triplet Ij = -1,-2 


18 


0.00000000000000 


1 


3 


triplet I 3 ; = 0, 3 


19 


0.00000000000000 


1 


1 


triplet = 2, -1 


20 


0.00000000000000 


1 


5 


triplet = -2, 1 


21 


0.474357244982949 





1 


real = -1.5,2.5 


22 


0.474357244982949 





5 


real I 3 = -2.5, 1.5 


23 


1.00000000000000 


1 


2 


triplet Ij = 3,-1 


24 


1.00000000000000 


1 


4 


triplet Ij =3,1 


25 


1.43079477929458 





4 


real Ij = 1.5,2.5 


26 


1.43079477929458 





2 


real Jj = -2.5,-1.5 


27 


1.50000000000000 





3 


77-pair 


28 


2.00000000000000 


1 





triplet Ij = 2,-2 


29 


2.49213401737360 








real Ij = -2.5, 2.5 


30 


2.52598233951011 





2 


complex to = 8, £ = 0.710224864788777 


31 


2.52598233951011 





4 


complex to = 4, £ = 0.710224864788777 


32 


3.00000000000000 


1 


5 


triplet Ij =2,3 


33 


3.00000000000000 


1 


1 


triplet Ij = —2,3 


34 


3.63524140640220 





1 


complex to = 7, £ = 0.313056827256169 


35 


3.63524140640220 





5 


complex to = 5, £ = 0.313056827256169 


36 


4.36743182146314 








complex to = 6, £ = 0.425405934759021 



Table 1: Classification of all energy levels for L = 6, N = 2, M = 1 and U = 1.5. 
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We consider case (iii). There are 15 states with 5 = 1 and S z = 0. We describe them as triplets in table 1. They are 
obtained by multiplying the spin-lowering operator to the regular Bethe states with 5 = 1 and S z = 1. (For the 
notation for the £0(4) symmetry see section II.B). The regular Bethe states with 5 = 1 and S z = 1 correspond to 
N = 2 and M = 0. We recall again that the Ij's are integer (half- integer) when ^ m M m + M m is even (odd). Thus, 
the Ij's which belong to the S z = states in spin triplets should be integer- valued. They should take one of the 
values —2, — 1, 0, 1, 2, or 3, which means that there are (^j = 15 states according to Takahashi's counting (p8|)-(p0|). 

There is one ^-pairing state. The energy of this state is equal to U, since the two electrons occupy the same site. 

Now, let us sum up all the numbers of the different types of eigenstates: 

15 + 5 + 15 + 1 = 36 . (124) 

Thus, we have shown that all the energy eigenstates obtained by Method 1 are confirmed by Method 2. In particular, 
we have confirmed numerically the completeness of the Bethe ansatz. 

Let us consider the total momentum. Using equations (ftq)-(|27l) we can express the total momentum P of the 
eigenstates with real charge momenta in terms of I\ and I2, 

P = ^-{h + I 2 ) mod27r. (125) 

This formula is consistent with table 1. 

Let us now discuss the [/-dependence of the spectrum. In figure 4 we show the spectral flow from strong-coupling 
to weak-coupling, where the Takahashi condition does not hold. 



10.0 




-5.0 1 1 

0.0 1.0 2.0 3.0 4.0 5.0 

U 

FIG. 4. The spectral flows for N = 2 and M = 1 for a 6-site lattice. Solid lines denote the k-A strings. Dashed lines denote 
the energy of real roots (5 = 0). Dotted lines denote the triplet states (5 = 1). The dash-dotted line denotes the energy of the 
ry-pair. At U5 = U7 = 1.1547, as indicated by the arrow, two fc-A two strings with m — 5, 7 collapse into real solutions. 



In figure 4 we show the redistribution phenomenon discussed in sections II. D, III. A and above. There are five fc-A 
two-string solutions in table 1. The entries 30, 31 and 36 correspond to even m. According to our discussion above, 
these k-A two strings are stable as U becomes small. Entry number 36 is the highest energy level in the figure. Entries 
number 30 and 31 are degenerate. They correspond to the third highest level at U = 5. The entries number 34 and 
35 correspond to odd m = 5,7. The corresponding levels are degenerate. At U5 = U7 = 1.1547 these fc-A two-string 
solutions collapse into pairs of real charge momenta. This is indicated by the arrow in figure 4. On the other hand, 
all five fc-A string solutions do exist as long as the Takahashi condition is satisfied. 
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C. Numerical solution for three electrons 



1. Equations for the k-A string 

We shall now consider the set of solutions to the Lieb-Wu equations for N = 3. The fc-A strings to be searched for 
have two complex kj, a real kj and a real A. We express these momenta and rapidities by 

k 1 =q-i£ , fc 2 =<z + i£ , k 3 , A . (126) 

The total momentum is 2irm/L, so that 

fci + k 2 + k 3 = 2vrm/i , 2tt(£ -1)/L<k 3 < 2n£/L . (127) 

We define Re 6 and Im 5 by 

sin ki = A + iU/4 + Re 6 + i Im 6 , (128) 

smfc 2 = A-i£//4 + Re£-iIm<5 • (129) 

We recall that Re<5 and lm5 describe the deviations from the ideal string solutions. Now we derive equations for 
four variables, k 3 , q : £ and A, starting from the Lieb-Wu equations (|82|), (|S3|). Taking logarithms, we obtain a set of 
equations of the form /j = (i = 1, . . .4). Here the fa's are defined by 

h =Lk 3 -2 arctan ( ^ZjE^ j _ 2 n(£ - 1/2) , (130) 



f 2 = -2iarcsin (^1/2^ (A + Re 6 + l) 2 + (U/4 + ImS) 2 

- V(A + RcS - l) 2 + (U/A + Ini(5) 2 } J + L/c 3 - 2tt(£ - 1/2 - J') , (131) 



A = ^ ( - 2i0 - 75 iS^ta^ , (132) 



(Re J) 2 + (U/2 + lmS) 2 
25\ ( Re 



. ' Re 5 \ / Re <5 \ / A - sin fc 3 , 
j4 = —2 arctan - — - I + 2 arctan — — - — - — - I — 2 arctan — — I — air , (133) 



where 



TLeS = singcosh£ — A , Im<5 = — U/A — cosgsinh^ , (134) 
and the parameter a is given by 

o- = l for k 3 <^(£-l/2) , <r=-l for k 3 >^-(£-l/2) . (135) 

Note that the equation fi = is equivalent to equation (Q) for j = 3. 

To solve this set of coupled equations by Method 2, we need a proper initial guess. We employ the ideal strings 
given by the discrete Takahashi equations as initial approximation. The fact that the ideal strings provide a good 
estimate for the true solution is crucial to reach the correct answer. This is because the /j's are rather singular 
functions having many diverging points. Very often, the true solution is very close to a divergent point. We cannot 
approach the solution from a point beyond a branch cut using an iterative way like the Brent method. 



2. Numerical examples of k-A string solutions 



We present some numerical solutions for one fc-A two string and one real k. The parameters in the examples below 
are N = 3 (three electrons), L = 10 (10 sites) and {7 = 5. 
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(1) m = 10, £ = 1 (7 3 = 0, J = 1/2) 

Householder 
Bethc Ansatz 



Q 

A 
Re <5 
Im<5 



£ = 4.46666961980768 
E = 4.46666961980768 
2.98893848049280 
1.05674954466496 
0.305308346193987 
0.245232889885225 
-6.42851441900183 x 10" 
2.84511548476196 x 10~ 6 



(2) m = 7, 1 = 1 (7 3 = 0, J = 7/2) 

Householder £ = 3.48250616148668 

Bethe ansatz E = 3.48250616148668 

q = 1.92454676428806 

£ = 1.99506896270533 

k 3 = 0.549136186449599 

A = 3.51250692222577 

Re5 = 2.08765360554253 x 10~ 9 

ImS = 4.99459629210719 x 10~ 9 

Let us discuss possible numerical errors for the above solutions. Their numerical errors may be evaluated by the 
residual, f it given for these solutions as follows. 



(1) m = 10, £ = 1 (7 3 = 0, J = 1/2) 



with a = +1. 

(2) m = 7, t = 1 (7 3 = 0, J = 7/2) 



h 
h 
H 
U 



fi 
h 

h 
U 



-1.24900090270330 x 10~ 15 
-3.99680288865056 x 10~ 15 
5.25259688971173 x 10" 22 



4.53592718940854 x 10 - 



-12 



-1.33226762955019 x 10~ 15 
1.42108547152020 x 10" 15 
-4.66698631842825 x 10~ 25 
-7.67983898697366 x 10" 10 



with a = —1. 

In comparison to other error values /, the number /j has a rather large value. However, in f±, we have an expression 
like e/e' (e ~ 0, e' ~ 0). So this value contains a larger cancellation error for smaller Re<5 and ImS. Note again that 
relative errors for the energy are always 0(1O~ 15 ). 
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3. A complete list of solutions for N = 3 and M = 1 

We can present complete lists of eigenstates for all finite systems tractable by our numerical technique. As a further 
example, we consider all eigenstates for N = 3 and M = 1 in a 6-site system (L = 6). The list is shown in appendix C. 
It confirms again the completeness of the Bethe ansatz. 

Let us briefly discuss the numbers of eigenstates of different types. First, we note that there are in total 6(2) = 90 
eigenstates. Inspection of the tables in appendix C shows that they can be classified into the following four types. 

(i) 40 eigenstates with three real charge momenta fci, fc2, and one real spin rapidity A. 

(ii) 24 eigenstates with one k-A string with k\ = fc_ = q — i£, k^ = fc + — q + i£, and A and k% real. 

(iii) 20 eigenstates with S — 3/2 and S z = 1/2 belonging to spin quartets. 

(iv) 6 eigenstates with one 77-pair and one real charge momentum. 

Let us now confirm that these numbers agree with Takahashi's counting, (pq)-(pO|): Case (i) was considered in 
section III. A below Lemma III. 2. There we showed that Takahashi's counting predicts a number of 2 (3) real solutions 
for three electrons and one down spin. Inserting L = 6 we have 2(g) = 40 eigenstates, which is in accordance with 
our numerical calculation. Case (ii) was considered below Lemma III. 5. The number of eigenstates obtained there by 
Takahashi's counting was L(L — 2), which for L = 6 gives as desired 6 • 4 = 24. Let us consider case (iii). These states 
are the second highest states (S z = 1/2) in spin quartets. They arc obtained from regular Bethe states with N = 3, 
M = by multiplication with the spin lowering operator £t. Since N = 3 and M = 0, we have u) = 20 states of this 
type. The ?7-pair in case (iv) is obtained by acting with rj + on regular Bethe states with N = 1 and M = 0. Hence, 
Takahashi's counting gives 6 eigenstates of this type on a 6-site lattice. 



D. Summary 



In this section we have presented a thorough numerical study of the Hubbard model. We calculated, in particular, 
all eigenstates and eigenvectors for a six-site lattice with two and three electrons and one down spin by direct numerical 
diagonalization of the Hamiltonian. These data were compared with data obtained by numerical solution of the Lieb- 
Wu equations. Both sets of data are in perfect numerical agreement and confirm once again the results of our analytical 
investigation in the previous sections. The structure of our numerical data is fully consistent with Takahashi's string 
hypothesis. The number and classification of the eigenstates is consistent with Takahashi's counting (p8|)-(^o[). Thus, 
our numerical data confirm not only the existence of fc-A strings but also the completeness of the Bethe ansatz. 
Without fc-A strings the Bethe ansatz would be incomplete. @ 



VI. THERMODYNAMICS IN THE YANG- YANG APPROACH AND EXCITATION SPECTRUM IN THE 

INFINITE VOLUME 



Let us now turn to the determination of thermodynamic quantities and the zero-temperature excitation spectrum 
in the infinite volume. A convenient way to construct the spectrum was pioneered by C.N. Yang and CP. Yang 
for the case of the delta-function Bose gas [Q. The starting point are the Bethe Ansatz equations in the finite 
volume. They are used to derive a set of coupled, nonlinear integral equations called thermodynamic Bethe Ansatz 
(TBA) equations, which describe the thermodynamics of the model at finite temperatures. The quantities entering 
these equations have a natural interpretation in terms of dressed energies of elementary excitations. Yang and Yang's 
formalism is a natural generalization of the thermodynamics of the free Fermi gas to interacting systems. 

In what follows we review Takahashi's derivation of the TBA equations for the case of the repulsive Hubbard model 
fUl . The analogous calculations for the attractive case can be found in . 

Our starting point are the discrete Takahashi equations (^B|)-f|2^) and expressions for energy (|3^) and momentum 
( |32f) for very large but finite L. A very important property of (P5f)-(p^) is that as we approach the thermodynamic 
limit L — > 00, N/L and M/L fixed (finite densities of electrons and spin down electrons), the roots of (|25|)-(27) become 
dense 



k j+1 - kj = OiL- 1 ), A n a+1 -A n a = 0(L- 1 ), A' n a+1 -A'l = 0{L- 1 ). 



(136) 
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We now define so-called counting functions y, z n , z' n as follows 



n—1 a—1 



j=l v ' 7 m=l /3=1 v ' 7 

L[arcsin(A' + nf) + arcsin(A'-nf)]=Lz; i (A')+ E g J!"' +EE 6 " ^T 1 ) ■ (139) 

j=l ^ ' 7 ro=l/3=l \ / / 

By definition the counting functions satisfy the following equations when evaluated for a given solution of the discrete 
Takahashi equations 

y(k j )=2nI j /L, z' n (A' n a ) = 2-kJ'I/L , z n {A n a ) = 2^/ L . (140) 

In the next step we define the so-called root densities, which are related to the counting functions as follows. By 
definition the counting functions "enumerate" the Bethe Ansatz roots e.g. 

L[y(k,) - y(k n )} = 27T& - I„). (141) 



For a given solution of (25)-(g7j) certain of the (half-odd) integers between Ij and /„ will be "occupied" i.e. there will 
be a corresponding root k, whereas others will be omitted. We describe the corresponding fc-values in terms of a root 
density p(k) for "particles" and a density p h (k) for "holes". In a very large system we then have by definition (here 
the property ( 136 ) is very important) 

Lp(k) dk — number of k's in dk , 
Lp h (k) dk = number of holes in dk . (142) 
It is then clear that in the thermodynamic limit we have 

2n[p{k) + p h (k)} = ^-. (143) 



The analogous equations for the other roots of (p5|)-(|27D are 

2n[a n (A) + o*(A)] = , JhrKCA) + = • (144) 

In the thermodynamic limit the discrete Takahashi equations can now be turned into coupled integral equations 
involving both counting functions and root densities 

k = y(k) -YlJ J Ad { ZuJa J K(A) + CT " (A)1 ' (M5) 

fjk 9 ( A ~J^ fc ) P(k) = z n (A) + J2 j™JA> B nm (j^jf) v m (A'), (146) 

arcsin(A + n^) + arcsin(A - n^-) = z' n (A) + J dk 9 (J~^Jj^J P( k ) 

+ £ / dA ' Q «™ \~UJ4~) CT ™ (A,) - (W7) 



m— 1 



As we are interested in the Hubbard model at finite temperatures we need to express the entropy in terms of the 
root densities. This is achieved by observing that e.g. the number of vacancies for fc's in the interval [k, k + dk] is 
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simply L(p(k) + p h (k))dk. Of these "vacancies" Lp(k)dk are occupied. The corresponding contribution dS to the 
entropy is thus formally 



ds _ (L[p{k)+p h {k)]dk)\ 
{Lp{k)dk)\{Lp h {k)dk)\ 



(148) 



where ! denotes the factorial. The differential dS is obtained via Stirling's formula. After integration we obtain the 
following expression for the total entropy density of the Hubbard model 



S/L= / dk{[p(k)+p h (k)}ln[p(k)+p h (k)}-p(k)lnp(k)-p h (k)lnp h (k)} 

J — 7T 

+ dA{[a n (A) + a' r l (A)] ln[<7„ (A) + a h n (A)] - a n (A) In a n (A) - a h n (A) In a\\ (A) } 

{K(A) + a' n h (A)\ lnK(A) + ^(A)] - ^(A) ln<(A) - a' n h (A) \na' n h (A)} 



E 



dA 



The Gibbs free energy per site is thus 

E - pN — 2BS Z - TS 



f = 



/7T 00 POO 

dk [-2 cos k - ^ - C//2 - S] p(fc) + E / dA 2nS CTr 



(A) 



E 



dA 



4ReVl - (A - inU/4) 2 - 2np - nU a' n (A) - T S/L . 



(149) 



(150) 



Here p is a chemical potential, B is a magnetic field and T is the temperature. The alert reader will have realized 
that we are still missing a set of equ atio ns tha t allows us to completely determine the root densities and counting 
functions (the thermodynamic limit ( |l45| )-(147) of the discrete Takahashi equations are clearly insufficient). This is 
the topic of the following subsection. 



A. Takahashi's Thermodynamic Equations 

Let us start by differentiating ( |!45[ )-( pT7| ), which yields 

nU/4 <(A) + a n {A) 



1 f°° nU 

p(k) + p h (k) = — +cosfcy^ / dA '■ , . 

KW ' w 2tt _ x J-az ti" (™£//4) 2 + (smfc - A) 2 

cr£(A) = - A nm * cr m + / dfc !— 

m=l A • / -' r 7r 



p(fc) 



<"(A) = ±Re 



1 



(n[//4) 2 + (sinfc-A) 2 ' 
nt//4 



5> 



tt ^1 - (A - in£7/4) 5 ,„_, 
Here A nm is an integral operator acting on a function / as 

d 



dk- 



p(k) 



7T (nt//4) 2 + (sinfc - A) 5 



5nmf{x) + 



dx 



dy n ( x-y 
2tt nm \ U/A 



(151) 



(152) 



Equations (151) can be used to express the densities of holes in terms of densities of pa rticles. 
A second set of equations is obtained by considering the Gibbs free energy density (150) as a functional of the root 
densities. In thermal equilibrium / is stationary with respect to variations of the root densities 
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= Sf 



JL 

5p(k)' 



(153) 



where we need to take into account (151) as constraint equations. In this way one obtains a set of equations for the 
ratios C = p h /p, Vn = o"£/<7„ and rf n = a' n ja' n 



lnC(fc) 



-2cosfc - fx - [7/2 - B 



„=l J -o° 



1 



4tt (n£//4) 2 + (sinfc-A) 2 



In I H ttt \ - hi f 1 ' 



<(A) 



»7n(A) 



(154) 



ln(l + 7? n (A))+ / dfc 



cosfc 



nU/A 



it (n?7/4) 2 + (sinfc-A) 2 



In 1 



C(k) 



2nB 



oo / 1 

A„ m * In I H 



(155) 



ln(l + 7j;(A)) + 



dk- 



cos k 



nU/A 



7T (n?7/4) 2 + (sinfc - A) 2 



In 1 



1 



4Rev/l - (A - inU/A) 2 - 2n/x - nC7 



In 1 



m= 1 



C(fc) 
1 



(156) 



Note that (151) together with (|1 5 4[) - (|156|) completely determine the densities of holes and particles in the state of 



thermal equilibrium. 



The Gibbs free energy per site is given in terms of solutions of (154)-(156) as 



(IA 



In 1 



1 



Re- 



« v <(A)7 y/l-(A- inU/A) 2 
Following Takahashi we define 

n(k) = Tln(C(*)) .en(A) = Tm(7j„(A)) .^(A) = Tln(^(A)) . 



(157) 



(158) 



As was first shown by Yang and Yang for the delta-function Bose gas |5(J , the quantities defined in this way describe 
the dressed energies of elementary excitations in the zero temperature limit. Before turning to this we will give a brief 
summary on how to calculate thermodynamic quantities in the framework of Takahashi's approach. 



B. Thermodynamics 



The expression for the Gibbs free energy density (157) can be simplified |22| 

poo 

dk Po (k) ln(l + C(fc))+ / dA ctq(A) In (1 + 771(A)) 



where 



f = E -Li- U/2 - T 



to (A) 



dk- 



1 



1 



U cosh f- (A -sinfc) 



1 , f°° dA 

PoW = — + cosfc / 

2tt J_ 00 tt 



Po(k) 
U/A 



En = -2 



(J7/4) 2 + (sinfc-A) 2 
dk cos(fc) po(k) = —4 



<ro(A) 



00 to. WW 



o 1 + exp(vU/2) 



(159) 



(160) 
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We note that po, ao and Eq are the root density for real fc's, the root density for real A's and the ground state energy 
for the half-filled repulsive Hubbard model, respectively. Since the occurrence of quantities related to the half-filled 



Hubbard model in (159) may be surprising, we would like to emphasize that (159), (160) holds for all negative values 
of the chemical potential fi i.e. for all particle densities between zero and one. 



The representation (159) is convenient as it shows that the Gibbs free energy is determined by the dressed energies 
for real fc's and real A's only. In order to derive (|159|) the following identities are useful 



dk 



dk 



nU 



Po(fc) 



1 



-Re- 



1 



4tt (nC//4) 2 + (sinfc-A) 2 n \/\ _ (A - inU/A) 2 
nU po(k) 



An (nC//4) 2 + (sinfc-A) 2 



A n i * (Tq 



(161) 



At very low temperatures T <C B it is possible to determine the Gibbs free energy by using an expansion of the 
TBA equations ( |l 5 4| ) - ( J]~56| ) for small T ]2^]. The TBA equations essentially reduce two only two coupled equations 
for £ and 771 in this limit. For generic values of B and arbitrary temperatures one needs to resort to a numerical 
solution of ( [15 4j ) - ( |l~56| ) . In order to do so, one needs to truncate the infinite towers of equations for A and fc-A strings 
at some finite value of their respective lengths. In |6^j6^] such a truncated set of equations was solved by iteration. 
The integrals were discretized by using of the order of 50 (100) points for the fc (A) integrations. The results of these 
computations are compared to the results of the Quantum Transfer Matrix approach in the next section. 



C. Zero Temperature Limit 



For the remainder of this section we set t he m a gneti c field to zero B = 0. 
For T — > the thermodynamic equations (154)-(|l56|) then reduce to [^8) 



f°° dA 

K (k) = -2cosfc - \i - U/2+ / 

J-00 n 



U/A 



(t//4) 2 + (sinfc-A) : 



rei(A) , 



(162) 



(l-5 nil )e n (A) 



^ cosfc 
dk 



nU/A 



it {nU/A) 2 + (sinfc - Af 



n(k) - A n i * ei 



(163) 



4(A) = 4ReVl-(A-znC//4) 2 - 2np - nU + t dk C ° Sk ^ «(fc) . (164) 

J-Q 7r {nU + (smfc — Ay 



Note that k(±Q) = 0, re(fc) < for fc| < Q and ei(A) < 0. Using Fourier transform, equations (|l62j) and ( |l63| ) can 
be simplified further with the result 

rQ 

K {k) = -2 cos k - [i - U/2 + / dk' cos fc' i?(sin fc' - sin fc) n(k') , 

J-Q 



Q 



cos fc 



ei(A) - / dk— -3— — k(^) , 

7_q U cosh jf- (A - smfc) 

e„(A)=0 n = 2,3,... (165) 

where 

= du exp^x) 

1 ; 2^1 + cxp(C/| W |/2) 1 ' 

The vanishing of the dressed energies of A-strings of lengths greater than one i.e. e n (A) = for n > 2 is due to the 
absence of a magnetic field. For finite magnetic fields all e n (A) will be nontrivial functions. 

In order to characterize the excitation spectrum we need to determine the dressed momenta in addition to the 
dressed energies. This can be done by considering the zero temperature limit or (151) 
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p{k) 



i r Q 

h / dk'cosk R(smk' - sinfc) p(fc') , Ifel < Q 

27r 7-Q 



'-<2 



P h ( k ) = 7T + / cos fc i?(sin k' - sin fc) p(fc') , Ifel > Q 

2tt 7_ q 



= , 
P h (k)=0 



\k\>Q, 
\k\ < Q 



CTl (A) 



dk 



-Q 
1 



_q [/ cosh ^ (A - sin fc) 



P(A) 



<r„(A) = n = 2,3,. 



1 



(A) = = <4(A) m = l,2, 



< n (A) = -Re 



7T ^1 - (A - mf7/4) 5 



Q „ n£//4 
dk 



p(k) 



7T (nJ7/4) 2 + (sinfc-A) 2 



(167) 



Equations (167) describe the ground state of the repulsive Hubbard model at zero temperature and zero magnetic 
field. There is one Fermi sea of fc's (charge degrees of freedom) with Fermi rapidity ±Q and a second Fermi sea of 
A^s (spin degrees of freedom), which arc filled on the entire real axis. 

The total momentum (p3) can be rewritten by using (|25|)-([27|) in the following useful manner 



P 



2tt 



' N—2M' oo Mm oo M' n 

j=l m=l p=l 



n—1 a— 1 

oo ML 



N-2M' oo M m oo M; oo M' n 

E vfo) + E E z ^f) -EE <( A «) + ^EE^ + 1 )- 

j=l m=l /9=1 



(168) 



n— 1 a— 1 



n—1 a— 1 



Using this expression for the total momentum we now can identify the dressed momenta of various types of excitations. 
We find that an additional real k with |A;| > Q ("particle") or hole in the sea of fc's (\k\ < Q) carry momentum ±p(fc) 
respectively, where 

(169) 



p(fc) = y(k) = 2tt / dk' [p(k') + p h (k')} . 
Jo 

Similarly, the dressed momentum of a hole in the sea of A 1 's is 

/•oo f>Q / \ pj 

Pi(A) = --Zi(A) = 2n dA'cri(A') - Zi(oo) = 2 / dk arctan[exp ( (A - sinfc) )] p(k) - tt — . 

Ja J-q V U J 2L 

This result was first obtained by Coll [^8) . Finally, a fc-A string of length n has dressed momentum 



(170) 



-2Rearcsin(A - inU/A) + 



dk 2 arctan 



A — sin fc 
nU/4 



p(k) +n(n + 1), 



(171) 



where the second line is obtained from the T — > limit of (147). We are now in a position to completely class if y the 
excitation spectrum at zero temperature. The dispersion relations of all elementary excitations follow from (164)-( |l65| ) 
and ( 169 )-( 171 ). These equations involve only the two unknown functions, p{k) and Ac(fc), which are solutions of linear 
Fredholm integral equations. 



D. Ground state for a less than half- filled band 



The integral equations describing the root densities of the ground state are (167) 



i r Q 

p{k) = h / dk cos k R(sin k - sin k) p(k') , \k\ 

27r J-Q 



< 



ax(A) 



dk- 



Q U cosh (A - sinfc) 



:P(k) 



(172) 



35 



where 

/Q poo 

dkp{k) = N/L , / dAcri(A) = M\/L = N/2L . (173) 

-Q J-oo 

The ground state energy per site is given by 

rQ 

(E-iiN)/L= dk (-2cosfc-/i-[//2) p(k) . (174) 
J-Q 

E. Excitations for a less than half-filled band 

There are three different kinds of elementary excitations. 

• The first type of elementary excitation is gapless and involves the charge degrees of freedom only. It corresponds 
to adding a particle to or making a hole in the distribution of fc's. Such excitations have dressed energy ^n(k) 
and dressed momentum p J , : / l (fc). 

• The second type of elementary excitation is gapless and carries spin but no charge. It corresponds to a hole 
in the distribution of As. Such excitations are called spinous. They have dressed energy — ei(A) and dressed 
momentum pi(A). 

• There is an infinite number of different types of gapped excitations that carry charge but no spin, they correspond 
to adding a fc-A string of length n to the ground state. Their dressed energies are e^(A), their dressed momenta 
are p^(A). 

Let us emphasize that this is a classification of elementary excitations in the repulsive Hubbard model below half- 
filling. It is important to distinguish these from "physical" excitations, which are the permitted combinations of 
elementary excitations. In other words, not any combination of particle-hole excitations, spinons and fc-A strings is 
allowed, but only those consistent with the selection rules (p8|). To illustrate how this works and relate our findings 
to known results in the literature we consider several examples. We introduce the following terminology: we call the 
set {N, M n , M' n \n = 1 . . . oo} of the numbers of real fc's, A-strings of length n and fc-A-strings of length n occupation 
numbers of the corresponding excitation. This is in contrast to our usage of the term quantum numbers, which is 
reserved for the eigenvalues of energy, momentum, S z , S ■ S,i] z and if ■ if. 



Example 1 : Particle-hole excitation. 

This is a two-parametric gapless physical excitation with spin and charge zero, i.e. its quantum numbers as well as 
its occupation numbers are the same as the ones of the ground state. It is obtained by removing a spectral parameter 
fc/i with \kh\ < Q from the ground-state distribution of fc's and adding a spectral parameter k p with \k p \ > Q. Its 
energy and momentum are 

E ph = n{k p ) - n(k h ) , 

P P h = Pfe) - p(fc/J • (175) 

This excitation is allowed by the selection rules ( p8|) as in the ground state only the (half-odd) integers \Ij\ < (N — 1)/2 
are occupied and thus the possibility of removing a root corresponding to < (N — l)/2 and adding a root 
corresponding to \I p \ > (N — l)/2 exists. This excitation was first studied by a different approach in [^8). In Fig. || 
we show the particle-hole spectrum for densities n — 0.6 and n = 0.8. As we approach half-filling the phase-space for 
particles shrinks to zero. 
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Example 2a : Spin triplet excitation. 

Let us consider an excitation involving the spin degrees of freedom next. If we change the number of down spins 
by one while keeping the number of electrons fixed we obtain an excitation with spin 1 . Recalling that in the ground 
state we have N electrons out of which Mi = N/2 have spin down, the excited state will have occupation numbers N 
and Mi = N/2 - 1. The selection rules @ then read 

-\<h<\, l£l<£. (176) 

The first condition is irrelevant as we are below half filling, but the second one tells us that there are two more 
vacancies than there are roots. In other words, flipping one spin leads to two holes in the distribution of As. There 
is one more subtlety we have to take care of: changing the number of down spins by one, while keeping the number 
of electrons fixed leads to a shift of all Ij in (|25| ) by either | or — h. The consequence of this shift is a constant 



JV 
L 

The physical excitation obtained in this way is a gapless two-spinon scattering state with energy and momentum 



contribution of ±71" -j- to the momentum of the excited state. This leads to two "branches" of the same excitation 



Strip = -ei(Ai) - ei(A 2 ) , 

TV 

Ptrip=Pl(Ai)+ P i(A 2 )±7T- . (177) 



In Fig. ^| we show the spin-triplet spectrum for densities n = 0.6 and n = 0.8. 
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In the Hubbard model the spin-triplet excitations were first studied by Ovchinnikov |27]] and by Coll p8[ . The 
situation encountered here is similar to the spin- 1/2 Heisenberg chain |2(J in the sense that the spin-triplet excitation 
is a scattering continuum of two spin-1/2 objects. Furthermore there is a spin-singlet excitation, which is precisely 
degenerate with the triplet (see Example 2b below) . This fits nicely into a picture based on spin-1/2 objects: scattering 
states of two spinons give precisely one spin 1 and one spin multiplet | <E> \ — 1 © 0. Finally, when we approach 
half-filling, the spin-triplet continuum constructed above goes over into the 5 = 1 two-spinon scattering continuum 
of the half-filled Hubbard model |l!|§0|. 

On the other hand there are differences as well: in the less than half-filled Hubbard model the Fermi momentum is 
generally incommensurate, which leads to incommensurabilities in the spin excitations (see Fig. |^). 

More importantly, it is always possible to combine any type of excitation with a particle-hole excitation. It is 
therefore not possible to distinguish the two-parametric spin-triplet excitation constructed above from the special 
case of a four-parametric excitation, where a particle-hole excitation sits "on top" of the spin-triplet excitation and 
where the momenta of the particle and the hole are fixed at the Fermi rapidity. In other words, due to the presence 
of gapless particle-hole excitations there is an inherent ambiguity in the interpretation of the excitation spectrum on 
the basis of an 0(1) calculation of energy eigenvalues of the Hamiltonian. 

Example 2b : Spin singlet excitation. 

Let us now choose the occupation numbers as TV, Mi = y — 2, M% = 1. The corresponding state has the same 
quantum numbers as the ground state. From (|2^) we find that there are N/2 vacancies for real A's and thus two 
holes corresponding to rapidities Ai and A2. In other words the excitation considered involves two spinons. As far as 
the 2-string is concerned we find that the associated integer must be zero J\ — 0. The same shift as in example 2a 



occurs for the Ij's. Using (165) and ( 16S ) we obtain the energy and momentum of the associated excitation 



E sing = -ei(Ai) -ei(A 2 ) , 

N 

P sing =pi(Ai)+pi(A 2 )±7r— . (178) 

We see that the spin singlet is precisely degenerate with the spin triplet considered above. This is a consequence of 
the spin SU (2) symmetry of the Hamiltonian in zero magnetic field. 



Example 3 : k-A string of length 2. 
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Let us consider the simplest excitation involving a fc-A-string. One possibility is to choose the occupation numbers 
as N, Mi = N/2 — 1, M[ = 1. In addition we keep the distribution of Ij fixed in such a way that Ij+i — Ij = 1. It is 
easily checked that this excitation is allowed by (ps|). Its energy and momentum are 

E k _ A = e[(A) , P fe _ A =pi(A), (179) 

where A G (—00,00). 




FIG. 7. Dispersion of a fe-A string excitation of length 2 for several values of U and density n. 

In Fig. [7] the dispersion of a fe-A string of length 2 is shown for U — 0.5 and U = 2.0 and several values of density n. 
We see that the range of momenta collapses to zero as we approach half-filling. At the same time the dressed energy 
approaches zero. This is in agreement with the results for a half-filled band [Q, where both the dressed energy and 
the range of momentum are identically zero. 

In order to further exhibit this collapse we subtract the offset — 2/i — U. The resulting curves are displayed in Fig. @. 
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FIG. 8. Dispersion of a k-A string excitation of length 2 for several values of U and density n, where the contribution 
-2fi — U has been subtracted. 
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F. Excitations in the half-filled band 



In the limit of a half-filled band /i — > the Fermi-rapidity Q tends to n and the excitation spectrum simplifies 
drastically (30|. We find that e^(A) = Vn and the only non- vanishing dressed energies are 



ei(A) 



«(fc) = -2 cosfc -U/2-2 [ 

Jo 



duj J\{lo) cos(u;A) 
u> cosh(w[//4) ' 

dco Ji(w)cos(w sinfc)e - "' 7 / 4 



a; cosh(w?7/4) 

(180) 



where Jo,i are Bessel functions. The charge excitations are now gapped as a result of the Mott-Hubbard transition. 
The complete spectrum of physical excitations is derived in detail in [j30| (see in particular the appendix). It is given 
in terms of spinon and holon scattering states forming representations of SO (4). 



G. Relation to the root-density formalism 



There exists a second standard approach to the zero temperature excitation spectrum in Bethe Ansatz solvable 
models, which in the case of the Hubbard model has been employed in for example in p7j , p8| . We call this the 
root-density formalism (RDF). In this approach one specifies a distribution of (half-odd) integers Ij, J™, J'™ in the 
discrete Takahashi equations (p5|)- (27 ) . One then takes the thermodynamic limit using that the distribution of roots 
becomes dense, so that (p5[)-(p7f) turn into a set of coupled linear integral equations for root densities. From these 
one calculates energy and momentum via the thermodynamic limit of ([52|) and ( p3[ ) . We will now show how to relate 
the results of the Yang- Yang approach we implemented above to the RDF. For dcfinitcness we consider the example 
of a k-A string excitation of length n. 

We start by rewriting the integral equation for n(k) in the following way 



(1 - K) * K 



= / dk' [S(k-k')- cos k' i?(sin k' - sin jfc)] n(k') = -2 cos k - fx - U/2 
k J-Q 



(181) 



Here the kernel of K is given by K(x, y) — cosyR(smy — sin a;). Equation ( |l8l| ) is solved by the Neumann series 



00 i-Q 

/ rf fc'(2 cos k' + n + U/2)K l (k, k') , 
i=o J -Q 



(182) 



where K l (x,y) denotes the /-fold convolution of the kernel K(x,y). Using ( 182 ) in ( |l64| ) we obtain the following 
representation for the dressed energy of a fc-A string of length n 



e' n (A) = ARe^/l - (A - inU/4) 2 - 2nn - nU 
cosfc nU/ 4: 



-E 



= 



dk 



dk' 



Q 



7T (nil/ A) 2 + (sin A; - A) : 



7 K l (k,k') (2 cos k' + n + U/2) 



Let us now define a function p^(k\k) by 



dk' [5{k - jfc') - K(k',k)}p^ s {k'\A) = 



cos k 



nU/4 



it (nU /4) 2 + (sin k-A) 2 ' 
The solution of (184) is given by the Neumann series 

p* s (fc|A) = E [ Q —K\k>,k) cosfc' l U/4 , 

J-q 7T (nU/4) 2 + (A - sink') 2 



i=o ■ 



Using (|185| ) in (183) we finally obtain 

e' n (A) = 4Rev/l - (A - inU/4) 2 - 2np - nU 



dk(2 cos k + fi + U/2) p b n s (fc| A) 



(183) 



(184) 



(185) 



(186) 
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H. Summary 

In this section we have reviewed the Yang- Yang approach pQ] to the thermodynamics of the Hubbard model Hi 



We have shown how to express the Gibbs free energy per site (157) in terms of the solutions of the infinite set of 



coupled nonlinear integral equations (154)-( |l5q ). 

By taking the zero-temperature limit of the thermodynamic equations we have obtained the complete classification 
of the spectrum of elementary excitations of the Hubbard model below half-filling for vanishing magnetic field. 

We emphasize that this approach is based on Takahashi's string hypothesis jL8| (see section II. C). It does not only 
provide the dispersion curves of all elementary excitations in the zero temperature limit, but also a set of 'selection 
rules' (p5|)-(p7|). These rules determine the set of physical excitations, i.e. the allowed combinations of elementary 
excitations (see section VI. E and, for the half-filled case, |2jjj3(|). 

The numerical calculation of the Gibbs free energy from ( [157D o r fll5S| ) is difficult, since it involves the solution of 
an infinite number of coupled non-linear integral equations ( |154[ )-(L56). A truncation scheme has to be introduced 
[B]]|62] , which restricts the numerical accuracy of the calculation. 



In the next section we present a different approach to the thermodynamics of the Hubbard model, which circumvents 
such difficulties. 



VII. THERMODYNAMICS IN THE QUANTUM TRANSFER APPROACH 

In this section we present a treatment of the thermodynamic properties of the Hubbard chain in a lattice path 
integral formulation with a subsequent eigenvalue analysis of the matrix describing transfer along the chain direction 
(quantum transfer matrix, QTM). This approach has several advantages. First, the analysis is very simple as only the 
largest eigenvalue of the QTM is necessary in order to calculate the free energy. This has to be compared with the 
traditional TBA (see Sect. VI) where all eigenvalues of the Hamiltonian have to be taken into account. Second, 
the number of 'density functions' and integral equations obtained below is finite in contrast to [fyl (see Sect. VI A). 
Finally, the finite temperature correlation lengths can be derived through a calculation of next-largest eigenvalues of 
the QTM |^,[7(j. This however, will not be explored in this report. 

Within the QTM approach we obtain several results of which we want to point out the conceptual achievements. 
First, the data obtained for various physical quantities agree with those obtained in |6l]|32| based on Takahashi's 
string hypothesis. Judging from this we do not see any evidence for a failure of Takahashis's formulation based on 
strings. Second, in the low temperature limit our integral equations quite naturally yield the Tomonaga-Luttinger 
liquid picture of separate spin and charge contributions to the free energy as suggested in 0,^2) ■ Mathematically, 
the dressed energy formalism known from the ground state analysis is recovered. 



A. The classical counterpart 

There are many direct path integral formulations of the Hubbard model, see for instance f^Q. For our purposes we 
want to keep the integrability structure as far as possible. To this end it proved useful to employ the exactly solvable 
classical model corresponding to the Hubbard chain, namely Shastry's model which is a two-sublattice six-vertex 
model with decoration [f42| . 

More precisely, the vertex weights of the classical model for the case U — are given by the product of the vertex 
weights of two six-vertex models (with components a and r): l\^{u) = ^ 12 (w) ®^i2( u )> where 

t\A u ) = |( cos (") + sin(u)) + |(cos(u) - sm(ii))0i<7| + {a\a^ + <rf er^) . (187) 

Taking account of the U interactions, the following local vertex weight operator (denoted by S) was found |42] 

Si,2(v, u) = cos(u + v) cosh(h(v, U) - h(u, U)) li$(v — u) + cos(v - u) smk{h(v, U) - h(u, U)) l\,<z{u + v) <rf rf , (188) 

where sinh(2fe(u, U)) :— -jsin(2u). The Yang-Baxter equation for triple S matrices was conjectured p2| , but only 
recently proved in |40] . 

The so-called L-operator is related to S by 

Li, g {u) =S i<g {u,0) , (189) 
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where i and g are indices referring to the ith lattice site and the auxiliary space, respectively. For the L-operator 



the proof of the Yang-Baxter equation with S as intertwiner was already given in |42|. The commutativity of the 
row-to-row transfer matrix 

T{u) :=trjji ii9 (u) (190) 

i 

is a direct consequence. 

Next we define Ri^(u,v) = Tli^Si^iv, u)\u—>-u, where 11^2 is the permutation matrix. The matrix elements 
R^(u,v) will be considered as the local Boltzmann weights associated with vertex configurations a, (3, /i, v on the 
lower, upper, left, and right bond, where the spectra^ parameters u and v "live" on the vertical and horizontal bonds, 
respectively. For later use we introduce R(u, v) and R(u, v) (u and v associated with the vertical and horizontal bond) 
by clockwise and anticlockwise 90° rotations of R, or in matrix notation 

R^{u,v) = R^{v,u) , R^(u,v) = R^(v,u) . (191) 



Similar to (189) and (|190| ) we can associate a row-to-row transfer matrix with R. We note the Hamiltonian limits 
T(u) — exp(iP + uH + 0{u 2 )) and T{u) — exp(— iP + uH + 0(u 2 )). Consequentially, the partition function of the 
Hubbard chain at finite temperature T — 1/(3 is given by 

Z = lim tre~ /3ff = lim lim tr [T{u)T(u)] N/2 \ u=f}/N . (192) 

We regard the resulting system as a fictitious two-dimensional model on a L x TV square lattice. Here N is the 
extension in the fictitious (imaginary time) direction, sometimes referred to as the Trotter number. The lattice 
consists of alternating rows each being a product of only R weights or of only R weights, respectively. Now by looking 
at the system in a 90° rotated frame which turns R and R weights into R and R weights, it is natural to define the 
'quantum transfer matrix' (QTM) by 

N/2 

{T QTU ( U ,v)}\§ :=Enc;;Lu^:r(^) . ( i93 ) 

which is identical to the column-to-column transfer matrix of the square lattice for v = 0. The interchangeability of 



the two limits (L, N — > oo) |63 6J] leads to the following expression for the partition function, 

L 

(194) 



Z = lim lim tr 

N — >oo L — >oo 



Tqtm ( u = — ,0 



There is a gap between the largest and the second largest eigenvalues of 7qtm(u,0) for finite (3. Therefore the free 
energy per site is expressed just by the largest eigenvalue A max (u,0) of 7qtm(u,0), 

f = -T lim In A max (u = ^-,o] . (195) 



iiioa i — ■» r 

N^oo V N 



It is relatively simple to see that (193) is integrable, i.e. a family of commuting operators for variable v and fixed 
A non-vanishing chemical potential /i and magnetic field B can be incorporated^ as they merely lead to trivial 



u 

modifications due to twisted boundary conditions for the QTM (cf. J70|) 



B. Diagonalization of the Quantum Transfer Matrix 



Here we summarize the main results of |34{ where the diagonalization of (|193| ) on the basis of an algebraic Bethe 
ansatz was performed. Note that the general expression for the eigenvalue of the quantum transfer matrix is quite 
complicated |34], but simplifies considerably at v — and u — » 0, 



4 We note that our conventions for the magnetic field in this article are different from J34|. In |m] the magnetic field was 
denoted by H = 2B. 
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A(v = 0) = eP u / 4 (l + e^ +B) )(l + e^- B ))u N Y[ Zj . (196) 

3=1 

The numbers Zj are charge rapidities satisfying Bethe ansatz equations which are most transparently written in terms 
of the related quantities 

Sj = 1(^-1) . (197) 

For these rapidities Sj and additional rapidities w a the coupled eigenvalue equations read 

e -P{M-B) 6( \ _ I2(sj ~ \U/A) 20u q 2 (w a +iU/2) _ qi (w a +iU/A) 

9[J> q 2 (s j +iU/4) ' q 2 (w a -iU/2) qi (w a - iU/4) ' 



(198) 



where we have employed the abbreviations for products over rapidities 

<7i 00 =II( S -^) > l2(s)^H(s-w a ) . (199) 



The function <f> is defined by 



(1 -*_/*(*))(!- 24-M*)) • 



\ N/2 

\(l + z-/z(s))(l + z + /z(s))j V 7 

where z(s) possesses two branches. The standard ("first") branch is chosen by the requirement z(s) ~ 2is for 
large values of s, and the branch cut line [—1,1]. The numbers z± are defined by z± = exp(a)(tanM) ±:L , where 
sinh(a) = — -jsin2w. 



As usual, equations like (198) which are identical in structure to (Q) and (0) have many solutions. In our approach 



to the thermodynamics just the largest eigenvalue of the QTM matters. The corresponding distribution of rapidities 
is relatively simple. For fj, = B = Q the rapidities are all situated on the real axis. Naturally, for finite fi and B we 
have modifications which, however, do not affect qualitative aspects of the distribution. 

In a similar way a path integral formulation of the Hubbard chain and a diagonalization of the corresponding QTM 



was employed in 75 33]. In ]75l the eigenvalue equations were studied numerically for finite Trotter number and the 



case of half-filling. In [p3| an analytic attempt was undertaken to study the limit of infinite Trotter number and to 
derive a set of non- linear integral equations. Unfortunately, these equations were rather ill posed with respect to 
numerical evaluations. In the next section we present a formulation of non-linear integral equations on the basis of 
recent work fill. 



C. Non-linear integral equations 



In this section we are concerned with the derivation of well posed integral equations equivalent to (19S) for the 
largest eigenvalue of the QTM and thus for the fre e en ergy per site (see (195)). We introduce a set of auxiliary 
functions (b, c, and c) described in more detail in ( 207 ) below. These auxiliary functions are complex functions, 
however mostly evaluated close to the real axis. 

In terms of the auxiliary functions the Gibbs free energy per site is expressed in several ways 



/ = -M - j ~ ^:J[hyz(s)}' In (1 + c + c) ds 



+ 



T 
47ri 



In 



z{s - iU/2) 



z(s) 



U 



T 



= T-^ZT Ms)]' In 



ln(l + b(s))ds 



T 



T 
47ri 



hi 



z{s + \U/2) 



z(s) 



4 2?ri 
For yet another expression see 



:ds ~7T^ / [lnz(s-if//2)]'ln(l + c(s))ds 

27T1 



ln(l + l/b(s))ds , (201) 

(202) 
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Equations (201) and fl202| ) have to be compared with equations (157) and (159) which give the free energy in the 
string based approach of Takahashi. In contrast to the string based approach, the auxiliary functions b, c and c 
entering (201) and (202) satisfy a closed set of finitely many (non-linear) integral equations, 



In b = -2/3B + K 2 D ln(l + b) - ~K[ o ln(l + 1/c) , 
In c = -0U/2 + i3(n + B) +(p — K[u ln(l + 1/b) - ~K{ o ln(l + c) 
In c = -/3Z7/2 - (3(fi + B) - <p + K x O ln(l + b) + K x o ln(l + c) . 
Here we have used the definition 



<p(x) = -2/3ixy/l - l/x 2 



and have introduced the notation 



(9°f)(s) 



g{s - t)f{t)dt 



(203) 



(204) 



(205) 



for the convolution of two functions g and / with contour C surrounding the real axis at infinitesimal distance above 
and below in anticlockwise manner. The definition of □ is similar, with integration contour surrounding the real axis 
at imaginary parts ±{7/4. 

The kernel functions are rational functions, 



Ki(s) 



U/Air 



s(s + iU/2) 



U/Att 



s(s - iU/2) 



K 2 (s) 



U/2n 



U 2 /A 



(206) 



Next, we want to point out that the function b will be evaluated on the lines Im s = ±17/ A. The functions c and 
c need only be evaluated on the real axis infinitesimally above and below the interval [—1,1]. Also the convolutions 
involving the "c functions" in (203) can be restricted to a contour surrounding [—1,1] as these functions are analytic 
outside. 



Lastly, we want to comment on the derivation of ( p02| ), (205). The explicit expressions of the functions b, c, c are 

r h+h+h+h 



h+h+ h + k 
h + h 



h±l2±h± h 



h + h 
__ h + h 



h+h+h+h+h+h+h+ h 
h+h+h+ h 



(207) 



where the functions lj and lj are 



h +h h+h+h+h+h+h+ h + h 



h (s) = Xj (s - iU/A) ■ e 2 ^cj) + (s)<j)- (s) , lj (s) = Xj (s + iU/A) , 
and the Xj are defined in terms of the q± and q 2 functions, i.e. in terms of the Bethe ansatz rapidities 



(208) 



X 1 (s) 
A 3 (s) 



c 0(u+B) til - 1[/ / 4 ) 

qi(s-iU/4) 

q 2 (s + iU/2) 
q 2 (s) qi ( S + iU/A) ' 



X 2 (s)=e^ ^ S - {U ' 2 ) 
M[S) q 2 (s) qi ( S -iU/A) 



A 4 (s) 



0(s + i?7/4) qi {s + iU/A) 



(209) 



The functions defined in (207) are proven to satisfy a set of closed functional equations which can be transformed 
into integral form ( ^03[ ), cf. [34|. Also ( ^02| ) follows from ( [207| ) after a lengthy yet direct calculation. The merit of 
(202), (203) is that this formulation does no longer make any reference to the Bethe ansatz equations! Hence the 
calculation of an infinite set of discrete rapidities is replaced by the computation of analytic functions for which much 
more powerful tools are available. 
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D. Analytical solutions of the integral equations 



Before numerically studying various thermodynamic properties for general temperatures and particle concentrations 
we want to give some analytic treatments of limiting cases of the Hubbard model. 



1. Strong- coupling limit 

In the strong-coupling limit U — > oo at half-filling (// = 0) the Hubbard model is expected to reduce to the Heisenberg 
chain. Indeed, in the strong-coupling limit we find that c, c — » 0. Hence, the only non-trivial function determ ining 
the eigenvalue of the QTM is b, see the first expression of (202). The integral equation for b as obtained from (203) 
is identical to that obtained directly for the thermodynamics of the Heisenberg model |34[] . 



2. Free-Fermion limit 



Next, let us consider the limit U — > leading to a free fermion model, however representing a non-trivial consistency 
check of the equations. Indeed, the auxiliary functions can be calculated explicitly. Finally, the free energy per site 
reads 

f = -^~ I ln{[l + exp((/i + B + 2cosfc)/r)][l + exp((/i-S + 2cosfc)/r)]}dfc , (210) 

2 71 " J-7T ' 

which, as desired, is the result for free tight binding electrons. 



3. Low-temperature asymptotics 



The low-temperature regime is the most interesting limit as the system shows Tomonaga-Luttinger liquid behavior. 



We want to describe the relation of the non-linear integral equations to the known dressed energy formalism 71 7^] 
of the Hubbard model. This represents a further and in fact the most interesting consistency check. 

For T = X/P — y we can simplify the non- linear integral equations as they turn into linear integral equations, 
however in the generic case {B and /i ^ 0) with finite integration contours. This can be seen as follows. To be specific 
let us adopt fields B > 0, fi < (particle density n < 1). In the low-temperature limit several auxiliary functions 
tend to zero, namely b(s) — » on the line Im s = —U/4, c(s) — > above and below the real axis, and l/c(s) — > just 
below the real axis, i.e. on the line Im s = — e. The remaining non-trivial functions are b on the line Im s = +U/4 
and 1/c just above the real axis for which we introduce the notation 



6(A) = b(A + W/4) , c(A) = l/c(A + ie) 



(211) 



We note that 



|c| » 1 for \x\ < Aq, o-q and \b\, |c| <C 1 for |sc| > Aq, er 



(212) 



for certain crossover values Ao, <Jq. The slopes for the crossover are steep, so that the following approximations to the 
integral equations (203) are valid, 



lnb = 



° fc 2 (A - A') lnfe(A') dA' + f ° ki (A - sin k') cos k' In c(fc') dk' , 

A J— k 



lnc = 4> c + / k\ (sin k — A') In b(X') dA' 



-An 



(213) 



(214) 



where fci(A) = K±(\ — it// 4) = K\(X + if// 4) and fc 2 (A) = i4T 2 (A). In order to facilitate comparison with the dressed 
energy formalism we also introduced a new i nteg ration var iable k leading to a change of the boundaries of integration, 
co — > fco = arcsinero- The driving terms in ( [213| ) and (214) are related to the bare energies 



B 



-2cosfc - /j, — U/2 — B 



(215) 
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by 

cp b = -f3e s + O(l/f3) and C = -/3 £ ° + 0(l/j3) . (216) 
Therefore, we find the following connections between auxiliary functions and the dressed energy functions, 

ln&= -/3s s + 0(l//3) , lnc= -0e c + O(l//3) . (217) 



For a comparison with |7l|]72| ] note the different normalization of the chemical potential. 

For a comparison with the results in section VI we set the magnetic field equal to zero. Then e Q s = B = 0, and it 
can be seen that Ao = oo. Setting Q — ko we can identify (213), (214) in the zero temperature limit with (165). Wc 
find 

- lim Tin 6(A) = lim ei(A) = lim Tlnm(X) , (218) 

- lim Tlnc(A) = lim n(k) = lim TlnC(fc) . (219) 

The free energy also admits the above approximation scheme, yielding up to C(T 2 )-terms the low-temperature 
expansion 

/ = eo-~ (- + -V 2 . (220) 

Here the definitions of the sound velocities and the ground state energy are standard. The additive occurrence of l/v c 
and 1 jv a on the right hand side of ( |22C| ) is a manifestation of spin-charge separation in the one-dimensional Hubbard 
model, due to which each elementary excitation contributes independently to ( |220| ) . The velocities v c and v s typically 
take different values. 



4- High-temperature limit 

Finally, we consider the high-temperature limit T — > oo with B, U as well as 0fi fixed ensuring a fixed particle density 
n. The integral equations turn into algebraic equations which are easily solved, resulting in the high-temperature 
limit 

5 = 21a(^)-nln(^) (221) 

for the entropy, as expected by counting the degrees of freedom per lattice site. Especially at half- filling, n — 1, this 
equals to S = ln(4). 



E. Numerical Results 



Here we show numerical results for specific heat C, magnetic susceptibility Xm, and charge susceptibility \ c for 
half-filling and small doping, see Figs. ^|, JT^, O. In addition, we aim at a comparison of results |^l|,^2| obtained within 
the thermodynamic Bethe ansatz |l6| ] based on the string hypothesis, and results |34| ] obtained within the quantum 
transfer matrix approach. 

In essence, we observe convincing agreement of the data obtained within the two absolutely different approaches. 
We conclude that there is no indication for any failure of Takahashi's formulation of thermodynamics based on the 
string hypothesis. Of course rather small differences exist in the sets of data. A more thorough comparison shows 
deviations of the data presented in Fig. ^| and [l^. This is nothing but expected due to the truncation procedure 
adopted in ]6l]j62] ]. Instead of dealing with an infinite set of integral equations for density functions (151,154-156) 
for string excitations of spin rapidities (|2^) and charge rapidities (p3|), a finite subset was taken into account. In the 
case of complex charge rapidities, strings describe excitations with gap. These degrees of freedom are less sensitive to 
errors introduced by the truncation than strings involving only spin rapidities which describe gapless excitations. In 
fact, the agreement of the data for the charge susceptibility Xc is best, small deviations are observed in C and Xm- 



4G 



The advantage of the QTM approach over the traditional TBA is threefold. First, in |Q the Hubbard chain with 
very strong doping was analyzed (apparently not possible in the traditional TBA) and quite unexpected structures in 
the susceptibility data were found. In the magnetic susceptibility only traces of spinon excitations were visible. The 
charge susceptibility, however, exposed holon signatures and additional maxima due to spinon excitations indicating 
deviations from the concept of spin-charge separation. Second, the accuracy of numerical data obtained within the 
QTM approach is much higher as it is more efficient to deal with a set of integral equations which is strictly finite 
from beginning. Finally, even correlation lengths can be calculated within the QTM approach although not presented 
in detail yet. 




FIG. 9. Specific heat C versus temperature for U = 8 and particle densities n = 1, n = 0.9, n = 0.8 and n = 0.7. The left 
graph shows results obtained within the string based approach, the right graph shows results obtained by the quantum transfer 
matrix method. 
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FIG. 10. Magnetic susceptibility \ m versus temperature for U — 8 and particle densities n = 1, n = 0.9, . . . , n — 0.5. The 
left graph shows results obtained within the string based approach, the right graph shows results obtained by the quantum 
transfer matrix method. 
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FIG. 11. Charge susceptibility Xc versus temperature for U — 8 and particle densities n = 1, n = 0.9, n = 0.8 and n = 0.7. 
The left graph shows results obtained within the string based approach, the right graph shows results obtained by the quantum 
transfer matrix method. 



F. Equivalence of string based thermodynamics with QTM approach 

Here we want to address the fundamental and perhaps puzzling question how two exact, however completely different 
formulations of the thermodynamics of the Hubbard model as presented in sections VI and VII can exist. It is very 
important to understand this problem as the means of derivation and the mathematical properties of both approaches 
are seemingly lacking any similarities. 

In the remainder of this section we restrict ourselves to a sketch of the general mathematical structures underlying 
the relation of the traditional thermodynamical formulation based on typically infinitely many density functions and 
the new formulation with strictly finitely many auxiliary functions. This relation has been worked out for a number 
of models including spin chains [[39| and correlated fermion models like t-J systems (7{J, not yet however for the 
Hubbard chain. In the following we want to introduce the general technique for spin-1/2 Heisenberg models and 
related systems. 

The starting point is an integrablc Hamiltonian corresponding to an exactly solvable classical model as discussed 
in section VILA. For this model a path integral representation is formulated leading to an integrable QTM T(v) with 
eigenvalue equation 

rpi , f(v-X)q(v + X) + f(v)q(v-X) _ 9 . 

T{v) = W) ' (222) 

where f(v) is a known function and q(v) is a polynomial or a product over trigonometric functions with zeros exactly 
at the Bethe ansatz rapidities. By a line of reasoning similar to that presented in VII. C one non-linear integral 
equation can be derived which yields full information on the free energy of the model. 

The key to understanding the relation of the standard TBA to the above approach is the notion of fusion algebras 
and inversion identities. By repeated use of the fusion procedure a set of transfer matrices T(v) = T 1 ^), T 2 (v), 
T 3 (v),..., is generated satisfying 

T»T> + A) = /(«- X)f(v + X)I + T q+1 (v)T q -\v + X) , (223) 

where q takes integer values, q = 1, 2, . . ., and T («) is proportional to the identity operator. Introducing operators 

{V) f(v-X)f(v + qX) {22A) 
we find the extremely useful inversion identity hierarchy 

Y q (v)Y q (v + X) = [I + Y q - 1 (v + X)][I + Y q+1 (v)} . (225) 
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This set of functional equations can be transformed into a set of non-linear integral equations for the largest eigenvalue 
of T(v). These integral equations appear to be identical to the equations obtained along the traditional TBA approach! 



Several remarks are in order. Relations (225) were derived in [ |77|j7q ] starting with the TBA equations, however 
the connection to the fusion hierarchies was not made. In contrast to the TBA integral equations the functional 
equations ( |225| ) (as of course equation ( |222| )) admit more than one solution. The physical meaning of these solutions 
becomes clear only through the microscopic construction of transfer matrices. The next-leading eigenvalues describe 
the correlation lengths of static correlation functions at finite temperature. Such applications are investigated in 
current research. Some results have been reported in |7(i||79| . 

In summary, the traditional TBA equations are not in contradiction with the new thermodynamics based on the 
QTM. Rather on the contrary, the QTM approach represents a unified approach to both sets of integral equations. 
The central object is the study of the largest eigenvalue of the QTM. This can be achieved by two different methods, 



either by a Bethe ansatz (222) or by use of the fusion hierarchy (223). 

In physical terms, the appearance of the fusion hierarchy is related to the pole structure of the intertwiner, i.e. 
the i?-matrix satisfying the YBE. In turn the i?-matrix is related to the S'-matrix of particles, the poles are related 
to bound states (strings) which in turn are the central objects of the traditional analysis of the thermodynamics of 
integrable systems. In this way we observe mathematical and physical connections between the initially so differently 
looking approaches. 

VIII. CONCLUSIONS 

In this article we considered fundamental questions concerning the exact solution of the one-dimensional Hubbard 
model. In appendices A and B we presented a detailed derivation of the wave function. We studied solutions of the 
Lieb-Wu equations. We concentrated our attention on fc-A strings. Unlike A strings they were never carefully studied 
in the literature before. We arrived to the conclusion that, for fixed coupling and large lattice size, k-A strings deviate 
from their ideal positions in a controllable way. In the infinite chain limit k-A strings approach their ideal positions. 
We also reviewed the thermodynamics of the model. There are two different approaches to thermodynamics: one is 
based on strings, whereas the other one is not. Both approaches show convincing agreement for the calculation of 
the bulk thermodynamic properties of the Hubbard model. Passing to the zero temperature limit we obtained the 
dispersion curves of all elementary excitations at zero magnetic field, below half-filling. 

We would like to conclude with pointing out some interesting open problems: 

(i) Perhaps the most fundamental open problem is the calculation of the norm of the wave function (|J)-(^]). Formulae 
for norms of wave functions are known for Bethe ansatz solvable models with only one level of Bethe ansatz 



equations, like the Bose gas with delta interaction |36 82 or the XXX and XXZ spin-i chains [p2j (see also 



[p5[). For these models the norms are expressed as determinants of the Jacobians of the Bethe ansatz equations. 
The calculation of the norm of the Bethe ansatz wave functions may be considered as a first step towards the 
calculation of determinant representations of correlation functions f35| . 

(ii) Another interesting open problem which we already mentioned above is the calculation of the finite temperature 
correlation length within the quantum transfer matrix approach. This requires to calculate the second largest 
eigenvalue of the quantum transfer matrix. One has to overcome certain technical subtleties coming from the 
fact that the Hubbard model is a model of Fermions. A formalism capable of calculating the correlation length 
of integrable fermionic systems has recently been developed [p3],p0|]8l|l . 



(iii) At zero temperature some of the correlation functions of the Hubbard model show a power law decay. Confor- 
mal field theory^ naturally describes these powers (the set of conformal dimensions) j7lj . Still, there are open 
problems within the conformal approach. For example, since the expansion of the lattice operators in terms of 
conformal fields is not known explicitly, the resulting expressions contain unknown amplitudes. Some of these 
amplitudes may actually vanish. Recently |85| the vanishing of the amplitudes corresponding to density corre- 
lations for the half- filled model was shown by use of the SO (4) symmetry. For a more thorough understanding 
of correlation functions an identification of the operators which are of interest for the Hubbard model with the 



In the context of condensed matter physics conformal field theory is equivalent to Luttinger liquid theory 
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standard operators of conformal field theory would be important. Recent work on a scaling limit of the Hubbard 



model [ 86 - 88 may prove to be useful in this context. 



(iv) The predictions of the conformal approach to correlation functions are limited to large distances, corresponding 
to very low energies. Theoretically as well as from the point of view of recent experiments on quasi one- 



dimensional structures in solids (e.g. [|31 32 1) it would be highly appreciable to have a method for the calculation 
of correlation functions at all energy scales. This problem was recently tackled in |8^] within the form factor 
approach ]90|-|92]|, which was originally designed for integrable 1+1 dimensional quantum field theories. In ]89[ | it 
was argued that the form factor approach might as well apply to the half-filled Hubbard model, and a formula for 
the two-spinon form factor of the spin-operator Sj was presented. It would be interesting to extend the result 
to form factors of electronic creation and annihilation operators, as such kind of extension could be directly 
applied to the interpretation of the angle-resolved photo emission spectroscopy data of |3l]|32| . 

(v) Despite the progress in the understanding of the mathematical structure of the Hubbard model, which was 
achieved over the past few years and which we briefly discussed in the introduction, we still feel uncomfortable 
with the present stage of our knowledge. Shastry's i?-matrix |f4~i|-[43|| , which is the key for our present under- 
standing of the algebraic structure behind the Hubbard model, is unusual as compared to i?-matrices of other 
integrable models. It does not possess the so-called difference property, i.e. it is not a function of the difference 
of the spectral parameters alone. The S-matrix at half-filling [29;] , on the other hand, possesses the difference 
property and can therefore be associated with a Y(su(2))©Y(su(2)) Yangian jlT]]. The precise relation between 
R- matrix and 5*-matrix is only understood in the rather simple situation of an empty band f4S| , |4"9| . Because of 
the lack of the difference property we can neither find a boost operator for the Hubbard model by the reasoning 
of f93f| nor can we associate a spectral curve with it. 

Another problem is the dimension of the elementary L-operator related to Shastry's i?-matrix, which is 4 x 4 
(rather than 3 x 3 as one could guess naively from the fact that the Bethe ansatz for the Hubbard model has two 
levels). For this reason there are too many candidates for creation and annihilation operators in the algebraic 
Bcthe ansatz p2|^4]j4q] . Again this redundancy has only been partially understood in the empty band case Q . 
The known algebraic Bethe ansatz j|4|,[l5) is of involved structure and hopefully will be simplified in the future. 

Acknowledgment. We would like to thank C. M. Hung for drawing figures 1-3 for us. We are indebted to H. 
Frahm, N. Kawakami, B. M. McCoy, A. Schadschneider, J. Suzuki, A. M. Tsvelik and M. Takahashi for helpful and 
stimulating discussions. This work was supported by the EPSRC (F.H.L.E), by the Deutsche Forschungsgemcinschaft 
under grant numbers Go 825/2-1 (F.G.) and Kl 645/3 (A.K.) and by the National Science Foundation under grant 
number PHY-9605226 (V.E.K). A.K. acknowledges support by the Sonderforschungsbcreich 341, Koln- Aachen- Julich. 



50 



APPENDIX A: DERIVATION OF THE WAVE FUNCTION 



1. General setting 

The expression (||)-(||) for the Bethe ansatz wave function of the Hubbard model was first presented by Woynarovich 
p3| . The purpose of this appendix is to give a detailed derivation of Woynarovich's wave function. We will make use 
of results obtained in references O] , O] and of the quantum inverse scattering method p5|,^l|| . Other derivations of 



the Bethe Ansatz wave functions for the Hubbard model can be found for example in [|94[ and in |95| . 

The outline of this appendix is as follows. In section 1 we define the wave function ip(xx, . . . ,xn), for N electrons 
and derive the first quantized Schrodinger equation for the Hubbard model. 

In sections 2 and 3 we present the explicit solution of the Schrodinger equation with periodic boundary conditions 
for the cases of 2 and 3 electrons, respectively. We treat these cases in considerable detail for pedagogical reasons. 
Finally, in section 4 we discuss the general case of N electrons. 

An important ingredient in the construction of the wave function is the exact solution of an inhomogeneous spin-1/2 
Heisenberg model. The essence of the nested Bethe Ansatz procedure employed in constructing wave-functions for 
the Hubbard Hamiltonian is the reduction of this problem to a simpler one, which involves only the spin degrees of 
freedom. The dynamics of the spin degrees of freedom are described by an inhomogeneous Heisenberg model, and its 
exact solution constitutes the "nesting" of the Bethe Ansatz procedure. We summarize the algebraic Bethe Ansatz 
solution of the inhomogeneous Heisenberg model in Appendix [b[ 

Let us recall the explicit form of the Hamiltonian 

L L 

H = -J2J2 + 4+i,<^>) + u YX n * - - 1) ■ (ai) 

J'=l ff=U .7=1 

As the total number of electrons N and the number of electrons with spin down M are good quantum numbers, we 
can use them to label eigenstates of ( |Al| ) 

\N,M) = Y,^^ Xl --- XN ^^---' a ^ c L^--- c t N ,, N \0) . (A2) 

{<Jj} {x k } 

Here J2{a} denotes summation over all Nl/((N — M)\M\) possible spin-configurations with M down-spins. Due to 
the anticommutation relations between the Fermion operators, we may assume without loss of generality that the 
amplitudes i\) are totally antisymmetric 

ip{x Pl ,. . . , xp N ; op x , . . . ,crp N ) = sign(P)V>(a;i, . . . , x N ; o\, . . -,cr N ) , (A3) 

where P = (Pi, P2, . . . , P/v) is a permutation of the labels {1,2,..., N}, i.e. an element of the symmetric group Sn- 
The antisymmetry property (A3) implies that the summation over spin configurations in (A2) is redundant. Indeed 
one finds that 

1 N ' M) = (N -M)!M! ^ ^ [X1 > rx:a: ^)4 1 , ffl ' ' ' |0> , (A4) 

{ x j } 

where (<7i, . . . , <jn) S Sn is arbitrary. In order to derive the Schrodinger equations it is therefore convenient to work 
with the following simplified expression for general eigenstates of H 

W,M;a) = ip(xi, . . . , x n ;<ti, ■ ■ . , crjv)4 1(T1 • ■ -a^JO) . (A5) 

{xj} 

It is now straightforward to show, that the eigenvalue problem 

H\N, M ; a) = E\N, M; a) , (A6) 
implies the following Schrodinger equation for the wave function ip [0 

N 

-^2 X! VK^i, ...,Xj+s,.. .,x N ;v) + U ^S(xj,x k )ip(xi, . . . ,x N ;cri, ■ ■ - ,<Jn) 

j=l s=±l j<k 

UN UL... 

= (E + — —)ip(x 1 , ...,x N ;cri,...,a N ) . (A7) 



Here S(a, b) denotes the Kronecker delta. 
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2. Two electrons 



Let us now explicitly construct the wave function for the case of two electrons, N — 2. The Schrodinger equation 
is of the form 



-ip(xi - l,x 2 ;ai,a 2 ) - ip(xi + l,x 2 ;<Ji,a 2 ) - ij){x\,x 2 - 1; 0-1,02) - ip(xi,x 2 + 1; 01,0-2) 



+US(xi, x 2 )ip(xi, x 2l cti, a 2 ) = (E + U 



UL 



)t/j(xi,x 2 ;ai,a 2 ) 



(A8) 



As long as x\ < x 2 or x± > x 2 (A8) reduces to the Schrodinger equation for free electrons on a lattice and its solutions 
are therefore just superpositions of plane waves. When the electrons occupy the same site, they interact. This can be 
thought of in terms of a scattering process. Due to integrability this scattering is purely elastic, which means that 
the momenta of the two electrons are individually conserved. Thus, the most that can happen is that the electrons 
exchange their momenta. These considerations lead to the famous "nested" Bethe ansatz form for the wave functions, 
which we will discuss next. 

Let Q be a permutation of the labels of coordinates i.e. Q = {Qi,Q 2 ) G {(1, 2), (2, 1)}. In the "sector" Q defined 
by the condition xq x < xq 2 the nested Bethe Ansatz for the wave function is 



ip(xi,x 2 ; oi,o 2 

Pes 2 j=i 



Substituting (A9) into (AS) for the case x\ ^ x 2 we obtain 

E = -(2cosfci + 2cosfc 2 ) - U 



UL 



(A9) 



(A10) 



When xi = x 2 we have to "match" the wave function defined in the two sectors Q = (12) and Q = (21). This 
requires single valuedness 

ip(x,x;ai,a2) = [A ai(J2 (k%, k 2 ) - A ai<72 (k 2 , h)] exp(i[fci + k 2 ]x) 

= [A tT2ai (k2,ki)-A a2ai (ki,k2)]exp(i[ki + k 2 \x) . (All) 
In addition, the Schrodinger equation ( |A8| ) for x — x\ — x 2 needs to be fulfilled, which yields the condition 
-e- tkl A aia2 (h,k 2 ) + e-' lk2 A ai<T2 (k 2 ,h) + e lk2 A a2(Tl (k u k 2 ) - e lkl A„ 2ax (fc 2 , h) 
-e lk2 A ai (j 2 (fci , k 2 ) + e ikl A ai CT2 (k 2 ,h) + e" lfel A a . 2(Tl (fci , k 2 ) - e~ lk2 A a2<Tl (h 2 ,fa) 
+U [A aia2 (ki.k 2 ) - A aia2 (k 2 , fci)] 
= -2(cosfci + cosfc 2 ) [A aia2 (ki, k 2 ) - A,j ia2 (k 2 , fci)] . 



(A12) 



By means of (A12) and (All) we can express two of the four amplitudes A aQiCrQ2 (kp 1 , fcp 2 ) in terms of the other two. 
A short calculation gives 



A aiC r 2 (k 2 , fci) 



-U/2i 



sin fci — sinfc 2 — U/2i 



A ai a 2 (ki, k 2 ) 



sin ki — sin fc 2 
sin fci — sin fc 2 — U/2i 



Aa 2 n\ (ki , k 2 ) . 



(A13) 



Equation (A13) has a natural interpretation in terms of a scattering process of two particles. In order to see this we 
rewrite it as 



A a2ai (k2,k 1 ) = S a^(kx,k 2 ) A T1T2 (k u k 2 ) , 



(A14) 



where S(ki, fc 2 ) is the two-particle S-matrix with elements 

-U/2i 



sin fci — sinfc 2 —U/2i 



n <Tin , sin fci- sin fc 2 
CT2T2 sin fci - sinfc 2 - U/2i a2T2 



(A15) 
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Here I is the identity operator JflZJJ = £, 



£71 ,Ti U (T2 ,T2 



and II is a permutation operator = S aitT2 S a 



In the next step we want to impose periodic boundary conditions on the wave function 

ip(L+ l,x 2 ;ai,a 2 ) = ip(l,x 2 ;ai,a 2 ) , 
-0(0, x 2 ; 0-1,0-2) = ■0(L,x 2 ;oi,o 2 ) , 
i/}(xi,L + l;oi,o 2 ) = V'Oei, 1; 0-1,0-2) , 
-0(a;i,O;o-i,O2) = ^(afi, £; 01, cr 2 ) • 
A short calculation shows that this imposes the following conditions on the amplitudes 

A <?Qi°Q2 ( kp i > kp 2 ) = ex P ( ik Pi L ) A <?Q 2 °Qi ( kp 2 > k Pi) 1 



(A16) 



(A17) 



where P,QeS 2 are arbitrary. In principle we now could simply solve ( A17 ) and thus determine the quantization 
conditions for the momenta fci.2- Rather than proceeding in this way, we will introduce some seemingly unnecessary 
formalism, which however will be very useful for treating the case of more than two electrons. 

We define an auxiliary spin model on a lattice with TV sites i.e. a lattice formed by the electrons. On every site j 
there are two allowed configurations | t)j and | corresponding to spin up and down respectively. Next we define 
spin operators S^' z acting on the resulting Hilbert space as follows 

S j\ I) 3 = = Sf\ T)j , S$1 T>i = U>j . S+\ I) j = it),, s;\i)j =~\ih , S?\ t)j = ^| t>j • (A18) 
We now define a particular set of states in the spin model in the following way |5j] ] 

N 



\k Pl ,...,kp N ) = A (ri...(T N (kp 1 ,...,kp fr )Y[(S 1 

a 1 ...a N =±l t=\ 



,(l-a e )/ 2 



(A19) 



Here we have used conventions where Oj —] corresponds to 1 and Oj =[ corresponds to —1. The set of equations 
(A.14), due to the Schrodinger equation, now induces the following between states in the spin model 

\k Pll kp 2 ) = Y 1 ' 2 (smkp 2 ,smkp 1 )\kp 2 ,kp 1 ) , (A20) 

where the operator Y"- 7 ' is given by 

—U/2i Vi — v 2 



Y^(v u v 2 ) 



v\ — v 2 — U/2i v\ — v 2 — U /2i 
Here / is the identity operator over the Hilbert space of the spin model and n^ ,fe ' is a permutation operator 



4% • s k 



(A21) 



(A22) 



Here S^' v,z are the spin operators defined in (Al£), where = ± iSi. In order to derive ( A20) one simply uses 



the explicit form (A15) of the S-matrix and the following properties of permutation operators 

n^io) = |o> , nt^srnt 1 ' 2 ) = , n^s^n^ = . 



(A23) 



The Y-operators (which are related to the S-matrix (A1E ) via multiplication by a permutation operator) were first 
introduced by C.N. Yang in his seminal 1967 paper p2L 



On the level of the auxiliary spin model the periodic boundary conditions (A17) translate into 

\k Pl ,kp 2 ) =exp(ik Pl L)IlW\kp 2 ,kp 1 ) . 



With the use of (A20) this is then transformed into 

\k Pl ,kp 2 ) — exp (ikp 1 L) X ' (sin k Pl , sin kp 2 )\kp 1 , kp 2 ) , 
where the operator X^ ,k is given by 



(A24) 
(A25) 
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X^ k (v 1 ,v 2 ) =ILWYi> k (v x ,v 2 ) = 



-U/2i 



v\ — V2 — U/2i 



vi - v 2 



v\ — v 2 — U /2i 



(A26) 



In other words, X 1,2 is the S-matrix ( Ala ) viewed as an operator in the auxiliary spin model. We now make the 
crucial observation that the operator X 1 ' 2 is precisely the transfer matrix of an inhomogeneous spin-1/2 Heisenberg 
model on a 2-site lattice 



X ' (sin kp 1 , sin kp 2 ) — r(sin kp ± \{suikp 1 , sin kp 2 }) 



(A27) 



where r is given by ( BIO ) with N = 2. This is shown in full generality in (B16). 

The periodic boundary conditions (A17) can thus be rewritten as an eigenvalue problem for the transfer matrix 
t(A = sin kp t |{sin kp t , sin kp 2 }) of an inhomogeneous Heisenberg model on 2-site lattice 



\kp, , kp„) = 



iLk 



Pl r(sinfcp 1 |{sinfcp 1 ,sinfcp 2 }) \kp 1 ,kp 2 



(A28) 

The diagonalization of the transfer matrix t(A| {sin kp 17 sin kp 2 }) is carried out in Appendix^. From the point of 
view of constructing eigenstates of the Hubbard Hamiltonian with periodic boundary conditions, we have succeeded in 
reducing the problem to a much simpler one, namely diagonalizing the transfer matrix of an inhomogeneous Heisenberg 
model. This is the essence of C.N. Yang's nested Bethe Ansatz procedure. 

For the problem at hand we need to distinguish two cases, depending on the spins of the two electrons. 

• two electrons with spin up 

Here the appropriate eigenstate of r(A|{sin fop, , sin kp 9 }) is found in the sector with no overturned spins, i.e. it 
is the ferromagnetic state with all spins up. From ( B15 ) we see that its eigenvalue is equal to 1. The periodic 
boundary conditions for the case of two electrons with spin up thus take the simple form 



JkjL _ ^ 



J = 1,2. 



(A29) 



These are precisely the Lieb-Wu equations (0),(||) for the case N = 2, M = 0. Using that the appropriate 
eigenstate of r(A|{sinfcp 1 ,sinfcp 2 }) is equal to |0), we infer by comparison with (A19) that all amplitudes are 
equal to 1. This then implies the following explicit form for the wave function in sector Q 



ip(xi,x 2 ;(Ti,<J2) = ^ si g n (-PQ) ex P(^ ^ kpjXQj) 
Pes 2 j=i 



(A30) 



This agrees with Woynarovich's result 

one electron with spin up, one with spin down 

Now the appropriate eigenstate of r is found in the sector with one overturned spin. Using dB15| ) we obtain the 
corresponding eigenvalue 



1 - iU/[2(smk Pl - Ai + iU/4)] 



sin k 



Pi 



Ai - iU/4 



sin kp 1 — Ai + iU / 4 ' 

where Ai fulfills the Bethe Ansatz equations of the inhomogeneous Heisenberg model on a 2-site lattice 



(A31) 



Ai — sin kj + iU /4 
L J j Ai — sin kj — iU /4 



n 



(A32) 



Inserting (A31) into (A28) we obtain the following quantization conditions for the momenta kj due to periodic 
boundary conditions 



exp(iLA;p 1 ) 



Ai — sin kp 1 — iU/ '4 



forP e 



Ai - sin k Pl + iU/4 

Equations (A32) and (A33) coincide with the Lieb-Wu equations (0),(||) for the case N = 2 and M = 1. 



(A33) 
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Let us also determine an explicit expression for the amplitudes A (71(T2 (ki,k2)- Comparing the result (B35) for 
the eigenstate of r(sinfcp 1 |{sinfcp 1 ,sinfcp 2 }) with ( A19) we see that 



A^ 1(72 (kp 1 ,kp 2 ) = \ \ 



Ai - sin k Pj - ^ 



=l" V Al ~ Shl fcp 3 + X J A l - Sm k Py + X 



it/ ' 



(A34) 



where y is the position of the down spin in the sequence cr±a 2 . The wave functions ( [A9| ) with amplitudes (A34) 
coincide with Woynarovich's result (H) for the case N = 2, M = 1. 



3. Three electrons 



Let us now explicitly construct the wave function for the case of three electrons N = 3. The Schrodinger equation 
is of the form 

- ^2 [ip(xi + s,x 2 ,x 3 ;a 1 ,a 2 ,cr 3 ) + ip(xi,x 2 + s, x 3 ; <7i, a 2) 03) + ip(xi,x 2 ,x 3 + s; <7i, a 2) cr 3 )] 



s=±l 



-?7 ^^(xjjXfe)^^!,^,^;^!, 0-2,0-3) = - —r)^^!, ^2, 2:3; 0-1, 0-2,0-3) • 



3^7 t/L, 
~2 4~ 



(A35) 



In the sector defined by xq 1 < xq 2 < xq 3 , where Q is a permutation of the labels {1,2,3}, the Bethe Ansatz 
wavefunction reads 



ip(xi, x 2 , x 3 ; ctx, cr 2 , 0-3) = ^2 sign(PQ)A aQiaQ2aQ3 (k Pl ,kp 2 ,kp 3 )cxp(i^2k P: jXQ j ) . (A36) 

Pes 3 j=i 



Substituting ( |A36j ) into ( A35 ) for the case x% ^ x 2 7^ x 3 ^ x\ we obtain the following expression for the energies 

3U UL 



E = — (2 cos k\ + 2 cos k 2 + 2 cos £3) — 



(A37) 



Single-valuedness of the wave-function now leads to a larger number of relations between the amplitudes. We have 
to consider the three cases ip(x, x, x 3 ; a%, a 2 , cr 3 ), ip(x,x 2 ,x;ai,a 2 ,a 3 ), ip(x\,x, x; a±, cr 2 , <r 3 ). A simple calculation 
yields the following conditions on the amplitudes 

A °Ql °Q 2 °Q 3 ( k Pi > k P2 > k Ps ) - Ar Ql a Q2 a Q3 {k P[ , ^ , k P > ) = A a a a (fcp, , kp^ , fc P ^ ) - A (fc Pl , /cp 2 , fcp 3 ) , 



(A38) 



where Q and P are arbitrary permutations of {1, 2, 3} and Q' — Q{j,j + 1), P' = P(j,j + 1)- Our notations are such 
that for any permutation of N elements S — (Si, . . . , Sjv) 



+ 1) — (5*1, ■ ■ ■ , Sj-i, Sj+i, Sj, Sj +2 , . . . , Sjv). 



(A39) 



Let us consider a specific example of (A38) in more detail. The wave-function i/}(x, x, x 3 ; 01, a 2 , a 3 ) with x 3 > x can 
be expressed alternatively in sector Q = (1, 2, 3) or in sector Q' = (2, 1, 3). Equating the respective expressions ( A36) 
and using the orthogonality property of plane waves we obtain 

^ai(Tj(T3 ( k Pl ' k Pl 1 kPa) ~ <j 2 <j 3 (kp 2 , kp 1 , kp 3 ) — AaQ 2 a Ql aQ 3 {kp 2 , kp 1 , fcp 3 ) — A <rQ 2 <r Ql a Q3 {kp 1 , fcp 2 , fcp 3 ) . (A40) 



This indeed coincides with the general result (A38). 

We note that the Bethe ansatz wave function (A36) by construction is antisymmetric under simultaneous exchange 
of spin and space variables. It is easy to see that this fact assures the Schrodinger equation (A35) to be satisfied, 
when the three electron are occupying the same site. Moreover, this reasoning readily generalizes to the case of an 
arbitrary number of electrons in the following subsection. The only non-trivial case left to consider is the case of two 
electrons at the same site. 

Let us start with x\ — x 2 — x < x 3 . Using (A36) in (A35) we obtain the following condition on the amplitudes 
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[e-** + ] A CT1CT2CT3 (k Pl , k P . 2 ,k P3 ) + [e~^ 2 + ^ ]A aitT2a3 (fcp 2 , k Pl , kp 



+U [A ai (kp 1 ,kp 2 ,k P3 ) -A ai (kp 2 ,k Pl ,kp 3 )] 



, (kp 2 , k Pl , kp 3 ) 



= -2(cosfc Pl + cosfcpj [A (7ia2Cr3 (kp 1 ,kp 2 ,kp 3 ) - A ai(T2lT3 (kp 2 ,kp 1 ,kp 3 )] . 



(A41) 



Now making use of the fact that (A41) and ( |A38D are of the same structure as ( Al2| ) and (All), we conclude that 
the following relation between amplitudes holds 



At 2 £ti<x3 (kp 2 , kp 1 , kp 3 ) — £>a 2 T2 Q i p l , kp 2 ) A TlT2C j 3 (kp 1 , kp 2 , /cp 3 ) 



(A42) 



where S is the two-particle S-matrix (A15). All other cases of coinciding coordinates can be analysed in exactly the 
same way. The final result is 



A ^Q' 2 ^ Q ' 3 ( k pi' k p^ k n) = £ sl^ + \ T2 {kp v k Pj+1 ) A aQi ... aQj i 



TlT 2 <?Q j + 2 ---<JQ 3 



(kp 1 ,kp 2 ,kp 3 ) 



(A43) 



where Q and P are arbitrary permutations and Q' = Q(J,j + 1),P' = P(j,j + 1)- Equation ( A43| ) has important 
consequences. In order to exhibit these clearly, it is convenient to express ( A43| ) in the framework of the auxiliary 
spin model introduced above. Inserting (A43) into (A19) we obtain 



\kp>,kp^,kp^) = r- 7J+1 (sinfcp 3 ,sinfcp 3+1 )|fcp 1 ,fcp 2 ,A:p3} , 



(A44) 



where Y^ k is given by ( A21 ). There are altogether (N — 1)N\ = 12 equations ( A44 ) and all of them have to be 
consistent with one another! This puts severe constraints on the operators Y^ k . 

Let us explain this in more detail. We recall that the symmetric group Sn is generated by the identity id and the 
transpositions of nearest neigbours (j, j + 1), j = 1, . . . , N — 1, modulo the relations 



U,J + + 1, j + 2)(i,j + 1) = (j + 1,3 + 2)0', j + l)(j + 1, 3 + 2) , 

(j,j + l)(k,k+l) = (k,k + l)(j,j + l) for |j - fc| > 1 , 

(j'.i + l)(7,j + l) = id • 



(A45) 
(A46) 
(A47) 



Equation (A45) is called the braid relation. Note that (A46) is non-trivial only for N > 3. 

By inspection of ( All ) we see that the F-operators act on the states \ki,k2,k^) by exchanging neighbouring 
components of the vector k = (fci, fe, £3), that determines the state \k±, &2, fe) of our auxiliary spin system. Hence all 
states \kp 1 ,kp 2 ,kp 3 ), P G S3, can be obtained from \ki,k2,k s ) by repeated use of ( A44 ). Equivalently, ( A44 ) allows 
us to obtain the state corresponding to any permutation P G S3 from a state corresponding to any other permutation 

PeS 3 . 

It follows from ( A45| )-( A47) that a representation of a permutation as a product of transpositions of nearest neigh- 
bours is not unique. For the case at hand, N = 3, we have for instance, (1, 3) = (1, 2)(2, 3)(1, 2) = (2, 3)(1, 2)(2, 3). We 
are thus facing a consistency problem for equations (A44): The relations ( A45| )-( A47) impose consistency conditions 
on the y-operators. 

Let us study these consistency conditions. For N = 3 we only have to consider ( A45| ) for the case j = 1 and (A47). 
Thus QA45D implies that 



\kp i , kp 2 , kp ) — Y 1,2 (sin kp 2 , sinfcp 3 )F 2 ' 3 (sin kp ± , sin fcp 3 )y 1,2 (sin kp x , sinfep 2 )|fep 1 ,kp 27 kp 3 ) , 
= Y 2,3 (sin fcpj , sin kp 2 )Y 1 ' 2 (sa\ kp 1 , sin fcp 3 )l" 2 ' 3 (sin fcp 2 , sinfcp 3 )|fcp 1 , kp 2 , kp 3 ), 



(A48) 



where P = P(l,3) = P(l, 2) (2, 3)(1, 2) = P(2,3)(l,2)(2,3) = (P 3 ,P 2 ,Pi). Assuming \ki,k 2 ,k 3 ) to be arbitrary, we 
conclude that 



Y 1 ' 2 ^, s 3 )Y 2 ' 3 ( Sl , S3 )Y^ 2 ( Sl , S2 ) = Y 2 ' 3 ( Sl , s 2 )Y 1 > 2 (s 1 ,s 3 )Y 2 > 3 (s 2 ,s 3 ) , 



(A49) 



where Si,2,3 are arbitrary complex numbers. This is the famous Yang-Baxter equation. It is now crucial that (A4E) 
can be verified by direct calculation. Therefore (A44) is consistent with (A45). 



56 



Similarly, the use of equation (A47) in ( A44j ) yields the requirement 

(y H ' +1 («,u)) _1 = Y j ^ +1 (v,u) 



(A50) 



It is easy to verify by direct calculation that the operators Y^ k defined in(A21) indeed fulfil (A50). Hence (A44) is 
consistent with (A47), and we conclude that the entire set of equations (A44) is consistent. 

Note that the above considerations easily generalize to the case of arbitrary TV, which will be treated in the next 
subsection. In addition to ( A45| ) and (A47) we then will have to consider (A46) which leads to the condition 



Yj,j+i(si, s 2 )Y k ,k+i(s3, Si) = Yfc, fc+ i(s 3 , s 4 )Yj J+ i(si, s 2 ) for \j - k\ > 1 

which is trivially satisfied. 

Let us now impose periodic boundary conditions on the wave function 



(A51) 



?/> (0,2:2,2:3; o-i, 0-2,0-3) = tp(L, x 2 , x 3 ; 0-1, 0-2,0-3) , 

ip(l,x 2 , 2:3; 0-1, CT 2 , 0-3) = tp(L + 1,2:2, 2:3; 0-1, CT 2 , 0-3) , 

ip(xi, 0,2:3; 0-1,0-2, 0-3) = -0(2:1, £,2:3; 0-1, 02, 0-3) , 

ipfeu 1,^3! 01, 0-2, 0-3) =ip(xi,L + l,x 3 ;ai,a 2 ,a 3 ) , 

ip(xi, x 2 , 0; 0-1, 0-2, 0-3) = -0(2:1, 2; 2 , i; 0-1,0-2, 0-3) , 

0(2;i, 2:2, 1; 0-1,0-2,0-3) = 0(2:1, x 2 , L + 1; 0-1,0-2, 0-3) . 
Inserting ( |A36| ) into ( |A52D yields 

A <?Q 1 <?Q 2 <?Q ;i ( k Pn k P2, k P 3 ) = e xp(ik Pl L)A aQ2aQ3<TQl (kp 21 kp 3 ,kp 1 ) , 
where Q and P are arbitrary permutations of {1, 2, 3}. In terms of the auxiliary spin model ( A53 ) is expressed as 

\k Pl , kp 2 , kp 3 ) = expiikp.Ljn^n^ikp^kp.^p,) 

= exp(ikp 1 L)n {1 ' 2) n {2 ' 3) Y 2 ' 3 (smkp 1 ,smkp 3 )Y 1 ' 2 (smkp 1 ,smkp 2 )\k Pl ,kp 2 ,kp 3 ) 
= exp(?;/cp 1 L)n (1 ' 2) X 2 ' 3 (sinfcp 1 ,sinfcp 3 )n (1 ^ 2) A 1 ' 2 (sinfcp 1 ,sinA:pJ|fcp 1 ,A:p 2 ,fcp3) 
= cxp(z/cp 1 L)A 1:3 (sin kp 1 , sin k p 3 ) X 1,2 (sin kp ± , sin kp 2 )\kp 1 , kp 2 , fcp 3 ) , 



(A52) 
(A53) 



(A54) 



where we have used the identities (B17). Using (B16|) we now obtain in complete analogy with the 2-electron case 



\k Pl ,kp 2 ,kp 3 ) = e % p i T(sinkp 1 \{sink Pj ;j = l 1 ... 1 3})\k Pl ,kp 2 ,kp 3 ) . 



(A55) 



We again can use the results for the diagonalization of the inhomogeneous transfer matrix r(sin kp 1 \ {sin kp. ; j = 
1, . . . ,3}) derived in Appendix [B|. We need to distinguish two cases, depending on the spins of the two electrons 
(recall that we consider only states for which the number of down spins not larger than the number of up spins, all 
other states are obtained by using the spin-reversal symmetry). 

• Three electrons with spin up 

Here the appropriate eigenvector of r(A|{sinfcp j ;j = 1, . . . ,3}) is the ferromagnetic state, with eigenvalue 1. 
The periodic boundary conditions for the case of three electrons with spin up thus take the form 



.7 = 1,2,3. 



(A56) 



These are precisely the Lieb-Wu equations (@),(||) for the case N = 3, M = 0. 
Using ( A19| ) we see that all amplitudes are trivial 



t (k Pl ,kp 2 ,kp 3 ) = 1 . 



(A57) 



The corresponding wave- function (A36) coincides with Woynarovich's result 
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• two electrons with spin up, one with spin down 



Now the appropriate eigenvector of t is found in the sector with one overturned spin. From (B15) its eigenvalue 
is given by 



TT ;r„/ . ■ > j j \ n SU1 kp, — Ai — lU I 4 

1 - iU 2 (smk Pl - Ax + iU A)] = / , 



(A58) 



where Ai fulfills 



3=1 



Xi - sin fej + iU/4 
Ai — sin kj — iU /4 



(A59) 



Inserting reval into (A55) we obtain the following quantization conditions due to periodic boundary conditions 

Ai — sinfcpj — iU/A 



exp(iL/cp 1 ) = 



Ai — sin fcpj + iU/4 



forP G S3 



(A60) 



Equations ( A59| ) and ( A60| ) are precisely the Lieb-Wu equations for the case TV = 3 and M = 1. 

An explicit expression for the amplitudes again is obtained from ( ]A1S| ) and (B35), with the result 



y-i 



ArLCT2<73 {kPl 1 , kp 3 ) — J | 

3 = 1 



Ai — sin kp. 



iU 



Ai -sinfc P . +f 



Ai-sinfcp H +f 



(A61) 



where j/ is the position of the down spin in the sequence a\02<Jz- Inserting ( A61 ) into ( A36 ) we obtain (^) for 
the case N = 3, M = 1. 



4. AT electrons 



It is now clear how to generalize the above results to the case of N electrons. The Bethe Ansatz for the solution rp 



of the Schrodiger equation (A7) in the sector Q with xq 1 < xq 2 < . . . < xq n is 



N 



rp(xi,...,x N ;ai,...,(7 N )= sign(PQ)A aQi ... CTQjv (k Pl , . . . , kp N ) exp(i k Pj x Qj ) . (A62) 

Pes N 3=1 



Substituting (A62) into (A7) for the case Xj ^ Xk, (j,k = 1, .. . , TV; j ^ k) we obtain the following expression for the 
energies 



N 



E = — 2^ cosfcj 



iVE7 C/L 



(A63) 



3=1 



Using the single valuedness of the wave function and solving the matching conditions at the sector boundaries i.e. the 
Schrodinger equation for the cases where two of the coordinates coincide, we obtain the following set of equations 



-4, 



a Q' 1 - CT c}' N (^ p i'' ■ ' ■ ' ~~ ^ CT Q3 + i T 2( fcp j'^ p j + i) Aa Ql ... (TQ ._ lTl T2aQ :j+2 ...c?Q N (kp 1 , . . . ,kp N ) 



(A64) 



where Q and P are arbitrary permutations and Q' — Q(j, j + 1),P' = P(j, j + 1)- In terms of the auxiliary spin model 
(|A64[) reads 



\kpi,- ■ -,k P ' N ) = Y J (sin k Pj , sin k Pj+1 )\k Pl , . . . ,kp N ) 



(A65) 



The mutual consistency of equations ( |A65| ) follows from ( A49 ) and (A50). We now impose periodic boundary condi- 
tions on the wave function 
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i/>(xi, . . . ,Xj-i,0,Xj+i, . ■ ■ ,Xn;cti, . ■ . ,ctjv) = ip{xi, . ■ ■ ,Xj-i,L,Xj+i, . . . ,x N ;ai, . . . ,<r N ) , 
tp(xx, . . . ,Xj-i,l, Xj+i, . . . ,x N ; (7i, ... , a N ) = ip(xi, . . . ,Xj-i,L + 1, Xj+%, . . . ,x N ; <7i, . . . , a N ) 
where j — 1, . . . , N. Inserting ( |A62 ) into ( A66| ) yields 



A "q x -tq n ( k Pn ■ ■ ■ > k P N ) = exp(ik Pl L)A 



CQ 2 --- cr Q N ' T Qi ( kp z ' ' • ■ ' Kp N ' 



k Pi) , 

where Q, P G Sm are arbitrary. In terms of the auxiliary spin model ( |A67| ) is expressed as 

\k Pl , . . . , kp N ) = exp(ifcp 1 i)n( 1 ' 2 )n( 2 ' 3 ' . ..n^-^ifcp,,. . . , k PN , k Pl ) 



(A66) 



(A67) 



exp^L)^ 2 )^ 3 )...^- 1 ^) 

_m— 

-l,Nt- ;„_ u- wi,N-i 



'jV-2 



-1,N- 



'(sinfcpj, sin kp N 



^Pn) 



exp(ik Pl L)X • (smkp 1 ,smkp N )X ' (sinkp 1 ,smkp N _ 1 ) . . 
x . . . AT 1,3 (sin fcpj , sin fcp 3 )X ' (sin kp x , sin Jfc.pi, ) | fc^ , . . . , kp N ) 



(A68) 



where we have used the identities (B17). Using (B16) we now obtain in complete analogy with the 2 and 3 electron 

cases 



Ife 



Pi: 



k P N )=e Pl T{wnk Pl \{8mkp j ;j = l,...,N}) \k Pl , . . . , fc Pjv ) 



(A69) 



Next we again use the results for the diagonalization of the inhomogeneous transfer matrix r(sin kp t |{sin kp. ; j — 
1, . . . , N}) derived in Appendix ||. We now need to distinguish [N/2] + l cases ([x] is the integer part of x), correspond- 
ing to the possible values of M. In the sector with M down spins the eigenvalue of r(sin kp 1 |{sin kp j ; j = 1, . . . , N}) 
is given by (B15) 



M 

n 

3=1 



sin fcpj — Aj — «/7 /4 
+ i[//4 



sin k 



Pi 



(A70) 



where the rapidities A, fulfill 



TV 



M 



n 

fc#3" 



A,- — \i 



iU/2 



n\j — sin ki + ill/ 4 
^ Xj - sin ki — ill /A ^ A, - A* - iU/2 



(A71) 



Inserting this result into ( A69| ) we finally obtain 



exp(iLfc Pl ) = Yl 



Xj — sin kp 1 — iU /4 
Aj — sin fepj + zt//4 



forP e 5jv 



(A72) 



Equations (A71) and (A72) are precisely the Lieb-Wu equations 



1) for the case of N electrons, M of which have 



spin down. In order to obtain an explicit expression for the amplitudes we now need to make use of the general result 
( B18) for cigcnstates of the transfer matrix r(sin kp 1 1 {sin kp. ; j = 1, . . . , N}) of the inhomogeneous Heisenberg model. 
We find 



M 



A ai ... aN (k Pl ,...,kp N ) = A ^Y[ 

ir£S M 



1 



Vt-l 

n 



X 7Tt — sin k p s 



ML 

4 



iU 



t=l (. ~ Sin kp yt + V s=l ~ Sin k P* + "4 



(A73) 



where 1 < y 1 < y 2 < ■ ■ ■ < y M < N are the positions of the down spins in the sequence a± . . . crjv and is given by 



Ajj 



n 



Att, — Att^ 



2 



A^i A Wfc 



(A74) 



l<i<fc<A/ 

The resulting explicit expression for the wave function (A62) coincides with Woynarovich's result (H). 
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APPENDIX B: INHOMOGENEOUS HEISENBERG MODEL 



1. Algebraic Bethe Ansatz 

Our starting point is a lattice of N spin- 1/2 's. The corresponding Hilbert space is V\ ® Va 
isomorphic to C 2 . We define the Pauli matrices f = (t x ,t v ,t z ) by 




-i 

1 



1 

-1 



® Vn , where V, is 



(Bl) 



and 



(t x ±ir»). 



The central object of the Quantum Inverse Scattering Method is the i?-matrix, which is a solution of the Yang- 
Baxter equation. For the case of the spin- 1/2 Heisenberg model it is of the form 



R(X,fj.) 



//(At, A) 

g(jjt,\) 1 

1 g(n,\) 

V /(a*, A). 



(B2) 



where 



/(M) = i 



2(/z-A) 



9 , (i">A) 



it/ 



2(A* - A) 



(B3) 



R(X, fj.) acts on the tensor product space Vb <8> Vb, where Vq is isomorphic to C 2 . The Yang-Baxter equation for R is 
an equation on the space Vb <8> Vb <S> Vb and can be written as 



R 23 (X,li)R 12 (X,iy)R 23 (^iy) = R 12 (p,is)R 23 (\,v)R 12 (\, n) , 



(B4) 



where the superscript indicates in which spaces the R matrix acts nontrivially. We now define an L-operator acting 
on the tensor product between a "matrix-space" Vq and a "quantum-space" V n , which is identified with the Hilbert 
space over the n th site of our lattice of spins, by 



L n (X) 



X 



-I 



iU/2 



n (o.i 



A + iU/2 X + iU/2 

1 I X + (1 + KW* r~ iU/2 
X + iU/2 I r+iU/2 X + (1 - T*)iU/A 



The Yang-Baxter equation (B4) implies the following intertwining relations for the L-operator 

R(X, /i) (L n (X) (g) L n (fi)) = (L n (fi) ® L„(A)) i?(A, a*) , 



(B5) 



(B6) 



where the tensor product is between matrix spaces, i.e. (B6) is a relation over the space Vq ® Vb <8> V n . 

Next we note, that the intertwiner for the L-operator (B6) still holds, if we shift both spectral parameters A and fj, 
by an arbitrary amount v n , i.e. 



#(A,At) (L„(A - v n ) ® L„(a* - v n )) = {L n {(i - u n ) ® L„(A - u n )) R{X,fj) . 
We now construct an inhomogeneous monodromy matrix in the following way 

T(n\{a :j })=L N (n-a N )L N - 1 (n-a N - 1 )...L 1 (n-a 1 )= ^(a*)) 



(B7) 



(B8) 



Here a\, . . . , aw are N arbitrary complex constants. The intertwiner (B7) can be lifted to the level of the monodromy 
matrix 



i?(A lA t) (T(\\{ aj })®T(vL\{ aj })) = (T^\{ aj })®T(X\{a 3 }))R(X^) 



(B9) 



GO 



By tracing ( p38| ) over the matrix space Vq one then finds that the transfer matrices 

r(fx\{ aj }) = Tr (T( M |{ % })) = A(jjl) + £>(//) 



(BIO) 



commute for any values of spectral parameter fi, i.e. [r(/i|{aj ■}), T(i/|{aj})] = 0. That implies that the transfer matrix 
is the generating functional of an infinite number of mutually commuting conserved quantum operators (via expansion 
in powers of /i). 

Eigenstates of the transfer matrix are constructed by means of the Algebraic Bethe Ansatz. Starting point is the 
choice of a reference state, which is a trivial eigenstate of t(/i|{oj}). We choose the saturated ferromagnetic state 



|0> = |m...T> = ®£=ilT>n 



(Bll) 



The action of the L-operator (B5) on | |}„ can be easily calculated and implies the following actions of the matrix 
elements of the monodromy matrix 

A(ji)\Q) =a(/x)|0) , o(m) = 1, 



N 



D(jj)\Q) = d(/i)|0) , rf(/i) = H 



)I — Oj 



3=1 



ji — cij + iU/2 



C(/x)|0} = 0, B(p)\0)^0. 



(B12) 



From (B12) we see that B(X) play the role of creation operators. Acting with B corresponds to flipping a spin. 
States with M down spins can be constructed as 



M 



F{\x, A M ) - Y[(-2i/U)B(Xj - iU/4)\0) , 

3=1 



(B13) 

where we have shifted the spectral parameters and introduced a particular normalization for later convenience. The 
requirement that the states (B13) ought to be eigenstates of the transfer matrix puts constraints on the values A n : 
the set {Xj} must be a solution of the following system of Bethe Ansatz equations 



i r Xj -a k + iUjA _ t-j Xj — A; + iU/2 

11 Xj - o fc - iU/4 11 A, - Xi - iU/2 ' 3 ' 

*#3 

The corresponding eigenvalues of the transfer matrix are 

M M 



(B14) 



^\{a 3 }) F(X 1 ,...,X M ) = 



oo*) n /(a*. a ; - if/ / 4 ) + ^ n ^ - if/ / 4 > ^) 

3=1 3=1 



F(Xi, . . . , Am) 



(B15) 



For our present purposes we need to consider the transfer matrix evaluated at the first inhomogeneity. We find 
t(oi|{Oj}) = Tr L N (ai - a N )L N -i(ai - O/v-i) . . . L 2 (a 1 - a 2 )n (0,1) 



= Tr f 



= Tr f 



n(°. 1 )n(°' 1 )L^(a 1 - a N )n^n^L N ^( ai - ajv-On^^n^ 1 ) . . . n^L 2 ( ai - a 2 )n(°- 1 ) 

IlfMX 1 '"^, arfX^-Hai, a N -i) . . . X l *(ai, a 2 ) 



X*- N {ax, a^X^-^ai, a N ^) . . . X^( ai ,a 2 ) , 



-l,iV-l 



(B16) 



where X h is defined in (A26). Here we have first used the explicit form of the L-operator (B5) and then the identities 



jj(3.fc)jj(3»jj(3.fc) _ n(k,n) ^ n(j> k )jjU>k) — j 



(B17) 
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2. Explicit expressions for the eigenstates 



In this subsection we derive the following expression for the eigenstates ( BIS ) 

M ( 



1 



{•y.}^GS flf t=l r Wt ° P w + 4 s =l + 4 J j=l 



II ^!?T5 IRi°> 



(B18) 



where the summation extends over 1 < y 1 < y 2 < y 3 < ... < y M < N, P is a permutation of iV elements, and the A Tl 
is given by 



A — A —EL 

i<;<fe<M ' fc 



(B19) 



Our discussion closely parallels J51J . The basic tool for proving ( p318| ) is the "n-site generalized model" J96|,[)7j (sec 
also @ p.l51f and p.171). 

Generalised two-site model. 



Let us first divide the product of L-operators in (B8) into two parts 

T J (A) = J L n (A-a n )... J L 2 (A-a 2 ) Li(A-ai) , 
^xj(A) = L N (A - a N ) . . .L n+2 (A - a n+2 ) L n+ i(A - a n+ i) 

Clearly we have 

T(A) = T„(A) T,(A) . 
Both T n (A) and T I (A) are 2 x 2 matrices 



r,(A) 



A, (A) B,(A) 
C X (A) D,(A) 



r H (A) 



A„(A) B u (A) 
C„(A) D U (A) 



(B20) 



(B21) 



(B22) 



By construction the matrix elements of T f (A) commute with the matrix elements of T u (A) . The commutation 
relations of the matrix elements of the same T-operator are as in (B9) 



i?(Ai,A 2 )(T Q (Ai)®T a (A 2 )) - (r Q (A 2 )®T Q (Ai))i?(Ai,A 2 ) , a = 1,11. 
The matrix elements of T can be expressed in terms of the matrix elements of T 1 and T XI , e.g., 

B(A) = A U (A) B,(A) + B U (A) D T (A) . 



(B23) 



(B24) 



It is also possible to express the vectors Yln=i B{Aj) |0) in terms of Bj and B u . In order to do this we will use that 



C a (A)|0) = , A Q (A)|0) =o Q (A)|0) , £> Q (A)|0) = d Q (A)|0} , a = 1,11. 



Here 



»,(A) = 1, d I (A) = f[ 



A — a,- 



A - aj + iU/2 ' 



(B25) 



(B26) 



and 



N 



(B27) 



In pflM it was proved that 



G2 



M 



n s ( A i)io>=E n n ^(^^(^/(a^ob^ob^a^o) 



(B28) 



On the right hand side of ( |B28 ) we have summations with respect to partitions of the set of all Aj 's into two subsets 
Sj = {Aj,} and <S/j = {A^f}. Here k labels different A in the subset Si and m labels different A in the subset 5//. 
The above model is called generalised two-site model because T is represented as a product of two factors. This is 
not sufficient for our purposes however. Let us therefore now consider the so-called 



Generalised fc-site model. 



Let us represent T as a product of k factors 

T(A) = T fc (A) ■ • • T 2 (A) Tx(A) 



(B29) 



Here each T a (A) is a string of L-operators like in (B20). The commutation relations of the matrix elements of each 
of the T a (A) is given by the same intertwiner as in (B9). 

The matrix elements of Tp (A) commute with the matrix elements of T a (A) if a ^ /3 and like for the two-site model 
we have 



T a (A) 



A* (A) B a (A) 
C a (A) D a (A) 



(B30) 



C a (A)|0)=0 , A Q (A)|0) =a Q (A)|0) , D a (A)\0) = d a (A)\0) . (B31) 

By explicitly multiplying the matrices in ( |B29| ) we can express B(A) in terms of matrix elements of the T a (A). 
Iteration of (B2S) leads to the following expression for the eigenf unctions of the transfer matrix T(A) 



M I k 

HB(A j )\o)= e n n s .( a »ji o > 

J=l S lt ...,S h \a=l A ma eS a 



n n n M^j^Kp/^Ay . 



(B32) 



1<a</3<fc A m „e5c, A fe „£cS fl 



Here the summation is with respect to the partitions of the set of all Aj's into k disjoint subsets Sp, f3 = 1, . . . , k. 
The index m a enumerates different A in the subset S a and the index kp enumerates different A in the subset Sp. 



Equation ( B32 ) was first proved in 97|. 

We now consider the special case of the generalized A^-site model, where N is the length of the underlying lattice. 
This means that ea ch fac tor T a (A) in our generalised Af-site model is identified with an individual L-operator in ( |B^ ) . 
The eigenvalues in ( B31 ) are given by 



o Q (A) = 1 , d a (A) = 



A - a a + ill/2 



The i3-operators are given by 



B a (A) = 



iU/2 



A - a a + iU/2 



(B33) 



(B34) 



and have the important feature that -B a (Ai)£? Q (A2) = 0. This implies that each set S a in ( B32 ) consi sts of maximally 
one element. We are now in the position to write down explicit expressions for the eigenstates ( |B13|) . 

a) One overturned spin. 



Let us first consider the eigenfunctions in the sector with one overturned spin. Application of (p32j) yields 
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JV 



F(A) = (-2i/U)B(A iU/A)\0) = £ ^ |0) 



Here all sets <S Q in (B32) except one are empty. This one (one-element) set is S y = {A — iU/A}. 
a) Two overturned spins. 



Now application of (B32) leads to 



F(Ax,A 2 )= J2 E^^l°) 



1 



1 



l<i/ 1 <y 2 <iV7reS2 



n 



A -a--^> /y2_1 



= 1 



Ajri — <V +iU/4 A W2 — a y2 + ill/ A 

f(A^,A„ 2 ) 



3 = 1 ">■ 3 i J \ 1=1 Il7r 2 

Here 7r is a permutation of two elements 1, 2 and 



A „ lU 

W A T2 - a, + f 



f(A vl ,A n% ) = 



A - A - — 



A-n-i A^ 



(B35) 



(B36) 



(B37) 



In this case only two subsets 5 a are nonempty. Each of them consists of one element (S yi = {A 7ri — iU/A} and 
S V2 = {Air 2 — iU/A}). Equation (B36) is of the desired form (B18) if we identify 



A - A —EL 



A-ji 1 A<x 2 



(B38) 



which is in complete agreement with (B19). 
c) M overturned spins. 



The result for 2 overturned spins generalizes straightforwardly to M overturned spins. The nonempty subsets S a in 



= A„ 



( B32 ) (each of which consists of exactly one element) are S Vl = A Xl — ill /A, S V2 = A 7I2 —i U/A, . 
where it is some permutation of M elements. A straightforward calculation then yields ( B15 ) and ( B19 ). 



■iU/A, 
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APPENDIX C: THE SPECTRUM OF THREE ELECTRONS WITH ONE-DOWN SPIN FOR L = 6 



We present a complete list of eigenstates of the Hubbard Hamiltonian (119) in section V for the case N — 3 and 
M = 1 for a 6-site system (L = 6) and U = 5. The energy levels are listed in increasing order. 

The energy eigenvalues obtained by direct numerical diagonalization of the Hamiltonian (the Houscholder-QR 
method) and by the Bethe ansatz method coincide within an error of 0(l(r 15 ). In the table only the energy 
eigenvalues obtained by Bethe ansatz are listed. The last digit for each numerical value has a rounding error. The 
symbols S and P denote the spin and momentum of the eigenstate, respectively. 

There are 90 eigenstates for N = 3 and M = 1, as we enumerated in section V.C.3. The numerical results shown 
in the table confirm the completeness of the Bethe ansatz. 



No. 


Energy 


s 


P/(tt/3) 


types 


1 


-4.19862084914891 


1/2 


5 


real Ij = 0.5, -0.5, -1.5, J = 0.5 


k 3 = 0.542662224387082, -0.315837723840216, -1.27402205174346, A a = 0.587983554411128 


2 


-4.19862084914891 


1/2 


1 


real Ij = 1.5, 0.5, -0.5, J = -0.5 


k 3 = 1.27402205174346, 0.315837723840216, -0.542662224387082, A Q = -0.587983554411128 


3 


-4.00000000000000 


3/2 





quartet Ij — 1, 0, —1 


4 


-3.35402752146807 


1/2 


4 


real Ij = 0.5, -0.5, -1.5, J = -0.5 


k,j = 0.231206350487730, -0.682137186480516, -1.64346426640041, A Q = -1.27426628748753 


5 


-3.35402752146807 


1/2 


2 


real/j = 1.5, 0.5, -0.5, J = 0.5 


kj = 1.64346426640041,0.682137186480516, -0.231206350487730, A Q = 1.27426628748753 


6 


-2.63449090541439 


1/2 





real Ij = 1.5, -0.5, -1.5, J = 0.5 


kj = 1.52991315867874, -0.279871354037999, —1.25004180464074, A Q = 0.845079113394529 


7 


-2.63449090541439 


1/2 





real/j = 1.5,0.5,-1.5,7= -0.5 


kj = 1.25004180464074, 0.279871354037999, -1.52991315867874, A a = -0.845079113394529 


8 


-2.21221243379665 


1/2 


4 


real Ij = 0.5, -0.5, -2.5, J = 0.5 


k 3 = 0.554943548673639, -0.307417851034944, -2.34192080003189, A Q = 0.644785906059544 


9 


-2.21221243379665 


1/2 


2 


real Ij = 2.5, 0.5, -0.5, J = -0.5 


k 3 = 2.34192080003189, 0.307417851034944, -0.554943548673640, A a = -0.644785906059544 


10 


-2.17072407234717 


1/2 


5 


real Ij = 1.5, -0.5, -1.5, J = -0.5 


hj = 1.22622874924835, -0.660511074050219, -1.61291522639472, A a = -1.15790499577381 


11 


-2.17072407234717 


1/2 


1 


really = 1.5, 0.5, -1.5, J = 0.5 


k 3 = 1.61291522639472, 0.660511074050219, -1.22622874924835, A Q = 1.15790499577381 


12 


-2.00000000000000 


3/2 


1 


quartet Ij = 2, 0, —1 


13 


-2.00000000000000 


3/2 


5 


quartet Ij — 1,0,-2 


14 


-2.00000000000000 


3/2 


3 


quartet Ij = 2,1,0 


15 


-2.00000000000000 


3/2 


3 


quartet Ij = 0, —1, —2 



No. 


Energy 


S 


P/(tt/3) 


types 


16 


-1.68312858421806 


1/2 


3 


real Ij = 0.5, -0.5, -2.5, J = -0.5 


k 3 = 0.261468569600159, -0.628328469155396, -2.77473275403456, A Q = -0.993984117066770 


17 


-1.68312858421806 


1/2 


3 


real Ij = 2.5, 0.5, -0.5, J = 0.5 


kj = 2.77473275403456,0.628328469155396, -0.261468569600159, A Q = 0.993984117066770 


18 


-1.00000000000000 


3/2 


2 


quartet Ij = 2, 1, — 1 


19 


-1.00000000000000 


3/2 


4 


quartet Ij = 1, —1, —2 


20 


-1.00000000000000 


3/2 


2 


quartet Ij — 3, 0, —1 


21 


-1.00000000000000 


3/2 


4 


quartet Ij = 3, 1, 


22 


-0.867143471169408 


1/2 


3 


real Ij = 0.5, -1.5, -2.5, J = 0.5 


kj = 0.515010200239397, -1.28773972583561, -2.36886312799358, A a = 0.460329439565352 


23 


-0.867143471169409 


1/2 


3 


real Ij = 2.5, 1.5, -0.5, J = -0.5 


kj = 2.36886312799358, 1.28773972583561, -0.515010200239397, A Q = -0.460329439565352 


24 


-0.763936521983257 


1/2 


2 


real Ij = -0.5, -1.5, -2.5, J = 0.5 


^ = -0.404815536059732, -1.33874801493405, -2.44522665379261, A Q = 0.071452164223840 


25 


-0.763936521983257 


1/2 


4 


real Ij = 2.5, 1.5, 0.5, J = -0.5 


kj = 2.44522665379261, 1.33874801493405, 0.404815536059732, A a = -0.07.1452164223840 


26 


-0.724935621196484 


1/2 


5 


real Ij = 1.5, 0.5, -2.5, J = -0.5 


kj = 1.27598401405461, 0.318718889893781, -2.64190045514499, A a = -0.568959086784317 


27 


-0.724935621196484 


1/2 


1 


real Ij = 2.5, -0.5, -1.5, J = 0.5 


kj = 2.64190045514499, -0.318718889893781, -1.27598401405461, A a = 0.568959086784317 


28 


-0.633335023683704 


1/2 


5 


real Ij = 1.5, -0.5, -2.5, J = 0.5 


kj = 1.54060965024643, -0.274638622892597, -2.31316857855043, A a = 0.886033905628400 


29 


-0.633335023683705 


1/2 


1 


real Ij = 2.5, 0.5, -1.5, J = -0.5 


kj = 2.31316857855043,0.274638622892597, —1.54060965024643, A a = -0.886033905628400 


30 


-0.441363635441783 


1/2 


4 


real Ij = 1.5, -0.5, -2.5, J = -0.5 


kj = 1.24802157741378, -0.602680979199234, -2.73973570060774, A Q = -0.869103625350605 
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No. 


Energy 


s 


P/(tt/3) 


types 


31 


-0.441363635441782 


1/2 


2 


real Ij = 2.5, 0.5, -1.5, J = 0.5 




kj = 2.73973570060774, 0.602680979199234, -1.24802157741378, A a = 0.869103625350605 


32 


-0.164932552352930 


1/2 





real Ij = 1.5, 0.5, -2.5, J = 0.5 




kj = 1.61980144017864, 0.665436431311280, -2.28523787148992, A a = 1.18390417778921 


33 


-0.164932552352929 


1/2 





real Ij = 2.5, -0.5, -1.5, J = -0.5 




kj = 2.28523787148992, -0.665436431311280, -1.61980144017864, A a = -1.18390417778921 


34 


0.00000000000000 


3/2 





quartet Ij = 2, 0, -2 


35 


0.00000000000000 


3/2 


1 


quartet Ij = 3, 1, — 1 


36 


0.041639350031250 


1/2 


2 


real Ij = 0.5, -1.5, -2.5, J = -0.5 




kj = 0.242333631055016, -1.61486765220468, —2.81625618363673, A a = -1.16526625596158 


37 


0.041639350031250 


1/2 


4 


veallj =2.5,1.5,-0.5, J = 0.5 




^ = 2.81625618363673, 1.61486765220468, -0.242333631055016, A a = 1.16526625596158 


38 


0.599066446144297 


1/2 


1 


real I 3 = -0.5, -1.5, -2.5, J = -0.5 




kj = -0.702579425111547, -1.67294123942593, -2.86046709144551, A« = -1.39028884293989 


39 


0.599066446144296 


1/2 


5 


real Ij = 2.5, 1.5, 0.5, J = 0.5 




kj = 2.86046709144551, 1.67294123942593, 0.702579425111547, A a = 1.39028884293989 


40 


0.625362070788712 


1/2 


4 


real I 3 = 1.5, -1.5, -2.5, J = 0.5 




kj = 1.49854259024225, -1.26094870462919, -2.33198898800626, A a = 0.722114430777064 


41 


0.625362070788713 


1/2 


2 


real Ij = 2.5, 1.5, -1.5, J = -0.5 




kj = 2.33198898800626, 1.26094870462919, —1.49854259024225, A a = -0.722114430777064 


42 


1.00000000000000 


3/2 


1 


quartet Ij = 3, 1, 


43 


1.00000000000000 


3/2 


5 


quartet Ij = 3, 0, —1 


44 


1.00000000000000 


3/2 


1 


quartet Ij — 2, 1, — 1 


45 


1.00000000000000 


3/2 


5 


quartet Ij = 1, —1, —2 



No. 


Energy 


S 


P/(tt/3) 


types 


46 


1.24466526322010 


1/2 


3 


real Ij = 1.5, -1.5, -2.5, J = -0.5 


k 3 = 1.23359323035324, -1.58513579306981, -2.79005009087323, A a = -1.05370338517746 


47 


1.24466526322010 


1/2 


3 


real/j = 2.5, 1.5, —1.5, J = 0.5 


kj = 2.79005009087323, 1.58513579306981, -1.23359323035324, A Q = 1.05370338517746 


48 


1.26651528335317 


1/2 





real Ij = 2.5, -0.5, -2.5, J = 0.5 


kj = 2.65775114820892, -0.311801389555702, -2.34594975865322, A a = 0.614983937128585 


49 


1.26651528335317 


1/2 





real Ij = 2.5, 0.5, -2.5, J = -0.5 


kj = 2.34594975865322, 0.311801389555702, -2.65775114820892, A a = -0.614983937128585 


50 


1.55774931070832 


1/2 


5 


real Ij = 2.5, -0.5, -2.5, J = -0.5 


kj = 2.31151260511968, -0.609456221165668, -2.74925393515061, A a = -0.901701357480205 


51 


1.55774931070832 


1/2 


1 


real Ij = 2.5, 0.5, -2.5, J = 0.5 


kj = 2.74925393515061, 0.609456221165668, -2.31151260511968, A a = 0.901701357480206 


52 


2.00000000000000 


3/2 





quartet Ij = 3, 2, 1 


53 


2.00000000000000 


3/2 





quartet Ij = 3, —1, —2 


54 


2.00000000000000 


3/2 


2 


quartet Ij = 3, 1, -2 


55 


2.00000000000000 


3/2 


4 


quartet Ij — 3, 2, —1 


56 


2.59087887855288 


1/2 


4 


real I 3 = 2.5, -1.5, -2.5, J = 0.5 


kj = 2.60714264000049, -1.28680359174393, -2.36753659945316, A Q = 0.468661303303178 


57 


2.59087887855288 


1/2 


2 


real Ij = 2.5, 1.5, -2.5, J = -0.5 


kj = 2.36753659945316, 1.28680359174393, -2.60714264000049, A a = -0.468661303303178 


58 


3.00000000000000 


1/2 


3 


77-pairing ®Ij = 


59 


3.24859982950560 


1/2 


4 


real Ij = 2.5, -1.5, -2.5, J = -0.5 


k 3 = 2.29420005680492, -1.59212356246732, -2.79647159673079, A a = -1.07985902893246 


60 


3.24859982950560 


1/2 


2 


real Ij = 2.5, 1.5, -2.5, J = 0.5 


k 3 = 2.79647159673079, 1.59212356246732, —2.29420005680492, A a = 1.07985902893246 
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No. 


Energy 


S 


P/(tt/3) 


types 


61 


3.53450070903252 


1/2 


4 


complex to = 10.0, £ = 6.0, f = 1.5, a = -1.0 


q = 2.18262982829371, £ = 1.51998108466021, k 3 = 6.10671585537857, A = 1.96074971503998 


Red = 0.138234567206297 x 10~ 3 , ImS = 0.236218958183487 x 10~ 3 


62 


3.53450070903252 


1/2 


2 


complex to = 8.0, 1 = 1.0, f = -1.5, a = 1.0 


q = 4.10055547888588, £ = 1.51998108466021, fc 3 = 0.176469451801019, A = -1.96074971503998 


Red = -0.138234567211404 x 10" 3 , Im<5 = 0.236218958176826 x 10" 3 


63 


4.00000000000000 


3/2 


3 


quartet Ij — 3, 2, —2 


64 


4.00000000000000 


1/2 


4 


77-pairing ®Ij = 1 


65 


4.00000000000000 


1/2 


2 


^-pairing ®Ij = — 1 


66 


4.13147838656489 


1/2 


5 


complex to = 5.0, 1 = 1.0, J = 1.5, a = —1.0 


q = 2.23719190707960, £ = 1.45368482839940, k 3 = 0.761603941823783, A = 1.77325189299191 


Re (5 = 0.307921735414718 x 10" 3 , 1m5 = 0.266758568093550 x 10" 3 


67 


4.13147838656489 


1/2 


1 


complex to = 13.0, 1 = 6.0, j' = -1.5, a = 1.0 


q = 4.04599340009998, £ = 1.45368482839940, k 3 = 5.52158136535580, A = -1.77325189299191 


RcS = -0.307921735405836 x 10~ 3 , 1m 6 = 0.266758568103764 x 10~ 3 


68 


4.33177213111493 


1/2 


5 


complex to = 11.0, £ = 6.0, J — 0.5, cr — —1.0 


q = 2.77491534511596, £ = 1.10073629130283, k 3 = 5.96934237293065, A = 0.601275696875128 


RcS = -0.273302194108371 x 10" 2 , ImS = -0.199420060157429 x 10~ 2 


69 


4.33177213111493 


1/2 


2 


complex to = 7.0, 1 = 1.0, J' = -0.5, a = 1.0 


q = 3.50826996206362, £ = 1.10073629130283, fc 3 = 0.313842934248937, A = 0.601275696875128 


Red = -0.119981837180916 x 10 1 , ImS = -0.199420060157118 x 10" 2 


70 


4.57467498439465 


1/2 





complex to = 6.0, £ = 1.0, j' = 0.5, a = 1.0 


q = 2.88935479499487, ^ = 1.07299907022393, fc 3 = 0.504475717189853, A = 0.411558289383258 


Red = -0.399317226400675 x 10" 2 , 1m5 = 0.222938603299694 x 10~ 3 



No. 


Energy 


S 


P/(tt/3) 


types 


71 


4.57467498439465 


1/2 





complex to = 12.0, £ = 6.0, f = -0.5, a = -1.0 


q = 3.39383051218472, £ = 1.07299907022393, fc 3 = 5.778709589989730, A = 0.411558289383258 


Re (5 = -0.819123406502514 x 10°, 1m 8 = 0.222938603298584 x 10~ 3 


72 


4.71197559752276 


1/2 


3 


complex to = 9.0, 1 = 5.0, J = 1.5, a = —1.0 


q = 2.16082673394224, 1.54894763250639, k 3 = 5.10312449288489, A = 2.04356190829635 


Re 5 = 0.892857558656424 x 10~ 4 , ImS = 0.211992582481946 x 10~ 3 


73 


4.71197559752276 


1/2 


3 


complex to = 9.0, 1 = 2.0, f = -1.5, a = 1.0 


q = 4.12235857323734, £ = 1.54894763250639, k 3 = 1.18006081429470, A = -2.04356190829635 


Re ,5 = -0.892857558563165 x 10" 4 , Im<5 = 0.211992582494824 x 10" 3 


74 


5.58600548103459 


1/2 


4 


complex to = 10.0, 1 = 5.0, f = 0.5, a = -1.0 


q = 2.72669142318074, £ = 1.11641164129461, k 3 = 5.01859266560450, A = 0.683372926842314 


Rc5 = -0.186687577849687 x 10" 2 , ha 5 = -0.244844007284084 x 10~ 2 


75 


5.58600548103459 


1/2 


2 


complex to = 8.0, £ = 2.0, f = -0.5, a = 1.0 


q = 3.55649388399885, £ = 1.11641164129461, k 3 = 1.26459264157509, A = 0.683372926842313 


ReS = -0.136487897790614 x 10 1 , Im5 = -0.244844007284262 x 10" 2 


76 


5.95823319001949 


1/2 





complex to = 6.0, £ = 2.0, J = 1.5, a = —1.0 


q = 2.27376322866332, ^= 1.41368501399761, k 3 = 1.73565884985295, A = 1.66056016294642 


RcS = 0.455885386176691 x 10" 3 , ha 8 = 0.245749596857303 x 10" 3 


77 


5.95823319001949 


1/2 





complex to = 12.0, £ = 5.0, f = -1.5, a = 1.0 


q = 4.00942207851627, £ = 1.41368501399761, fc 3 = 4.54752645732663, A = -1.660560162946420 


RcS = -0.455885386181354 x 10" 3 , ha 8 = 0.245749596852640 x 10" 3 


78 


6.00000000000000 


1/2 


5 


ry-pairing ®Ij = 2 


79 


6.00000000000000 


1/2 


1 


?7-pairing ®/j = — 2 


80 


6.03334376215914 


1/2 


1 


complex to = 7.0, £ = 2.0, j' = 0.5, a = 1.0 


q = 2.96456730023144, £ = 1.06124089913747, k 3 = 1.40124825791330, A = 0.288686041886590 


RcS = -0.375435197170576 x 10" 2 , ha 6 = 0.208597324485038 x 10" 2 



70 



No. 


Energy 


S 


P/(tt/3) 


types 


81 


6.03334376215914 


1/2 


5 


complex to = 11.0, £ = 5.0, f = -0.5, a = -1.0 


q = 3.31861800694814, £ = 1.06124089913747, k 3 = 4.88193704926629, A = 0.288686041886588 


Re<5 = -0.573617731801462, ImS = 0.208597324485194 x 10~ 2 


82 


6.70841301401077 


1/2 


2 


complex to = 8.0, 1 = 4.0, J = 1.5, a = -1.0 


q = 2.16393969303817, £ = 1.54471884698899, k 3 = 4.04970102349645, A = 2.03142493862260 


Red = 0.956259360065381 x 10~ 4 , ImS = 0.215691965483877 x 10~ 3 


83 


6.70841301401077 


1/2 


4 


complex to = 10.0, 1 = 3.0, f = -1.5, a = 1.0 


q = 4.11924561414142, £ = 1.54471884698899, k 3 = 2.23348428368314, A = -2.031424938622590 


Red = -0.956259360207490 x 10" 4 , ImS = 0.215691965477882 x 10" 3 


84 


7.00000000000000 


1/2 





77-pairing (g>Ij = 3 


85 


7.48332665113181 


1/2 


1 


complex to = 7.0, 1 = 3.0, j' = 1.5, a = -1.0 


= 2.20081321450656, £ = 1.49694330556893, k 3 = 2.92875642936307, A = 1.89535079698079 


Red = 0.187321491183834 x 10~ 3 , ImS = 0.252337664304214 x 10~ 3 


86 


7.48332665113181 


1/2 


5 


complex to = 11.0, £ = 4.0, f = -1.5, a = 1.0 


q = 4.08237209267303, £ = 1.49694330556893, k 3 = 3.354428877816510, A = -1.89535079698079 


Red = -0.187321491188275 x 10" 3 , Im<5 = 0.252337664298441 x 10" 3 


87 


7.59363119464461 


1/2 


3 


complex to = 9.0, 1 = 4.0, f = 0.5, a = -1.0 


q = 2.74080575567842, £ = 1.11156342668968, k 3 = 3.94316644941255, A = 0.659158265098385 


Red = -0.212814421543628 x 10" 2 , Imd = -0.234942161399343 x 10" 2 


88 


7.59363119464461 


1/2 


3 


complex to = 9.0, i = 3.0, f = -0.5, a = 1.0 


q = 3.54237955150117, £ = 1.11156342668968, k 3 = 2.34001885776704, A = 0.659158265098384 


Re,5 = -0.131618838598134 x 10 1 , Im5 = -0.234942161399543 x 10~ 2 


89 


8.02701965828632 


1/2 


2 


complex to = 8.0, 1 = 3.0, f = 0.5, a = 1.0 


q = 2.89722498813356, £ = 1.07154335884761, k 3 = 2.58313043330567, A = 0.398665945793610 


Red = -0.401340764685176 x 10" 2 , Im<5 = 0.414789194807641 x 10" 3 


90 


8.02701965828632 


1/2 


4 


complex to = 10.0, 1 = 4.0, / = -0.5, a = -1.0 


q = 3.38596031904603, £ = 1.07154335884761, k 3 = 3.70005487387392, A = 0.398665945793609 


Red = -0.793318483940372, ImS = 0.414789194806753 x 10~ 3 
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